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ABSTRACT 


Ihe  experimental  properties  of  tin  and  lead  variable- 
thickness microbridges  have  been  investigated.  We  find  that 
the  dc  I-V  characteristics  and  the  high  frequency  response 
of  these  Josephson  devices  are  determined  or  modified  by 
nonequilibrium  effects. 

We  find  that  the  improved  cooling  of  the  variable- 
thickness geometry  allows  us  to  observe  microwave-induced 
Josephson  steps  to  voltages  as  high  as  4,2  mV,  a factor  of  4 
improvement  over  conventional  thin-film  microbridges,  This 
improved  performance  is  analyzed  in  terms  of  a simple  heat- 
ing approximation  for  three-dimensional  weak  links.  From 

this  analysis  we  find  that  the  critical  current  decreases 

-P/P 

with  the  dissipated  power  P as  e ° where  P^  is  a material- 
and  temperature-dependent  parameter.  This  prediction  of  the 
heating  model  is  found  to  be  in  very  good  agreement  with  the 
reduction  of  the  amplitude  of  the  microwave-induced  ac  steps 
at  all  voltages  as  well  as  the  temperature  dependence  of  the 
maximum  voltage  at  which  these  steps  are  observed. 

We  use  the  gap  and  the  subharmonic  gap  structure  observed 
in  the  dc  I-V  characteristics  of  our  microbridges  as  a local 
thermometer  in  order  to  demonstrate  the  correlation  between 
low  dissipation,  reduced  hysteresis,  and  the  good  high 
frequency  performance  of  our  variable -thickness  microbridges. 


i 


ii 


Measurement  of  low  voltage  foot-like  features  of  the 
I-V  characteristics  are  found  to  be  consistent  with  a 
voltage-dependent  enhancement  of  the  supercurrent  induced  by 
the  disequilibrium  of  the  quasi-particles  during  the 
Josephson  cycle.  Comparison  of  the  experimental  data  with 
the  predictions  of  the  models  of  Golub  and  of  Aslamazov  and 
Larkin  show  many  points  of  agreement  in  the  regimes  where 
these  theories  are  applicable. 

We  also  report  on  the  general  characteristics  of  these 
devices  such  as  the  magnitude  and  temperature  dependence  on 
the  critical  currents,  products,  and  the  Wyatt-Dayem 

effect.  Exploratory  experiments  on  series  arrays  of  micro- 
bridges and  detection  of  far-infrared  radiation  are  also 


discussed 
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CHAPTER  ONE 


INTRODUCTION 

Superconducting  microbridges,  in  which  a superconducting 
film  is  narrowed  down  to  a small  constriction,  are  ideal  for 
the  study  of  a large  variety  of  physical  phenomena.  In  very 
long,  narrow  microbridges  (length  L larger  than  a few  urn, 
width  w less  than  5 urn),  as  shown  in  Fig.  1.1  (a),  one  has  an 
essentially  one-dimensional  system  in  which  the  superconducting 
order  param.eter  varies  only  along  the  length  of  the  bridge. 

Thus  such  bridges  provide  a simple  system  for  the  study  of 
time-dependent,  dissipative,  nonequilibrium  processes  in 
superconductors.  For  short  (L  < 1 iim)  and  narrow  (w  < 1 urn) 
microbridges,  in  either  a two-dimensional  (Fig.  1.1  (b) ) or 
variable-thickness  geometry  (Fig.  1.1  (c) ) , one  can  study 
the  conditions  for  the  existence  of  ideal  Josephson  behaviour 
and  how  this  behaviour  is  modified  or  limited  by  nonequilibrium 
processes . 

Despite  a considerable  amount  of  work,  the  understand- 
ing of  these  weak  superconducting  systems  has  progressed 
slowly  mainly  due  to  the  complex  behaviour  of  nonequilibrium 
superconductors.  In  the  case  of  short  microbridges,  of 
interest  in  this  work,  it  became  clear  in  the  earlier  work 
of  Dayem  and  Anderson^  that  many  of  the  properties  of  bridges 
could  be  attributed  to  some  extent  to  Josephson  behaviour. 


1 


3 


2 3 

Subsequent  improvement  in  fabrication  techniques  ' , which 

reduced  the  size  of  the  microbridges  by  almost  an  order  of 
magnitude,  showed  that  at  least  near  the  properties  of 
microbridges  were  in  qualitative  and  to  some  extent  quanti- 
tative agreement  with  ideal  Josephson  behaviour. 

4 

It  was  not  until  1973  that  Skocpol,  Beasley  and  Tinkham 
showed  that  Joule  dissipation  played  a major  role  in  determin- 
ing the  properties  of  all  constriction  weak  links.  In  par- 
ticular,  this  idea  explained  the  pronounced  hysteresis  ob- 
served in  the  I-V  characteristics  of  planar  microbridges 
and  the  fact  that  the  ac  Josephson  effect  was  only  observed 
up  to  voltages  of  the  order  of  0.5  mV,  corresponding  to 
frequencies  of  500  GHz.  This  voltage  limit  is  significantly 
lower  than  that  of  superconducting  point  contacts,  in  which 
ac  effects  are  observed  up  to  voltages  of  the  order  of  15  mV. 

It  was  suggested  that  the  superiority  of  the  point  contact 
geometry  arises  from  its  three  dimensionality,  which  more 
efficiently  removes  the  heat  generated  in  the  contact.  Fur- 
thermore, point  contacts  have  higher  resistances,  which 

implies  lower  dissipation  at  the  same  voltage  since  the  dis- 

2 

sipation  scales  as  V /R^j.  Thus,  it  was  recognized  that  a 
reduced-heating  or  variable-thickness  geometry,  in  which  a 
thin  bridge  is  connected  to  much  thicker  Ijanks,  would  signi- 
ficantly reduce  heating  effects  due  to  the  improved  cooling 
provided  by  the  thicker  banks.  At  the  same  time,  Likharev^ 
and  coworkers  in  the  Soviet  Union  considered  var iable-thickness 
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microbridges,  not  only  because  of  their  better  cooling,  but 
also  because  the  lower  current  density  in  the  banks  should 
restrict  nonlinearities  to  the  bridge  region,  so  that  the 
microbridges  should  behave  more  like  ideal  Josephson  junc- 
tions . 

At  that  point  we  began  to  try  to  make  lead  and  tin 
variable-thickness  superconducting  microbridges.  The  advan- 
tage of  a thin  film  geometry  for  high  frequency  applications 
is  that  it  is  a permanent  structure  which  is  compatible  with 
thin-film  techniques  for  the  fabrication  of  antennas,  arrays, 
and  other  complex  systems.  It  was  our  goal  to  fabricate 
variable-thickness  bridges  in  which  heating  effects  would  be 
significantly  reduced  and,  if  possible,  in  which  the  ac 
Josephson  effect  could  be  observed  at  far-infrared  fre- 
quencies . 

Once  a reliable  fabrication  technique  was  developed,  as 
discussed  in  Chapter  Two,  it  was  found  that  the  properties 
of  our  microbridges  were  remarkably  different  from  those 
of  two-dimensional  microbridges.  The  most  significant  dif- 
ference was  the  observation  of  the  ac  effect  up  to  voltages 
a factor  of  four  higher  than  had  previously  been  observed 
in  any  thin-film  weak  link.  This  implied  the  presence  of 
ac  supercurrents  flowing  through  the  weak  links  at  far- 
infrared  frequencies  and  suggested  the  possibility  of  using 
our  variable-thickness  microbridges  for  the  detection  of 
far-infrared  radiation. 
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In  Figure  1.2,  the  I-V  characteristics  of  two  of  our 
tin  microbridges  are  shown.  The  bridge  in  Fig.  1.2  (a) 
corresponds  to  a well-cooled  bridge  in  which  the  ac  Joseph- 
son  effect  was  observed  up  to  2,3  mV.  Note  that  at  the  low- 
est temperatures  a multivalued  curve  develops.  (The  seg- 
ments of  negative  slope  are  the  result  of  averaging  over 
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circuit-controlled  relaxation  oscillations  ).  Furthermore, 
very  sharp  "bumps"  corresponding  to  the  gap  and  subharmonic 
gap  structures,  discussed  in  Chapter  Five,  are  observed. 

While  the  performance  of  this  bridge  is  already  a considerable 
improvement  over  planar  microbridges  it  should  be  contrasted 
with  the  I-V  characteristics  shown  in  Fig.  1.2  (b) . This 
corresponds  to  a bridge  of  resistance  similar  to  the  pre- 
vious one,  but  which  is  better  cooled  because  of  much 
thicker  banks.  In  this  case,  the  I-V  characteristics  are 
very  smooth,  with  no  multivaluedness  even  at  the  lowest  tem- 
peratures. When  microwave  radiation  is  applied  to  this 
bridge,  ac  Josephson  steps  are  observed  up  to  3.7  mV.  That 
this  improved  behaviour  is  due  to  better  cooling  is  confirmed 
by  the  presence  of  the  gap  structure  at  much  higher  voltages 
corresponding  to  lower  temperatures  than  in  the  first  bridge. 
As  discussed  in  Chapter  Five,  we  have  used  this  gap  struc- 
ture as  a local  thermometer  of  the  average  temperature  rise 
in  the  bridge  region  in  order  to  quantitatively  demonstrate 
the  connection  between  high  frequency  performance,  hysteresis, 
and  heating. 
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Because  of  the  novel  nature  of  our  microbridges,  we 
have  attempted  to  characterize  their  overall  properties.  In 
the  course  of  our  work,  we  accumulated  a large  amount  of 
data  on  the  critical  currents  of  bridges  with  a wide  distri- 
bution of  resistances.  It  was  found  that  our  results  dif- 
fered quantitatively  from  theory  both  near  and  at  low 
temperatures,  as  discussed  in  Section  3.1.  On  the  other 
hand,  the  temperature  dependence  of  the  critical  current 
was  found  to  be  Josephson-like  (linear)  for  both  tin  and 

lead  bridges  near  T^,  despite  their  differences  in  coherence  I 

length.  , 

Due  to  the  reduced  hysteresis,  some  of  the  features  of 
the  I-V  characteristics  became  very  pronounced  and  could 
be  followed  over  a much  wider  temperature  range  than  had 
previously  been  possible.  We  became  very  interested  in 
the  foot-like  features  at  low  voltages  which  can  be  seen  in 
the  I-V  characteristics  of  both  bridges  shown  in  Fig.  1.2. 

Note  that  in  the  case  of  the  better-cooled  bridge  this 
feature  can  be  observed  all  the  way  down  to  2.1  K,  because 
of  the  absence  of  hysteresis.  As  discussed  in  Chapter  Four, 
these  features  are  found  to  be  consistent  with  a voltage- 
dependent  enhancement  of  the  supercurrent  induced  by  the  dis- 
equilibrium of  the  quasiparticles  during  the  Josephson  cycle. 

Comparison  of  our  experimental  results  with  the  predictions 

. 7 . 8 

of  the  theories  of  Golub  and  of  Aslamasov  and  Larkin  yield 


good  agreement  with  many  aspects  of  our  data,  although  some 
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questions  remain  open. 

Earlier  indications  were  that  the  voltage  limit  for  the 
ac  Josephson  effect  scaled  with  the  superconducting  energy 
gap  2A/e,  which  favors  high  materials."^  However,  our 
experimental  results  on  lead  microbridges  differed  little 
from  those  on  tin,  suggesting  that  a more  careful  analysis 
of  the  three-dimensional  problem  was  required.  This,  to- 
gether with  our  experimental  observations  indicating  that  at 
high  dissipation  levels  the  basic  Bessel- function-like  peri- 
odicity of  the  amplitudes  of  the  ac  steps  remained,  led  us 
to  a simple  heating  model  for  three  dimensional  geometries, 
discussed  in  Chapter  Five.  This  model  predicts  that  the 

critical  current  decreases  exponentially  with  power  as 
-P/P 

e where  P is  the  power  dissipated  in  the  junction  and 

P^  depends  on  temperature,  material  parameters,  and  junction 
geometry.  From  this  model  it  is  possible  to  infer  the  best 
possible  materials  and  geometries  for  good  high-frequency 
performance.  The  predictions  of  the  model  are  in  very  good 
agreement  with  our  results  in  variable-thickness  micro- 
bridges as  well  as  with  other  results  found  in  the  literature 
for  both  variable- thickness  microbridges  and  point  contacts. 

The  observation  of  ac  steps  at  voltages  corresponding 
to  far-infrared  frequencies  encouraged  us  to  attempt  the 
observation  of  an  ac  Josephson  step  induced  by  a far-infrared 
laser.  This  unsuccessful  effort  is  described  in  Section  6.1. 
The  main  difficulty  was  our  inability  to  couple  enough  far- 
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infrared  radiation  to  our  low  impedance  microbridges  (R^  < 

1^^)  . While  a very  small  decrease  of  the  critical  current 
was  observed,  it  was  not  possible  to  determine  whether  this 
was  due  to  heating  or  to  the  Josephson  effect.  While  it 
might  be  possible  to  increase  the  impedance  of  the  micro- 
bridges by  proper  choice  of  materials  and  design,  another 
possibility  is  to  use  arrays  of  microbridges  in  series.  In 
Section  6.2,  we  include  some  of  our  experimental  results 
with  series  pairs  of  microbridges.  While  no  coupling  effects 
were  observed  at  high  voltages,  some  interesting  effects 
were  observed  at  low  voltages,  indicating  that  further  study 
of  these  systems  is  needed. 

The  organization  of  this  report  is  as  follows;  In 
Chapter  Two  the  experimental  techniques  used  in  the  fabrica- 
tion and  measurement  of  the  properties  of  our  bridges  are 
presented.  Chapter  Three  discusses  the  general  character- 
istics of  our  microbridges,  such  as  critical  currents  and 
the  Wyatt-Dayem  effect.  Chapter  Four  deals  with  the  theory 
and  experimental  results  on  the  nonequilibrium-enhanced  super- 
currents at  low  voltages.  In  Chapter  Five  the  heating  model 
is  described  and  comparison  is  made  to  our  experimental 
results  both  at  high  and  low  voltages.  Chapter  Six  con- 
tains some  discussion  of  our  exploratory  experiments  with 
far-infrared  radiation  as  well  as  with  arrays  of  microbridges. 
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CHAPTER  TWO 


EXPERIMENTAL  TECHNIQUES 

2 . 1 Sample  Preparation 

Since  our  primary  goal  was  the  study  of  the  physical 
properties  of  variable-thickness  microbridges,  we  decided 
upon  a sin.ple  and  quick  microfabrication  technique  which 
allowed  us  to  fabricate  a complete  sample  using  only  dia- 
mond knives  and  microneedles.  While  this  technique  has 
drawbacks  in  terms  of  reproducibility  and  size  control, 
once  a film  was  evaporated  a good  quality  sample  could  be 
made  in  less  than  half  an  hour.  Furthermore,  this  technique 
allowed  us  to  make  very  thick  banks  without  the  limitations 
of  more  involved  techniques  using  resists,  in  which  the 
maximum  thickness  is  limited  by  the  resist  thickness. 

Because  of  the  stresses  associated  with  cutting  the 
films  with  a diamond  knife,  it  was  found  essential  that  the 
films  adhere  well  to  the  substrate.  To  achieve  this,  we  use 
a technique  similar  to  that  described  by  Chiao  et  al^^.  We 
first  clean  the  sapphire  substrates  ultrasonically  in  suc- 
cessive solutions  of  Sparkleen  detergent,  distilled  water 
and  isopropyl  alcohol.  The  substrates  are  then  wiped  clean 
with  Kodak  lens  paper  and  are  placed  on  a hot  plate  where 
small  indium  patches  are  attached  to  them.  Two  three-inch 
solder -coated  Nb^Sn  tapes  are  then  attached  to  the  indium 
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patches  in  order  to  provide  electrical  contact  to  the  films 
which  will  be  evaporated  on  top  of  the  substrates.  The  sub- 
strates arc  then  placed  on  a copper  plate  which  is  screwed 
to  a copper  block  in  the  evaporator.  This  copper  block  can  be 
cooled  to  liquid  nitrogen  temperatures. 

Once  in  the  evaporator,  the  substrates  are  cleaned  in 
an  oxygen  glow  discharge  for  periods  of  half  an  hour  to  an 
hour.  The  glow  discharge  is  produced  by  a pure  aluminum 
wire  at  a voltage  of  600  volts  and  currents  of  25-30  mA 
corresponding  to  pressures  of  0.03  Torr.  Care  is  taken  that 
there  is  no  dixect  line  of  sight  between  the  aluminum  wire 
and  the  substrates  as  this  yields  buckled  films,  presumably 
because  of  the  sputtering  of  aluminum  atoms  onto  the  sub- 
strates. After  this  step,  a thin  tin  or  lead  undercoat 
(50-100  A thick  depending  on  the  ultimate  thickness  of  the 
film  to  be  evaporated)  is  evaporated  and  allowed  to  oxidize 
in  the  glow  discharge  for  periods  of  two  to  four  hours. 

After  the  samples  are  cooled  to  liquid  nitrogen  temperatures, 
a thick  1-3  am  film  is  evaporated  at  pressures  of  2 x 10~ 

Torr  and  rates  of  90-200  ’/sec.  The  thickness  of  the  film 
is  monitored  with  a Sloan  200  deposition  thickness  monitor. 

Using  the  above  technique  we  obtain  very  well  adhered 
films  which  withstand  the  'tape  test",  in  which  a piece  of 
Scotch-Brand  tape  is  attached  to  the  film  and  pulled  without 
any  of  the  film  coming  off  the  substrate.  The  good  adhesion 
also  offers  two  additional  desirable  features.  First,  the 
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leads  can  be  attached  to  the  films  after  the  evaporation. 

This  is  very  convenient  if  a large  number  of  samples  are 
fabricated  or  if  a lead  accidentally  comes  off  the  substrate 
after  the  microbridge  has  been  fabricated.  Second,  the  good 
adhesion  significantly  reduces  the  burnout  rate  of  our 
samples.  In  those  cases  in  which  a bridge  with  good  adhesion 
was  accidentally  destroyed,  the  burnout  region  was  restricted 
to  the  bridge  region  only.  In  the  case  of  a film  with  poor 
adhesion  this  burnout  region  was  found  to  spread  into  the 
banks  for  several  urn. 

In  attempting  to  reduce  the  burnout  rate  of  our  bridges, 
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we  also  fabricated  our  microbridges  on  germanium  substrates 
to  provide  a high  temperature  short  to  the  bridge.  While 
this  was  found  to  be  helpful,  the  germanium  substrates  were 
found  to  be  inconvenient  to  handle  as  they  were  brittle  and 
had  a tendency  to  break  along  any  cut  made  on  the  substrate 
with  a diamond  knife.  Another  disadvantage  was  that  while 
our  sapphire  substrates  could  be  reused  repeatedly,  the 
germanium  substrates  could  not.  Cleaning  the  films  off  the 
germanium  substrates  produced  an  oxide  layer  on  the  surface 
of  the  substrate  and  thus  there  was  no  electrical  contact 
between  the  substrate  and  the  film  evaporated  on  it. 

Two  types  of  diamond  knives,  both  manufactured  by  the 
J.  Robert  Moore  Co.,  were  used  in  cutting  our  films.  Ihe 
first  type  was  a sharp  point  knife  formed  by  the  intersection 
of  three  cleaved  pianos  which  was  used  for  making  small 
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grooves  in  the  substrates,  Tlie  second  typo  was  a long  (2  rrm.) 
straight-edge  knife  which  was  used  for  the  removal  of  largo 
areas  of  film  and  in  the  fabrication  of  the  microbridges. 

The  diamond  knives  are  mounted  at  the  end  of  a rotat- 
able rod  attached  to  a micromanipulator.  This  micromanipu- 
lator provides  coarse  and  fine  control  in  the  vortical  di- 
rection as  well  as  angular  control.  Thus  the  edge  of  the 
knife  can  be  aligned  parallel  to  the  surface  of  the  film  to 
be  cut.  Tfie  substrates  are  placed  on  a differential  screw 
translation  stage  which  provides  translational  motion  in  the 
X-Y  directions.  The  knife  and  the  area  to  be  cut  can  be 
viewed  through  a Bausch  and  Lomb  StereoZoom  microscope  with 
a maximum  power  of  140X, 

The  grooves  on  the  sapphire  substrates  are  made  by 
sliding  the  substrates  beneath  the  single  point  knife  while 
periodically  adjusting  the  height  of  the  knife  to  ensure 
that  the  smallest  possible  groove  is  made.  In  this  manner 
we  obtained  grooves  typically  0.2-0. 3 pm  wide  and  0.1-0. 3 
pm  deep.  Typically  5-10  of  these  grooves  were  made  on  each 
substrate  so  that  a different  groove  could  be  used  each  time 
a substrate  was  reused.  This  was  done  to  avoid  possible 
shorts  caused  by  scratches  on  the  sapphire  substrate  over 
areas  which  were  previously  used. 

After  the  films  are  evaporated  onto  the  grooved  substrates, 
the  flat  long  knife  is  used  to  remove  the  film  from  several 
small  rectangular  regions  on  each  side  of  the  groove  as  shown 
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in  Fig.  2.1  (a)  and  (b)  . TlTon  the  flat  knife  is  aligned 
parallel  to  the  film  surface  by  rotating  the  knife  edge 
until  it  is  parallel  to  its  own  reflection  on  the  surface 
of  tlie  shiny  film.  The  long  flat  edge  is  then  gently 
pressed  against  the  film,  leaving  a small  thin  bridge  in  the 
groove,  but  separating  the  thick  banks  attached  to  it  as 
shown  in  Fig.  2.1  (c)  and  in  the  SEM  pictures  in  Fig.  2.2. 

Tlie  essential  difference  from  previous  "double-scratch" 
techniques  is  that  the  whole  film  is  cut  simultaneously  by 
a single  edge  rather  than  by  a point  dragged  through  the 
film.  This  has  the  advantage  that  a simple  manual  scheme  can 
be  used  since  the  sapphire  substrates  stop  the  diamond  knife 
from  penetrating  through  the  bridge.  (The  knife  actually 
penetrates  a little  in  order  to  separate  tlie  banks  but  this 
is  usually  not  enough  to  cut  completely  through  the  bridge 
in  tlie  groove.)  Using  this  fabrication  technique,  an  experi- 
enced operator  can  obtain  an  80-90%  yield;  the  main  problem 
is  not  in  obtaining  a microbridge,  but  in  aligning  the 
knife  edge  so  that  the  banks  are  completely  separated  in 
regions  far  from  the  bridge. 

After  15-30  bridges  are  made  along  a single  groove, 
they  arc  inspected  with  an  optical  microscope  at  400X  and 
the  best  bridge  is  selected  for  testing.  Best  is  defined 
as  the  optimum  combination  of  a small  bridge,  located  where 
the  groove  is  small,  with  the  banks  completely  separated. 

A Circon  microneedle  is  then  used  to  isolate  the  selected 
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Figure  2.2.  Scanning  electron  micrographs  of  a bridge. 
Top:  Overview  of  bridge  and  banks. 

Bottom;  Close-up  of  the  microbridge. 
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bridge  from  others  on  the  substrate  by  removing  the  film 
from  lines  connecting  the  small  rectangular  openings  on 
each  side  of  the  bridge  to  the  edges  of  the  substrate,  A 
small  short  is  always  left  near  the  leads  to  prevent  acci- 
dental burnout.  This  short  is  only  removed  once  the  leads 
are  soldered  (with  an  unplugged  and  grounded  soldering  iron) 
to  the  cryostat  and  all  electronics  are  turned  on.  If  this 
procedure  is  followed,  burnout  is  rare  as  long  as  all  the 
equipment  is  left  on.  Once  the  bridge  is  cooled  down  to 
helium  temperatures,  the  bridge  becomes  less  susceptible  to 
burnout  even  if  the  electronics  are  turned  off. 

Ihe  same  techniques  outlined  above  have  been  used  in 
the  fabrication  of  arrays  of  microbridges  in  series.  In 
this  case  we  found  it  necessary  to  keep  the  thickness  of  the 
banks  less  than  1 urn  or  the  fabrication  of  the  second  micro- 
bridge would  push  enough  material  onto  its  neighbor  either 
to  distort  it  or  to  completely  short  the  banks.  The  physi- 
cal similarity  between  adjacent  bridges  was  found  to  be 
primarily  determined  by  the  uniformity  of  the  groove?  if 
sufficient  care  was  taken,  a single  groove  could  be  made  uni 
form  over  lengths  of  the  order  of  100-200  um. 

2.2  Experimental  Details 

The  basic  circuit  used  in  the  measurement  of  the  I-V 
characteristics  of  the  microbridges  is  shown  in  Fig,  2.3  (a) 
The  microbridges  are  essentially  "voltage  biased"  using  a 
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0.1  ohm  resistor  shunting  the  series  combination  of  the 
bridge  and  a 0,1  ohm  current-sensing  resistor.  This  net- 
work of  resistors  is  inside  the  dewar,  submerged  in  liquid 
helium.  All  the  connecting  leads  are  Nb^Sn  tapes  which 
become  superconducting  at  17  K.  The  voltage  applied  to  the 
bridge  is  controlled  by  the  current,  which  is  supplied 

to  the  network  by  an  nP61774  power  supply.  The  voltage 
across  the  bridge  is  measured  directly  with  a Princeton 
Applied  Research  Model  113  Low  Noise  Amplifier.  The  leads 
across  the  0.1  ohm  current-sensing  resistor  are  a twisted 
pair  of  #33  copper  wires.  This  pair  goes  continuously  from 
the  resistor  to  a low-thermal  connector  outside  of  the 
dewar  through  a Stycast  seal.  The  voltage  across  this 
resistor  is  chopped  by  a low-thermal,  94  Hz  mechanical  chopper 
made  by  Stevens  Arnold  Inc.  The  chopped  signal  is  trans- 
formed with  a P.A.R.  AM-1  lOOX  transformer  and  measured  with 
a P.A.R.  124  lock-in  amplifier  in  tlie  direct  mode.  Except 
for  the  current-sensing  leads,  all  other  leads  are  filtered 
with  1 mH  commercial  rf  chokes. 

For  the  derivative  measurements  we  used  the  circuit 

shown  in  Fig.  2,3  (b)  , A small  modulation  current,  I sin(u)t), 

m 

is  added  to  the  dc  current  This  produces  a signal  at 

the  modulation  frequency  which  is  proportional  to  AI/AI^ 
which  is  approximately  dl/dV.  This  signal  is  measured  with 
a P.A.R,  124  lock-in  amplifier  in  the  XlOO  transformer  mode. 

While  we  have  described  our  measuring  circuit  as  a 
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voltage-biased  circuit,  this  is  not  strictly  true.  In  our 
earlier  experiments,  the  resistance  of  our  biasing  resis- 
tors was  0.01  ohm.  This  was  sufficiently  below  the  impedance 
of  our  microbridges  to  make  the  measuring  system  very  close 
to  a voltage-biased  one,  once  the  critical  current  was 
exceeded.  Later,  in  order  to  improve  the  sensitivity  of 
the  system,  the  resistors  were  replaced  by  0.1  ohm  resistors 
which  is  of  the  order  of  the  typical  impedance  of  our  micro- 
bridges. Thus,  the  biasing  is  not  strictly  voltage  biased 
and  is  dependent  on  the  impedance  of  each  particular  bridge. 

Most  of  our  experiments  on  tin  microbridges  were  per- 
formed with  the  samples  directly  exposed  to  the  helium  bath. 

Ihe  samples  are  mounted  on  a copper  block  and  their  tempera- 
ture is  monitored  by  both  a calibrated  Ge  thermometer  and 
the  vapor  pressure  of  the  helium  bath.  The  bridges  are 
mounted  so  that  they  stick  into  the  narrow  side  of  an  X- 
band  waveguide  a distance  A.g/4  from  its  closed  end.  The 
orientation  of  the  bridges  is  such  that  they  are  parallel  to 
the  E field,  but  little  change  was  observed  if  the  orientation 
was  perpendicular  to  the  E-field.  While  this  coupling  arrange- 
ment is  quite  crude  it  allowed  us  to  explore  the  ac  response 
of  our  samples  in  detail  as  long  as  enough  microwave  power 
was  available.  For  some  of  our  lead  microbridge  experiments 
the  tip  of  the  cryostat  was  surrounded  by  a brass  vacuum 
can  and  the  temperature  regulated  electronically.  The  entire 
dewar  was  surrounded  by  Mu-metal  shields  to  reduce  the  earth's 
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magnetic  field.  The  experiments  were  performed  in  a Ray- 
Proof  RF  shielded  room,  but  little  change  was  noticed  when 
tile  door  to  the  room  was  not  closed, 

Tlie  far-infrared  experiments  were  performed  in  an  un- 
shielded dewar  with  an  optical  window  in  front  of  the  micro- 
bridge. Tlie  cryostat  used  in  these  experiments  allowed  the 
rotation  of  the  sample  for  improved  coupling  to  the  far- 
infrared  radiation. 

Three  separate  microwave  sources  were  used  in  the 
course  of  our  experiments.  Initially  we  used  a Laboratory 
for  Electronics  Ultrastable  X-band  oscillator  v/ith  a maximum 
output  power  of  the  order  of  20  mW,  Sometimes,  depending 
on  the  sample,  this  was  not  sufficient  to  driv'e  the  micro- 
bridge completely  normal,  so  we  changed  to  a fia^ed  frequency 
(10.525  GHz)  Plessey  Semiconductor  Gunn  diode  oscillator 
with  a maximum  output  power  of  300  mW,  In  the  final  stages 
of  our  work  we  used  a Varian  Associates  32  GHz  Gunn  diode 
with  an  output  power  of  60  mW.  T!ie  32  GHz  radiation  was 
brought  to  the  sample  through  the  X-band  waveguide  so  that 
the  exact  mode  at  the  sample  was  not  known. 


CHAPTER  THREE 

STATIC  PROPERTIES  OF  VARIABLE-THICKNESS  SUPERCONDUCTING 

MICROBRIDGES:  CRIITCAL  CURRENTS  AND  I R„  PRODUCTS 

CO  N 

While  the  original  theoretical  predictions  of  B.  D. 

13 

Josephson  were  made  for  an  oxide  barrier  separating  two 
superconducting  films,  experimental  investigations  have  shown 
that,  in  some  general  sense,  Josephson-like  behavior  can  be 
observed  in  any  system  in  which  two  superconductors  are 
weakly  connected.  This  weakness  can  be  introduced  by  geo- 
metrically constricting  a superconducting  film,  as  in  the 

microbridges  of  interest  here,  by  overlaying  a superconduct- 

14 

ing  film  with  a normal  metal,  by  optical,  phonon  or  quasi- 
particle injection, by  ion  implantation^^  or  even  by 

locally  reducing  the  order  parameter  with  a small  magnetic 
17 

field.  While  all  of  these  systems  exhibit  some  form  of 
Josephson  behaviour,  their  properties  can  differ  markedly 
from  ideal  behaviour  either  because  they  are  large  compared 
to  the  coherence  length  or  because  they  do  not  have  three- 
dimensional  cooling  geometries  as  do  variable-thickness  micro- 
br idges . 

In  this  chapter,  we  characterize  and  discuss  the  basic 
properties  of  our  variable-thickness  microbridges  in  the 
context  of  a simple  Ginzburg-Landau  model  of  Josephson 
behaviour  applicable  to  short  weak  links.  The  observed 
critical  currents  near  T^  are  compared  to  the  predictions  of 
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this  model,  and  it  is  found,  for  both  low  and  high  resis- 
tance microbridges,  that  the  temperature  dependence  of  the 

critical  currents  agree  qualitatively  but  not  quantitatively  ! 

with  the  predictions  of  this  model.  At  low  temperatures  the 

i 

critical  currents  are  found  to  be  comparable  to  those  expec- 
ted from  microscopic  theory.  Finally,  some  aspects  of  the 
enhancement  of  the  critical  current  in  the  presence  of 
microwaves  are  discussed. 

3.1  The  Critical  Current 

3.1.1  Regime  Close  to 

In  the  zero  voltage  state  one  deals  with  the  equilibrium, 

/ 

time-independent  state  ahd,  at  least  near  T^,  it  is  possible 
to  use  tlie  two-fluid  picture  of  the  Ginzburg-Landau  theory. 

The  superconductor  is  then  represented  by  a complex  order 
parameter  V (r)  = where  cp  is  a phase  factor  which 

is  a function  of  both  position  and  time.  The  free  energy 
density  of  a superconducting  system  is  then  represented  by  a 
pov/er  series  expansion  in  | Y | and  the  gradient  |vY1  . The 
problem  is  then  reduced  to  minimizing  the  free  energy  in  a 
manner  consistent  with  the  boundary  conditions  of  the  system 
under  consideration.  The  first  application  of  this  concept 

to  the  properties  of  small  weak  links  was  made  by  Aslamazov  j 

18 

and  Larkin.  They  took  advantage  of  the  fact  that  if  the 


weak  link  is  small  compared  to  the  temperature-dependent 
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coherence  length  I (T) , then  the  gradient  term  will  dominate 

the  free  energy  density  inside  the  weak  linko  This 

2 

implies  that  ¥ satisfies  Laplace's  equation  v Y = 0 with 

appropriate  boundary  conditions.  In  the  case  of  the 

variable-thickness  or  point-contact  geometries,  the 

boundary  conditions  are  that  the  order  parameter  equals 

that  of  the  bulk  superconductor  at  either  side  of  the 

, icp  2 

weak  link,  Y e and  Y e where  cp,  ^ are  the  phases 
o o ^ 

at  each  side  of  the  weak  link.  Aslamazov  and  Larkin  con- 
sidered a particular  solution  to  Laplace's  equation  f(r) 
which  tends  asymptotically  to  unity  on  one  side  of  the 
weak  link  and  towards  zero  on  the  other.  The  order  para- 
meter is  then  given  by 


Y (?) 


To  jfff) 


1 


f (?) 


(3.1) 


where  Y^  is  the  equilibrium  value  of  the  order  parameter 
deep  in  the  superconductors  and  the  phases  are  inde- 

pendent  of  coordinates.  The  supercurrent  density  is  given 
by 
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-3 

s 


(¥*vY  - Y VY*)  = ^ |Y  1^  sin  cp  (3.2) 
2m’^i  m*  o 


where  is  the  total  phase  difference  across  the  weak  link, 
and  vf  is  averaged  over  the  cross-section  at  the  middle  of 

18 

tlie  bridge.  This  equation  is  equivalent  to  the  Josephson 
current-phase  relation  derived  from  microscopic  theory  for  a 
superconductor-oxide-supcrconductor  tunneling  structure  with 


CO 


A e*  ^ 
m* 


vf 


(3.3) 


where  I is  the  critical  current  of  the  weak  link  and  A is 

CO 

the  cross-sectional  area  of  the  weak  link. 

Hius  for  sufficiently  small  weak  links  one  expects  an 

ideal  sinusoidal  current-phase  relation.  Such  behaviour  has 

been  observed  in  niobium  point  contacts,  small  two- 

dimensional  microbridges  and  variable-thickness  micro- 
19  20  21 

bridges.  ' ' While  we  have  neglected  linear  terms  in 

the  free  energy  density  under  the  assumption  that  the  length 
of  the  link  L is  smaller  than  the  temperature -dependent  co- 
herence Icnght  I (T) , more  detailed  calculations  by  Likharev 
22 

and  Yakobson  have  shown  that  tlie  current-phase  relation 


will  be  essentially  sinusoidal  up  to  a critical  length 
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~ 3.5  5 (T) . As  the  length  L of  the  link  becomes  larger 
than  P (T) , the  current-phase  relation  becomes  skewed  until 
it  becomes  reentrant;  and  eventually  the  solution  to  the 
Ginzburg-Landau  equation  for  tlie  order  parameter  corresponds 
to  that  of  a long  bulk  filament.  One  of  the  advantages  of 
the  variable-thickness  geometry  becomes  apparent  from  this, 
because  the  thick  banks  have  a much  smaller  current  density 
than  the  bridge,  restricting  all  nonlinearities  to  the 
bridge  region.  Thus,  if  one  has  two  bridges  with  a well 
defined  length  L,  one  with  banks  with  the  same  thickness  as 
the  bridge  and  another  with  thicker  banks,  the  variable- 
thickness bridge  will  behave  more  like  an  ideal  Josephson 
element  than  will  the  two-dimensional  bridge.  Furthermore, 
the  ideal  behaviour  will  be  preserved  over  a wider  tempera- 
ture range  as  ^ (T)  becomes  smaller. 

In  order  to  calculate  certain  properties  of  our  bridges 
such  as  the  critical  current  or  their  approximate  size,  an 
appropriate  geometry  should  be  considered.  Aslamazov-Larkin 
consider  a hyperbolic  neck  of  radius  a and  apex  angle  6 as 
shown  in  Fig.  3.1  (a).  As  long  as  the  mean  free  path  £ « a, 
the  normal  resistance  of  the  link  is  given  by 

ctn(e/2)  (3.4) 

where  p is  the  normal  state  resistivity.  It  is  convenient 
for  calculations  to  rejjlace  this  geometry  by  the  more  arti- 
ficial one  shown  in  Fig.  3.1  (b) , defined  by  two  intersecting 
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j Figure  3. 

r 


Schematic  of  a)  hyperbolic  neck  geometry; 
b)  collinear  cone  geometry 
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collinear  cones  with  tlie  same  angle  20,  which  intersect  at  a 
radius  r^^.  Tlie  radius  r^^  is  chosen  so  that  the  normal  state 
resistance  for  this  geometry,  given  by 


^n  Q r- 


(3.5) 


where  = 4TTsin  (6/2)  is  the  solid  angle  for  cooling,  equals 
that  of  Eq.  (3.4).  Then  r^^  = 2a/rr  sin  0. 

Care  should  be  taken  in  attempting  to  extract  informa- 
tion about  the  physical  size  of  our  microbridges  from  Eq. 
(3.4)  or  (3.5),  since  they  were  derived  under  the  assumption 
that  the  mean  free  path  I « a,  the  neck  radius.  This  con- 
dition will  only  be  satisfied  for  very  low  resistance  micro- 
bridges. The  lowest  resistance  microbridge  studied  in  this 
work  had  = 0.06  Q;  from  which  a value  of  a = 0.14  qm  was 
obtained  using  Eq.  (3.5).  Since  our  films  have  typical  mean 
free  paths  I of  the  order  of  0.1  jam,  it  is  necessary  to 
relax  the  condition  £ « a by  explicitly  including  the 
effects  of  boundary  scattering.  In  that  case,  the  effective 
mean  free  path  can  be  approximated  as 


= i + 


2 r sin  6 


(3.6) 


making  the  resistivity  a position  dependent  quantity 


p ^-^(r)  = p (1  + 5 — ) 

^eff  ^o'  2 r sin  0 


(3.7) 
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The  normal  state  resistance  of  the  weak  link  will  then  be 
given  by 

2o„ 

(3.8) 


Using  Eq.  (3,8),  the  calculation  of  the  neck  radius  ^ for 
the  example  of  = 0.06  ohm  yields  a = 0.18  um.  Thus,  even 
for  our  lowest  resistance  microbridges,  the  corrections  due 
to  boundary  scattering  are  significant.  These  corrections 
will  be  quite  large  for  our  high  resistance  microbridges  or 
high  resistance  point  contacts.  The  collinear  cone  geometry 
should  be  a very  good  approximation  to  the  geometry  of  our 
variable-thickness  microbridge  as  evidenced  by  the  SEM  pic- 
tures shown  in  Chapter  Two.  In  some  cases  discussed  in 
Chapter  Five,  it  might  be  necessary  to  take  into  account  a 
small  finite  length  of  the  bridge  neck  between  the  cones,  but 
this  can  be  neglected  in  calculating  resistances  or  critical 
currents . 

In  Table  3.1,  some  of  the  basic  properties  of  our  tin 
microbridges  are  presented.  (A  more  complete  table  of  all 
the  properties  of  our  tin  and  lead  microbridges  is  included 
in  Appendix  A.)  The  bridges  are  arranged  in  Table  3.1  in 
order  of  increasing  resistance.  From  SEM  pictures  of  the 
microbridges  it  was  determined  that  they  do  not  vary  signi- 
ficantly in  length  or  width  at  the  narrowest  point  of  the 
neck.  Unfortunately,  the  thickness  of  the  bridge  depends  on 
the  groove  depth  which  can  vary  with  pressure  or  the  shape 
of  the  particular  diamond  point  used. 
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'Kible  3.1 

Characteristic  Properties  of  Some  Tin 
Variable-Thickness  Microbridges 


# 

Normal  State 
Resistance 
(Ohms) 

dl 

(mV/K) 

dT 

Exponent  of 

I oc(T-T  )°^ 

CO  c 

(mV) 

19 

0.06 

0.76 

0.94 

0.88 

21 

0.07 

0.77 

1.25 

1.10 

13 

0.09 

0.68 

0.87 

1.16 

14 

0.20 

1.16 

1.00 

1.38 

4 

0.24 

1.37 

1.07 

2.71 

11 

0.33 

0o56 

0.97 

1.32 

9 

0.36 

0.56 

1.07 

0.90 

16 

0.44 

0.80 

0.96 

1.19 

8 

0.44 

0.75 

1.09 

3.30 

7 

0.43 

1.02 

1.03 

1.63 

10 

0.85 

1.27 

0o87 

1.62 

12 

1.14 

1.05 

1.05 

1.49 

17 

1.25 

0.69 

- 

1.25 

20 

1.76 

0.42 

0,91 

1.44 
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The  range  of  resistances  in  Table  3.1  is  representative 
of  the  bridges  tested,  which  of  course  were  influenced  by 
the  author's  procedure  for  optically  selecting  oridges  as 
well  as  by  the  survival  rates  of  the  bridges  selected.  In 
• the  earlier  stages  of  this  work,  low  resistance  bridges 

were  used  because  of  their  larger  survival  rate.  As  the 
burnout  problem  was  brought  under  control,  the  smallest  pos- 
sible bridge  on  each  substrate  was  tested.  All  of  the  films 
from  which  the  bridges  were  made  had  mean  free  paths  of 
order  0.1-0.14  um,  so  that  differences  in  resistance  from 
1 bridge  to  bridge  are  mostly  due  to  size.  Table  3.1  illus- 

!'  trates  the  difficulty  in  making  microbridges  of  clean  tin 

' with  resistances  larger  than  1 ohm,  as  we  usually  discarded 

those  bridges  with  low  resistances  at  room  temperature  and 
I these  were  already  a selection  of  the  smallest  possible 

microbridges. 

[ In  some  cases  the  bridges  were  longer  than  0.5  |am. 

Tliis  occurred  when  the  bank  on  the  side  of  the  bridge  was 
lifted  off  the  substrate  by  the  diamond  knife  leaving  some 

I 

tin  in  the  groove  under  the  lifted  film  and  thus  creating  a 

4 longer  microbridge.  We  believe  this  causes  the  high  resis- 

tance observed  in  sane  bridges  such  as  #10,  12  and  20.  Thus, 
of  the  bridges  with  no  anomalous  properties,  bridge  #17  has 
the  highest  resistance;  using  Eq.  (3.8)  this  corresponds  to 
a neck  radius  of  0,02  ^m. 

Using  the  results  derived  for  the  collinear  cone  geometry 
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it  is  then  possible  to  calculate  the  magnitude  and  tempera- 
ture dependence  of  the  critical  current  Eq.  (3.3).  If  one 
assumes  the  electrical  potential  satisfies  Laplace's  equation 

tlien  V(r)  = Vj^+  Vf(r)  where  V is  the  total  potential  drop, 

2 2 

V^-  n-ien  since  I = rra  oWf  = ~ 1/ ovf)  and 

(3.3)  reduces  to 


e*^,  |Y 

I R s=  - ^ 

CO  N m*  a 

In  the  dirty  limit  (£  « P^)  this  reduces  to 


(3.9) 


c 

where  A (T)  is  the  energy  gap  and  A(0)  = 1.76  T^,  where 

temperature  is  measured  in  energy  units.  Eq.  (3.10)  agrees 

near  with  the  expression  for  the  critical  current  obtained 

23 

for  tunnel  junctions  by  Ambegaokar  and  Baratoff: 


I tanh[A  (T)/2Tl  (3.11) 

In  column  2 of  Table  3.1  the  experimental  results  for 

R^dl^^/dT  near  are  shown.  Note  that  these  values  are 

usually  higher  than  would  be  expected  from  Eq.  (3.10)  and 

may  exceed  the  expected  value  by  as  much  as  a factor  of 

four.  This  is  a common  discrepancy  found  in  two-dimensional 

24 

tin  and  indium  microbridges  as  well  as  in  variable- 

25 

thickness  microbridges.  The  origin  of  this  discrepancy  is 


not  presently  understood.  It  is  apparently  not  due  simply 
to  size  effect  corrections  to  the  resistivity,  since  these 
would  be  most  important  in  high  resistance  weak  links.  It 
is  also  not  due  to  uncertainty  in  determining  the  appropriate 
experimental  resistance  values,  because  the  theoretical 
is  the  full  normal  resistance  of  the  entire  weak  link, 
surely  an  upper  bound. 

The  more  detailed  temperature  dependence  of  the  criti- 
cal current  for  two  of  our  tin  microbridges  is  shown  in 

Fig.  3.2.  The  values  of  (dl  /dT)R„  shown  in  Table  3.1  are 

CO  N 

determined  from  the  linear  region  between  T/T^  = 0.8-0.98. 

The  value  of  the  exponent  obtained  from  a least  squares  fit 
over  this  temperature  range  is  shown  in  column  3 of  Table  3.1. 
These  values  are  in  very  good  agreement  with  what  is  expec- 
ted for  Josephson  behaviour  near  T^.  For  temperatures  below 
T/T^  ~ 0.8  the  critical  current  begins  to  grow  less  rapidly 
as  predicted  by  Eq.  (3.11) 

Near  T^ (T/T^  ^ 0.98)  the  situation  is  quite  different. 

Two  different  types  of  behaviour  can  be  observed  depending 

on  the  sample,  as  shown  in  the  inset  in  Fig.  3,2.  The  more 

common  behaviour  is  that  of  the  lower  curve  which  has  a small 

nonlinear  tail  very  close  to  T , where  T is  identified  as 

c c 

that  temperature  at  which  a supercurrent  is  first  observed. 

We  believe  this  behaviour  is  caused  by  the  bridge  region 
having  an  intrinsic  transition  temperature  T^^  below  that  of 
the  bank  material  This  is  supported  by  careful 
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measurements  of  the  resistive  transition  of  our  samples, 
in  which  a resistive  change  of  order  0.08  ohms  occurs  first 
as  the  bridge  is  cooled  down.  This  resistive  change  cor- 
responds to  that  expected  from  the  contribution  of  the  films 
between  the  superconducting  leads  and  the  microbridge.  The 
differences  between  the  transition  temperature  of  the  banks 
T^^  and  the  observed  T^  vary  from  sample  to  sample  and  can 

be  as  large  as  20  mK.  Using  the  results  obtained  by  Likharev 

23 

and  Yakobson  for  the  critical  currents  of  SNS  junctions  we 
calculate,  for  L/5 (0)  ~ 5,  that  the  observed  T^  will  be 
approximately  proximity  effect,  so 

that  the  intrinsic  transition  temperatures  of  the  bank  and 
bridge  materials  differ  by  as  much  as  40  mK.  Such  differ- 
ences are  not  at  all  surprising.  The  transition  temperature 
of  our  tin  films  (T^  ~ 3.8  K)  is  always  higher  than  that  of 
bulk  tin  (T^  = 3.722  K)  due  to  stresses.  In  view  of  the 
fabrication  technique  used  in  making  the  bridges,  it  is 
expected  that  the  whole  region  around  the  bridge  will  be 
under  very  nonuniform  stress  conditions.  This  leads  to  the 

observed  behaviour  which  is  typically  seen  in  tin  and  indium 

24  25 

microbridges  ' and  in  our  lead  microbridges. 

In  some  cases,  the  observed  behaviour  corresponds  to 

the  upper  curve  in  the  inset  of  Fig.  3.2.  We  interpret  this 

as  resulting  from  the  bridge  material  having  a higher  intrinsic 

T^  than  the  banks,  Ihis  is  less  common  in  tin  microbridges 

25 

but  is  usually  observed  in  aluminum  microbridges. 
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In  Fig.  3.3  the  temperature  dependence  of  for  a 

lead  variable-thickness  microbridge  is  shown.  Lead  micro- 
bridges have  received  little  attention  because  it  was  assumed 
that  the  short  coherence  length  [f  (0)  ~ 0.08  (jm]  of  lead 
would  make  them  behave  more  like  bulk  filaments.  However, 
the  temperature  dependence  of  the  critical  current  is 
essentially  the  same  as  that  shown  earlier  for  tin  micro- 
bridges. The  linear  " Josephson-like " behaviour  is  observed 
over  the  same  temperature  range  ~ 0.8  - 0.98.  This 

linear  temperature  dependence  of  is  to  be  contracted 

with  earlier  experimental  results  on  tin  and  indium  micro- 
bridges  ' in  two-dimensional  geometries  which  had  criti- 
cal currents  which  were  somewhere  between  Josephson  (T-T 

c 

1 5 

and  bulk-like  (T-T^)  * behaviour.  This  advantage  of  the 

variable-thickness  geometry  has  also  been  confirmed  recently 
27 

by  Wong  et  al,  who  observe  linear  Josephson-like  behaviour 
in  fairly  long  (~  1 pm)  variable-thickness  microbridges  of 


niobium  (which  also  has  a very  short  coherence  length). 
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Recent  calculations  by  Daalms  et  al  show  that  a small  two- 

1 25 

dimensional  hyperbolic  neck  should  have  I a (T-T  ) * and 

CO  c 

for  a three-dimensional  neck  they  obtain  that  will  vary 


linearly  as  outlined  above  or  as  in  calculations  which 

29-31 

assume  a short,  one-dimensional  bridge. 


While  the  overall  behaviour  of  the  critical  current 


for  both  our  tin  and  lead  microbridges  is  qualitatively 
Josephson-like,  some  serious  quantitative  discrepancies  in 
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40 


the  I R product  near  T remain.  In  Eci.  (3.9),  the  con- 

ductivity  should  bo  a position-dependent  quantity  when  the 

mean  free  path  is  of  the  order  of  the  nock  radius  How- 

ever,  |y^|  is  also  proportional  to  I when  I is  less  than 

tlio  coherence  length,  and  these  effects  should  cancel  to  a 

first  approximation.  Tliis  qualitative  argument  has  recently 

32 

been  confirmed  by  calculations  made  by  Kulik  et  al.  which 

show  that  the  behaviour  does  not  differ  from  Eq.  (3.10)  near 

T , for  either  f > ? or  f < - as  long  as  a « " . 

C 'O  'O  'O 

3.1.2  Low  Temperature  Regime:  Products 

Tfie  temperature  dependence  of  the  critical  current  of 
short  microbridges  far  from  T^  has  received  little  atten- 
tion both  theoretically  and  experimentally.  Experimentally 
it  is  always  problematic  to  measure  the  critical  current  in 
the  presence  of  hysteresis  since  premature  switching  to  the 
resistive  state  might  obscure  the  true  critical  current. 

This  temperature  regime  is,  however,  quite  important  in 
terms  of  practical  applications  at  high  frequencies  since 
the  product  defines  a characteristic  frequency 

<x  = I above  which  the  resistive  shunting  dominates 

the  Josephson  element,  resulting  in  poorer  performance.  It 
is  usual  to  compare  experimental  results  at  low  temperatures 
to  Eq.  (3.11),  because  of  its  agreement  with  the  A-L  result 
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low  tompora tures : 

Wliile  this  form  agrees  well  with  experimental  products 

33 

in  tunnel  junctions  and  high  resistance  point  contacts, 

the  same  is  not  true  of  microbridges  and  low  resistance 
33  34 

point  contacts.  ' In  Fig.  3.4  we  show  the  products 

at  2.1  K for  our  tin  variable-thickness  microbridges.  All 
of  our  experimental  results  arc  above  those  expected  from 
Eq.  (3.11),  which  is  shown  as  a solid  line  in  Fig.  3.4. 

These  discrepancies  appear  to  be  largest  for  our  low  resis- 
tance microbridges.  This  same  trend  is  also  observed  by 
34 

Yanson  in  tin  microbridges  fabricated  by  electrical  break- 
down of  metal-dielectric-metal  structures.  Calculations  by 
35  3 6 

Kulik  et  al.  and  Mitsai  from  microscopic  theory  in  the 

limits  I « 5 and  a > £ give  a value  of  I at  T = 0 

which  is  1.32  times  larger  than  that  given  by  Eq . (3.12). 

32 

More  recent  calculations  by  Kulik  and  Omelyanchuk  in  the 
limit  a « and  a « £ (for  either  £ > ?^  or  £ < e^) , 
which  correspond  to  limits  closer  to  the  experimental  ones, 
yield  a low  temperature  value  of  which  is  twice  as 

large  of  Eq.  (3.12).  As  seen  from  Fig.  3.4,  our  results 
agree  better  with  this  prediction  although  some  of  the 
experimental  results  are  still  somewhat  higher  than  the  highest 
expected  value.  While  those  predictions  arc  in  closer 


“2.1K  , _ tA(T)  , A(T) 


42 


Fip.  3.4  I R.,  oroducts  versus  R,,  for  tin  variable- thickness  microbridges. 
CO  f'  n 
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agreement  with  our  experimental  results,  the  question  re- 
mains of  why  this  last  prediction  should  be  valid  in  micro- 
bridges and  not  in  high  resistance  contacts  which  have  the 
same  ranges  of  neck  radius  a compared  to  and  I,  but 
nevertheless  show  I R,  products  which  are  below  or  in 
agreement  with  Eq,  (3.11).  It  should  be  pointed  out  that 
the  highest  values  at  low  temperatures  do  not  occur  in 

those  bridges  with  the  highest  I values  near  T as  shovm 

in  Table  3.1. 

It  is  clear  that  while  recent  theoretical  predictions 
agree  better  with  experimental  results  at  low  temperatures 
than  the  Ambegaokar  and  Baratoff  result,  further  attention 
to  this  problem  is  required  for  conditions  X ~ ~ a and 

in  geometries  closer  to  those  of  experiments. 

3.2  The  Wyatt-Davem  Effect 

One  of  the  most  interesting  effects  observed  in  most  typos 
of  Josephson  junctions  is  the  enhancement  of  the  critical  cur- 
rent I induced  by  microwave  radiation.  Tliis  effect  was 

CO 

first  noted  in  rather  large  (W  ~ 3 pm,  L ~ 3 pm)  two- 

37  38 

dimensional  microbridges  by  Wyatt  et  al.  and  Dayem  et  al. 

39 

though  it  had  previously  been  observed  by  Shapiro  in  tunnel 

junctions.  While  this  effect  has  now  been  observed  in  most  junc- 

40-42  43 

tion  geometries,  and  in  the  presence  of  phonon  radiation, 

no  complete  description  of  the  effect  for  all  geometries  and 

systems  is  available.  It  is  now  generally  accepted  that 

this  effect  originates  at  least  in  part  from  an  enhancement 
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of  tJie  gap  resulting  from  the  microwave  excitation  of 

quasi-particles  to  higher  energies.  This  explanation  was 

44 

first  proposed  by  Eliashberg.  However,  his  calculations 
were  made  for  a thin  superconducting  film  without  any  current 
flov.’ing  through  it  and  they  do  not  necessarily  apply  to  our 
variable-thickness  microbridges. 

Experimental  results  on  the  Wyatt-Dayem  effect  in  small 
weak  links  have  shown  a wide  variation  of  temperature,  power 
and  frequency  dependence.  The  results  in  our  microbridges 
are  no  exception,  as  they  differ  from  all  previously  published 
results  for  either  the  two-dimensional  or  variable-thickness 
geometry.  All  of  our  microbridges  show  an  enhancement  of 
I at  both  10  and  32  GHz.  At  low  microwave  powers  I 

CO  CO 

increases  slowly  with  power,  reaches  a maximum  and  then 

gradually  decreases  in  a manner  consistent  with  a decrease  of 
given  by  the  zero-order  Bessel  function  of  the  normalized 
microwave  voltage.  The  magnitude  of  the  enhancement  varies 
from  bridge  to  bridge  from  4-43%  and  appears  to  depend  on  the 
weak  link  resistance.  (In  Appendix  A,  the  maximum  enhance- 
ment observed  in  all  samples  is  included  in  the  table  of 
properties  of  all  bridges;  however,  this  is  the  largest 
enhancement  observed,  ratiier  than  an  absolute  maximum,  as  in 
many  bridges  the  Wyatt-Dayem  effect  was  not  studied  in 
detail. ) 

CXir  observation  of  the  Wyatt-Dayem  effect  in  our  bridges 
agrees  with  the  anticorrelation  found  by  some  workers  between 
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good  Josephson  behaviour  and  microwave  enhancement  of  the 

45 

critical  current.  Jillie  et  al.  found,  that  among  two- 
dimensional  microbridges,  bridges  with  the  smallest  enhance- 
ment exhibited  the  best  Josephson  properties.  Klapwijk  and 
46 

Mooij  did  not  observe  the  enhancement  in  their  high  per- 
formance tin  and  aluminum  variable-thickness  microbridges. 

In  our  case,  the  enhancement  is  largest  in  low  impedance 
microbridges  which,  as  shown  in  Chapter  Five,  have  their 
Josephson  properties  most  severely  affected  by  heating. 

Thus  good  diffusion  of  the  quasiparticles  appears  to  inhibit 
the  Wyatt-Dayem  effect. 

The  more  detailed  temperature  dependence  of  the  maxi- 
mum enhancement  observed  in  our  microbridges  is  shown  in 
Fig.  3„5  for  both  10  and  32  GHz.  These  results  are  similar 
to  those  observed  in  most  of  our  microbridges.  At  both 
frequencies,  the  enhancement  is  not  observed  very  close  to 
T^,  As  the  temperature  is  lowered,  it  appears  first  for  10 
GHz  and  later  for  32  GHz.  As  the  temperature  is  reduced 
further  the  enhancement  at  32  GHz  becomes  larger  than  that 
at  10  GHz.  The  enhancement  at  both  frequencies  reaches  a maxi- 
mum at  T/T^  Pa  0.97  and  decreases  and  disappears  at  lower 
temperatures,  with  the  enhancement  at  10  GHz  going  to  zero 
first.  However,  some  other  experiments  in  weak  links  quite 

similar  to  ours  have  yielded  significantly  different  results. 

46 

Latyshev  and  Nad  observed  the  appearance  of  a critical  cur- 
rent above  T at  both  10  and  23  GHz  in  cross-scratched  tin 
c 
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Fi(?.  3.5  Enhancement  of  the  critical  currents  at  10  and  32  GHz  for  bridcre  #21. 
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microbridges,  although  not  at  higher  frequencies.  Levin- 
47 

sen  observed  the  enhancement  in  cross-scratched  tin  micro- 
bridges only  below  at  both  10  and  32  GHz  but  with  a 
temperature  dependence  quite  different  than  ours.  At  10  GHz 
the  enhancement  began  very  close  to  T^,  reached  a maximum  at 
T/T^  = 0.997  and  then  rapidly  decreased  to  zero.  At  32  GHz 
the  enhancement  was  observed  first  below  T/T  = 0.99  , 
reached  a maximum  at  0.98  and  remained  many  orders  of  magni- 
tude larger  than  enhancement  at  10  GHz  at  least  down  to 
T/T^  = 0.95.  At  present  we  do  not  know  what  accounts  for 
these  discrepancies. 

While  none  of  the  experimental  results  can  be  said  to 
be  inconsistent  with  Eliasnberg's  idea  as  to  the  origin  of 
the  effect,  the  experimental  details  differ  from  theoretical 
predictions.  Ihe  model  predicts  the  enhancement  will  occur 
for  frequencies  of  order  uu  > 1/t^  whore  is  the  electron- 
phonon  inelastic-scattering  time.  Tliis  is  certainly  true  in 
our  tin  microbridges  ^ ~ 10“^®  sec.)  and  is  consistent 
with  the  absence  of  the  enhancement  in  our  lead  ('2  ~ 10~'^^ 

sec.)  variable-thickness  microbridges  and  with  the  onset 

—8 

of  the  effect  at  200  MHz  in  long  aluminum  (t ^ ~ 10  sec.) 

48 

microstrips.  However,  we  do  not  observe  the  increase  of 
the  critical  temperature  predicted  by  theory  at  the  fre- 
quencies and  temperatures  of  our  experiment.  Such  an  in- 

48 

crease  is  always  seen  in  long  A1  and  Sn  strips 
short  microbridge  experiments. 


and  in  some 
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Wliile  it  might  be  possible  that  a proper  theoretical 
treatment  of  the  problem  including  dc  and  rf  currents  in  a 
constricted  geometry  would  agree  with  most  of  the  experi- 
mental results,  it  appears  unlikely  that  it  will  resolve  the 
inconsistencies  between  experiments.  This  is  best  illustra- 
ted by  Ihble  3.2  where  we  show  the  combinations  of  material 
and  geometry  for  which  the  effect  has  been  observed  at 
different  frequencies.  Note  in  the  table,  that  the  two 
systems  one  expects  to  be  similar,  the  variable-thickness 
bridge  and  the  point  contact,  show  exactly  opposite  results. 

It  would  appear  from  the  above  results  that  any  further 
understanding  of  the  problem  will  require  theoretical  analysis 
of  more  realistic  configurations  in  closer  agreement  with 
the  experimental  conditions.  Experimental  differences  could 
be  identified  by  varying  the  details  of  the  experiment  such 
as  the  geometry  of  the  bridges  or  tlie  substrates  used.  It 
would  also  be  interesting  to  study  materials  with  very 
T^'s  at  much  higher  frequencies 


fast 
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Compar 

Aluminum ; 

Tin : 

Aluminum : 

Tin: 


Aluminum : 
Tin ; 

Lead ; 

Aluminum : 
Tin; 
Indium : 


Table  3,2 

ison  of  the  Wyatt-Dayem  Effect  Observed  in 
Different  Materials  and  Geometries 


Long  Microbridqcs 

Enhancement  up  to  300%,  Enhancements  from  200  MHz 

to  10  GHz,  Increase  of  T observed, 

c 

Small  (a  few  %)  enhancements  of  I at  10  GHz.^ 

CO 

Two-Dimensional  Short  Microbridges 


Small  enhancements  at  35  GHz  down  to 

49 

No  increase  of  T reported, 
c 


T/T 

c 


0,86, 


Various  results  reported.  Large  enhancements 
observed  in  some  cases.  Increase  of  T seen  a 


some  frequencies.  Enhancements  below  T seen 

. T 37,38,46,47  ^ 

from  1-76  GHz, 


t 


Variable-Thickness  Microbridges 

50 

No  enhancement  at  35  GHz, 

Up  to  50%  enhancements  at  10  and  32  GHz,,  No  in- 
crease of  T observed  (this  work)  but  no  enhance- 

^ 50 

ments  observed  by  others. 

No  enhancement  observed  at  10  and  32  GHz  (this  work). 


Point  Contacts 

5-10%  enhancements  from  1 to  10  GHz,^^ 

41 

No  enhancements  observed, 

41 

No  enhancements  observed. 


Niobium : Small  enhancements  reported 
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Ibble  3.2 


Aluminum ; 


50 


(continued) 


Tunnel  Junctions 

39 

Small  enhancements  of  I reported. 

CO 

Direct  enhancement  of  the  energy  gap  by  as  much 
as  100%. 

Increase  of  T reported, 
c 
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CHAPTER  FOUR 


NONEQUILIBRIUM  ENHANCED  SUPERCURRENTS 
AT  NONZERO  VOLTAGES 


4.1  Intrcxiuction 

Once  the  critical  current  of  a weak  link  is  exceeded, 

a voltage  develops  across  its  length  and  the  superconductor 

is  in  a time-dependent,  dissipative  state  in  which  non- 

etiuilibrium  effects  are  very  important.  Tliis  dynamic  state 

arises  from  the  superconductor  responding  to  the  voltage  in 

such  a way  that  phase  coherence  is  maintained  across  the 

18 

weak  link.  Then,  as  shown  by  Aslamazov  and  Larkin,  the 
order  parameter  in  the  center  of  the  weak  link  will  oscillate 
between  the  value  of  the  order  parameter  at  the  banks,  to 
zero  and  then  back  to  the  bank  value  in  a Josephson  period 
"j.  This  oscillating  order  parameter  or  energy  gap  will 
determine  the  energy  of  the  quasiparticles  and  if  t is 

U 

fast  compared  to  the  relaxation  time  of  the  quasiparticles 
^ 2'  the  quasiparticle  occupation  numbers  will  lag  behind 
the  energy  changes,  leading  to  changes  in  the  superconducting 
properties  of  the  weak  link.  At  higher  voltages  (and  thus 
frequencies)  the  order  parameter  will  be  limited  by  t and 

GLi 

will  not  be  able  to  relax  instantaneously  and  further  non- 
equilibrium corrections  to  the  superconducting  properties 
of  the  weak  link  will  occur.  At  still  higher  voltages, 

4 

Joule  heating  will  be  of  importance  and  will  eventually 
dominate  all  nonequilibrium  effects. 

51 
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In  this  chapter  we  discuss  the  low  voltage  regime  of 
the  I-V  characteristics  of  our  tin  variable-thickness  micro- 
bridges, where  heating  is  not  of  importance,  and  where  a 
temperature  dependent  foot-like  feature  is  observed.  We 
find  this  feature  to  be  consistent  with  the  idea  of  an 
enhancement  of  the  order  parameter,  and  thus  the  supercurrent 
flowing  through  the  weak  link,  due  to  the  long  relaxation 
time  of  the  quasiparticles.  In  section  4.2  we  describe  in 
detail  the  phenomena  observed  experimentally.  In  section 
4.3  we  present  first  a qualitative  discussion  of  the  origin 

of  this  effect,  before  discussing  in  detail  the  more  involved 

B 7 

models  of  Aslamasov  and  Larkin  and  of  Golub.  The  pre- 
dictions of  these  models  are  compared  to  the  experimental 
results  in  section  4.4  and  good  agreement  is  found  with  many 
aspects  of  the  experimental  results  on  our  tin  microbridges. 

4.2  Ttie  Experimental  Phenomena 

The  experimental  phenomena  to  be  explained  have  been 

observed  in  two-dimensional  and  variable-thickness  micro- 

28  5 2—54 

bridges  of  tin  and  indium  made  in  various  laboratories.  ' 

In  Fig.  4.1,  we  show  the  overall  I-V  characteristics  of  a 
tin  microbridge  near  T^.  At  high  voltages  there  is  a sub- 

4 

stantial  amount  of  curvature  due  to  heating.  At  lower 
voltages  one  observes  the  characteristic  features  of  gap 
and  subharmonic  gap  structures  which  will  be  discussed  in 
Chapter  Five.  Finally,  at  the  lowest  voltages  (0-50  qV) 


T(K) 


nredicted  by  Aslamazov  and  Larkin 
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there  is  a characteristic  "foot",  shown  in  detail  in  Fig. 

4.2.  Near  T , the  change  to  the  resistive  state  above  I 

C CO 

is  simply  marked  by  a steep  rise  in  voltage  and  a curve 

without  inflection  points.  The  curve  shows  an  apparent 

excess  current  defined  as  the  zero-voltage  intercept 

extrapolated  from  the  higher  voltage  parts  of  the  curve. 

Near  T , I is  a fraction  of  the  order  of  0.6-0. 8 times  the 
c s 

critical  current  At  somewhat  lower  temperatures  this 

feature  develops  continuously  into  an  initial  region  (up 

to  5 aV)  of  constant  differential  resistance,  which  is  much 

smaller  than  the  normal  resistance  of  the  weak  link, 

follov/ed  by  a rapid  growth  of  the  voltage  with  very  little 

current  change  at  a characteristic  current  I , . Above  this 

cl 

region,  there  is  an  apparent  excess  current  which  now 
appears  to  be  related  to  rather  than  ss  its  magni- 

tude is  larger  than  I . Because  of  the  well  cooled  nature 

CO 

of  our  variable-thickness  geometry,  heating-related  hyster- 
. 4 

esis  appears  only  at  low  temperatures,  which  allows  us  to 
follow  the  development  of  this  "foot"  over  a wide  temperature 
range.  We  find  that  the  ratio  I , /I  . increases  as  the  tem- 
perature  is  decreased,  reaching  a maximum  ranging  from  1.2-2 
in  our  microbridges  at  temperatures  of  the  order  of  3.2  K. 

At  still  lower  temperatures  this  ratio  decreases,  an  effect  that 
can  in  the  most  part  be  attributed  to  heating.  The  magnitude 
and  temperature  dependence  of  the  differential  resistance 
at  low  voltages,  of  the  ratio  ^q^^^co'  apparent 
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excess  supercurrent  are  the  major  experimental  facts  to  be 
accounted  for. 

4 . 3 Theoretical  Interpretation 

In  section  3.1.1  a simple  model  of  the  equilibrium 

behavior  of  weak  links  such  as  metallic  point  contacts  and 

variable-thickness  microbridges  was  introduced.  In  the 

form  first  proposed  by  Aslamasov  and  Larkin,  the  complex 

order  parameter  in  the  banks  is  assumed  to  have  the  values 
icp 

Y and  Ye  ° where  m is  the  phase  difference  between  the 
o o o 

two  sides  and  to  vary  linearly  across  the  one-dimensional 
bridge  of  length  L connecting  them: 

icp 

Y = Yo  [1  - f + G ° ^]  . (4.1) 

Thus  the  magnitude  of  the  order  parameter  in  the  very  center 
of  the  bridge  oscillates  at  the  Josephson  frequency  between 
a maximum  Y^,  equal  to  the  value  in  the  banks,  to  a minimum 
of  zero  and  then  back  to  Y^.  As  the  order  parameter  goes 
to  zero,  a phase  slip  of  2-^  occurs  at  the  center  of  the  bridge, 
but  no  physical  discontinuity  occurs  as  in  the  relaxation 

oscillation  model  or  the  phase  slip  models  applicable  to 
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long  micr obr idges . 

Ihe  above  model  has  served  as  the  starting  point  for 
several  previous  considerations  of  the  effect  of  the  relax- 
ation of  the  order  parameter  based  on  time-dependent  Ginzburg- 
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Landau  (TDGL)  theory  and  ad  hoc  extensions  of  it.  Likharev 

and  Yakobson^  and  Kramer  and  Bara toff  have  numerically 

solved  the  TDGL  equations  for  current-biased  microbridges  of 

different  lengths.  Relaxation  effects  become  important  for 

times  comparable  to  the  Ginzburg-Landau  time  = nh/8 (T^-T) , 

which  becomes  significant  primarily  at  high  voltages. 

Likharev  and  Yakobson  emphasized  that  end  effects  near  the 

banks  cause  an  "insufficient  voltage",  or  apparent  "excess 

current"  of  magnitude  0.75  at  high  voltages,  and  found 

the  appearance  of  an  additional  phase-dependent  conductance 

which  had  little  effect  on  the  I-V  curves.  Jensen  and 
5 2 

Lindelof  subsequently  noted  that  the  phase-dependent  con- 
ductance would  have  a much  greater  effect  on  the  I-V  charac- 
teristics if  the  relaxation  time  were  taken  to  be  two  orders 
of  magnitude  longer;  it  then  leads  to  features  in  the  I-V 
characteristics  which  are  qualitatively  similar  to  the  foot 
of  interest  here.  However,  the  use  of  such  long  times  was 
given  no  physical  justification,  and  Kramer  and  Baratoff  have 
noted  that  it  loads  to  inconsistencies  within  the  analytic 
approximations  that  Jensen  and  Lindelof  use.  Furthermore, 
the  use  of  TDGL  is  only  strictly  justified  in  the  gapless 
regime  which  is  not  appropriate  for  these  bridges  at  least 
near  equilibrium  at  low  current  densities.  Leaver  et  al.^^'^^ 
have  noted  similar  effects  on  the  I-V  characteristics  result- 
ing from  time  lags  in  a simple  phenomenological  nonequilibrium 


model 
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Recent  theories  which  go  beyond  the  simple  TDGL  des- 
cription take  explicit  account  of  the  relationship  between 
the  nonequilibrium  occupation  of  quasiparticle  states  and 
changes  of  the  superconducting  order  parameter.  We  will 
first  suggest  qualitatively  how  such  a relationship  can  lead 
to  an  enhancement  of  the  supercurrent  and  then  we  will  review 
the  more  specific  calculations  by  Golub  and  by  Aslamasov  and 
Larkin  on  which  this  picture  is  based. 


4.3.1  Qualitative  Picture 

To  start,  let  us  assume  that  the  simple  variation  of  the 

order  parameter  Eq.  (4.1)  also  describes  the  space  and  time 

variations  of  the  energy  gap,  whose  magnitude  determines 

the  energy  of  each  excitation  through  the  relation  E = 

2 2^ 

+ A ) , where  is  the  energy  relative  to  the  Fermi 
energy  E^  of  an  excitation  in  momentum  state  k in  the  absence 
of  the  pairing  interaction.  The  occupation  of  states  with  a 
given  energy  will  relax  toward  the  equilibrium  value  given 
by  the  Fermi  function  f^lE)  = [l  + expCsE)]"^  by  two  pro- 
cesses, inelastic  scattering  (primarily  by  phonons)  to  and 
from  states  of  other  energy,  and  spatial  diffusion  at  constant 
energy  to  and  from  adjacent  regions  with  different  occupa- 
tions. Diffusion  over  the  distance  L will  occur  over  a 
characteristic  time  = L /ttD,  where  D is  the  diffusion 

coefficient;  for  dimensions  of  order  one  micron,  Tjj  is  com- 

2 —12 
parable  to  = F /D,  i.e„,  of  order  10 


seconds  for  tin. 
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Thus  at  energies  which  exceed  the  order  parameter  in  the 

banks  the  diffusive  exchange  of  particles  with  the  banks 

will  fix  the  occupation  of  these  states  at  f^(E).  Here  we 

are  neglecting  the  effects  of  Joule  heating  which  would 

tend  to  increase  the  occupation  of  these  states.  Such 

effects  are  minimized  in  the  variable-thickness  geometry 

(which  allows  diffusion  in  three  dimensions  away  from  the 

bridge),  and  are  insignificant,  close  to  T^,  at  the  low 

voltage  levels  of  interest  here. 

When  the  magnitude  of  the  energy  gap  in  the  bridge  is 

depressed  below  A^,  then  can  be  below  A^,  and  except  for 

possible  level  broadening,  these  states  can  not  exchange 

particles  with  the  banks  as  shown  in  Fig.  4.3(a).  In  this 

case,  the  only  mechanism  for  relaxing  the  occupation  of 

those  states  with  energies  below  A^  involves  inelastic 

phonon  scattering,  which  is  characterized  by  a relatively 

-9  -10 

long  time  which  for  tin  is  of  order  10  - 10  seconds. 

Thus  we  should  expect  different  regimes  of  behaviour  depend- 
ing on  whether  the  Josephson  period  is  long  or  short  compared 

to  Tg. 

At  very  low  voltages  where  the  Josephson  period  T , is 
much  longer  than  Tg,  the  relaxation  will  occur  quickly  on 
the  scale  of  t , and  the  corrections  to  the  equilibrium  model 

U 

will  be  small.  This  is  shown  qualitatively  in  Fig.  4.3(b)- 
(d).  As  the  Josephson  period  decreases  with  increasing 
voltage,  the  rate  of  change  of  the  order  parameter  will 


Spatial  variation  of  the  gap  in  the  middle 
of  the  Josephson  cycle  showing  the  quasi- 
particle  states  that  can  not  diffuse  into 
the  banks. 

Temporal  variation  of  the  gap  in  equilibrium 
(solid  line)  and  in  the  presence  of  a 
lag  (dashed  line) . 

Temporal  variation  of  excitation  energy 

Temporal  variation  of  the  quasi-particles 
occupation  numbers  for  different  values 
of  the  Josephson  period  relative  to  t 
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increase  and  the  disequilibrium  will  become  more  important. 

Let  us  assume  that  the  gap  follows  the  equilibrium  behaviour 
shown  by  a solid  line  in  Fig.  4.3(b)  and  then  estimate  the 
correction  due  to  the  disequilibrium.  When  the  magnitude  of 
tlie  gap  is  decreasing  (which  occurs  while  the  bridge  is  carry- 
ing current  in  the  forward  direction)  the  energies  of  the 
quasiparticles  will  be  reduced  as  shown  for  = L^/2  in 
Fig.  4.3(c).  In  the  absence  of  inelastic  scattering  and  dif- 
fusion, the  occupation  fraction  of  any  state  will  remain 
unchanged  and  therefore  will  be  smaller  than  the  equilibrium 
occupation  at  the  reduced  energy  as  shown  in  Fig.  4.3(d). 
Inelastic  scattering  and  diffusion  will  tend  to  restore 
equilibrium,  but  to  the  extent  that  there  is  a lag  in  this 
process,  as  in  Fig.  4.3(d),  (dashed  line  ~ Tj),  there  will 
be  fewer  excitations  than  in  the  equilibrium  model,  as  if  the 
bridge  were  cooler.  Tlius,  the  magnitude  of  the  order  parameter 
will  be  somewhat  larger,  and  will  lead  to  a larger  forward 
supercurrent  during  this  half  of  the  cycle.  On  the  return 
part  of  the  cycle  [see  Fig.  4.3(b)-(d)],  the  gap  is  increas- 
ing, the  occupation  numbers  are  larger  than  thermal,  the 
magnitude  of  the  order  parameter  is  decreased,  and  the  super- 
current (which  is  negative  in  this  part  of  the  cycle)  is 
therefore  less  negative.  Thus,  during  both  parts  of  the  cycle 
an  extra  positive  supercurrent  occurs,  compared  to  the  equi- 
librium model,  and  its  magnitude  should  be  proportional  to 
the  rate  of  change  of  the  gap,  and  hence  the  voltage.  Since 
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a current-biased  junction  at  low  voltages  already  carries  a 

current  ~ in  tlie  equilibrium  model  (because  most  of 

the  time  is  spent  near  the  maximum  of  the  forward  super- 

current  part  of  the  cycle),  this  additional  supercurrent 

proportional  to  the  voltage  will  result  in  a characteristic 

which  starts  from  I and  rises  with  a slope  smaller  than 

CO 

the  normal  resistance  of  the  link,  up  to  voltages  of  order 

h/2e"g  (a  few  pV  for  tin).  This  describes  the  observed  foot 

at  low  voltages,  except  one  would  still  expect  an  initial 

vertical  rise  of  the  voltage  at  I similar  to  the  RSJ  model. 

CO 

As  the  voltage  begins  to  exceed  h/2eTg,  the  quasi- 
particle occupations  cease  to  follow  closely  the  equilibrium 
values  and  there  should  be  a deficiency  of  quasiparticles  in 
those  states  below  throughout  the  entire  Josephson  cycle. 
This  suggests  that  the  enhancement  of  the  magnitude  of  the 
order  parameter  will  reach  a limiting  value  (corresponding 
to  no  relaxation  at  all),  and  the  supercurrent  will  be 
enhanced  in  both  the  forward  and  reverse  part  of  the  cycle. 
Thus  there  should  be  a transition  to  a regime  with  the  usual 
behaviour  of  a weak  link  (imperfectly  understood  though  that 
may  be)  but  with  an  enhanced  critical  current  parameter  1^^* 
Ihis  is  the  region  on  each  I-V  curve  which  appears  to  have  a 
critical  current  and  a value  of  the  apparent  excess 

supercurrent  I which  is  a fraction  of  I , rather  than  I 

s cl  CO 

The  key  to  a more  quantitative  treatment  of  the  situa- 


tion described  above  is  the  derivation  of  generalized 
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Ginzburg-Landau  equations  for  the  order  parameter  which 
contain  terms  depending  explicitly  on  the  quasiparticle 
distribution  function,  together  with  Boltzman  equations  for 
the  quasiparticle  distribution  function.  Tliis  approach  is 
adopted  both  by  Aslamasov  and  Larkin  and  by  Golub,  but  they 
make  rablier  different  physical  and  mathematical  approxima- 
tions. As  a result,  Golub's  theory  is  more  appropriate  to 
describe  tlie  time  variation  within  the  Josephson  cycle  at 
low  voltages,  while  that  of  Aslamazov  and  Larkin  is  more 
appropriate  to  describe  the  time-averaged  limiting  behaviour 

near  I , . 
cl 

4.3.2  Golub's  Theory  (-j  > Tg) 

Golub  describes  the  derivation  of  a modified  form  of  the 
Ginzburg-Landau  equation  appropriate  to  the  short  weak  link 
geometry,  and  a perturbation  scheme  for  solving  it  at  low 
voltages.  The  final  equation  has  the  form 


1^ 


(— 

T'at 


+ 


1 oil! 


Y 


+ (1  - |Y 


(4.2) 


where  the  order  parameter  Y = i-'5  normalized  to  its 

magnitude  in  the  banks,  and  the  times  are  normalized  to  the 
current  relaxation  time  t = t^./u  where  u = n^/14C(3)  = 

O VjXj  O O 

5.79  in  the  dirty  limit.  Hie  usual  TDGL  equation  would  cor- 
respond to  the  values  u„  = u and  u^  = 

i o 


0 


Instead,  Golub 
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shows  that  is  very  small  and  the  new  term  is  the  most 
important  one  for  typical  materials.  As  we  shall  see,  the 
smallness  of  is  related  to  the  fact  that  the  quasiparticle 
diffusion  length  rather  than  the  coherence  length  governs 
the  variation  of  quasiparticle  electrochemical  potentials, 
and  tlie  more  important  u term  arises  from  the  slow  relaxa- 

Xj 

tion  of  disequilibrium  by  the  mechanisms  we  discussed  above. 

Golub's  equation  is  derived  within  the  particular  micro- 
scopic picture  described  by  Schmid  and  Schttn^^  (which  Golub 
related  to  the  work  of  Larkin  and  Ovchinnikov”"^ ) . In  this 
picture,  the  distribution  function  f(E)  of  quasi-electron 
states  as  a function  of  energy  E (with  the  convention  that 
states  inside  the  Fermi  surface  have  negative  energies)  is 
given  by  the  equilibrium  Fermi  function  f^(E),  plus  odd  and 
oven  correction  terms  associated  with  the  longitudinal  and 
transverse  modes,  respectively.  The  longitudinal  correc- 
tion f adds  equal  numbers  of  electron-like  and  hole-like 

i-i 

excitations  respectively  above  and  below  the  Fermi  surface, 

and  is  related  to  changes  in  the  magnitude  of  the  order 

parameter.  The  transverse  correction  f^  is  associated  with 

changes  of  the  gauge  invariant  electrochemical  potential 

difference  p - u = $ = -^  (200  + ^)  * where  cp  is  the  phase 
p 2 d t 

of  the  order  parameter  and  0 is  the  scalar  potential;  it 
adds  quasi-electrons  both  above  and  below  the  Fermi  surface 
corresponding  to  a net  charge. 


66 


In  computations  the  distribution  function  f(E)  is 
weighted  by  appropriate  functions  N^(E),  N^CE),  and  ^2^^^ 
shown  in  Fig.  4.4  for  three  different  values  of  N^(E) 

is  the  normalized  density  of  states,  rounded  off  by  phonon 
broadening  which  leads  to  the  existence  of  a small  but  non- 
zero density  of  states  between  -A  and  A.  R2  (E)  looks  simi- 
lar to  the  density  of  states,  but  it  is  an  odd  function  of 

E and  falls  as  a/E  far  from  the  Fermi  surface.  For  large 

2 2 ^2 

"gA*  ^2^^^  ~ a/ (E  - A ) which  is  the  usual  weighting  factor 

for  the  quasi-particle  occupations  in  the  BCS  self-consistent 

gap  equation;  moreover,  this  form  is  precisely  N (E)  — 

X SA 

2 2 ^ 

where  E = (e  + A ) The  physical  significance  of  N2{E), 
which  is  important  only  in  the  region  of  "unphysical"  states 
between  -A  and  A,  is  less  clear,  but  it  appears  to  pick  out 
states  for  which  the  transverse  imbalance  equilibrates  rapidly 
with  the  superconducting  condensate. 

The  steady-state  Boltzmann  equation  for  f^  balances 
inelastic  electron-phonon  scattering,  spatial  diffusion,  and 
a transverse-mode  quasi-electron  sink  rate  2AN2f^.  In  situa- 
tions where  the  spatial  variation  of  f^  is  slow,  the  equation 
has  the  approximate  solution 


(4.3) 


which  is  like  the  distribution  function  for  a normal  metal 
except  for  the  missing  energy  range  between  -A  and  A where 
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the  sink  rate  is  large.  This  solution  also  approximately 
satisfies  the  condition  of  electron-neutrality 


§ 


J dE  N^(E)f^  . 
^00 


(4.4) 


The  generalized  GL  equation  which  depends  on  f^  is 


nP  / 2 

8T  ^2 

c V Bx 


pCp  BA  \ 
Bx  3x  / 


= A f dE  N^f^ 

mmCD 


(4.5) 


wnich  for  f^  given  by  (4.3)  can  be  shown  to  reduce  approxi- 
mately to 


div  ^ I (4.6) 

S /,q2  e 

1 P 

where  >.g  = (DTq)  ^ = (Dt^T^/a)  ^ » F (T)  is  the  branch-mixing 

quasi-particle  diffusion  length^^'^^  and  a is  the  normal 

state  conductivity.  With  = - -^  and  the  continuity 

equation  divC?  + "3  ) =0,  this  establishes  as  the  length 
s n U ^ 

scale  for  variations  of  f . 

Tlie  steady-state  Boltzmann  equation  for  f balances 

Li 

electron-phonon  scattering,  diffusion,  and  the  generation  of 

3f 

disequilibrium  at  a rate  R | ~ which  for  large  t„A  re- 

^ ^ O'^ 

duces  to  N I . ' — I -1^  . Tills  justifies  our  interpretation 
of  the  origin  of  the  disequilibrium  as  arising  from  the 
changes  of  the  quasiparticle  energies  as  the  gap  changes. 
Golub  assumes  that  this  disequilibrium  relaxes  by  spatial 
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diffusion  alone  (although  the  effect  of  phonon  scattering 
does  enter  through  the  broadened  density  of  states),  which 
leads  to  the  solution 


I 

I 


M. 

bt 


(4„7) 


where  L is  the  length  of  the  weak  link.  Thus  f , shown  in 

Xj 

Fig.  4.5  is  a function  localized  between  -A  and  A,  where  the 
coefficient  (N^  **  ^2^  diffusion  term  in  the  Boltziitann 

equation  is  small.  The  insets  of  Fig.  4.5  show,  in  exag- 
gerated form,  the  shape  of  the  total  quasi-electron  occupa- 

I tion  f(E)  = f + f (with  f_  = 0),  Increased  quasi-electron 

O 1 

occupation  at  negative  energies  corresponds  to  fewer  hole- 
like excitations  below  the  Fermi  surface.  Thus  we  see  that 
a decreasing  gap  decreases  the  number  of  electron-like  and 
hole-like  excitations,  as  suggested  qualitatively  above, 
and  an  increasing  gap  increases  the  number  of  excitations, 
compared  to  the  equilibrium  model.  The  resulting  corrections 
to  the  gap  are  obtained  directly  from  the  generalized  Ginzburg- 
Landau  equation  (neglecting  unimportant  static  terms) 


2 

/ b 

8T  I ^ 2 

c \ ax 


(^)' 

'ax' 


^ oA 

8T  at 
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00 

J dE 
_00 


(4.8) 


The  two  real  equations  (4.5)  and  (4.8)  can  be  combined 
into  a single  complex  equation  of  the  form  (4.2)  where  the 
u_  and  u,.  terms  arise  from  the  integrals  of  f and  f , using 

X i-i  X Li 


-1.5  Shape  of  the  lonp.i tudinal  correction  f.  to  the 
equilibrium  distribution  function  f as  a func- 
tion of  K/A.  The  inset  shows  the  Sverall  dis- 
tribution function  with  the  nonequilibrium  cor- 
rection exanr^erated  for  the  cases  dA/dt<0  and 
dA/dt>0. 
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ms  (4.^)  and  (4.7).  At  the  level  of  approximation 


ic  eads  to  Eq.  (4.6)  s is  simply  and  is 

sma] 1 because  the  quasi-particle  diffusion  length  = 

i-  1 

(D-^  ^ is  much  longer  than  the  coherence  length  § = (Dt^j^)  ^ 

so  that  the  quasi-particle  current  and  hence  the  supercurrent 

must  be  nearly  uniform  on  the  length  scale  of  the  bridge. 

The  otlier  term  u is  given  by  u,.  = (1  + y - 6)u  where 
J-i  Lj  o 

Y = I^(TyA)  (L^/D)t^L 


I = — 
Y 3-t 


2 J d(E/A)R2(N^  - Rp 


.2,-1 


(4.9) 


37ie  variation  of  1 with  t„A  is  shown  in  Fig.  4.6.  For  large 

Y ^ 

~„A,  I is  relatively  independent  of  7„,  and  the  diffusion 

£.  Y ti 

time  L^/d  plays  a major  role. 

Even  for  u^  = 0 it  is  difficult  to  obtain  an  exact 

solution  to  Eq.  (4.2)  but  Golub  has  obtained  a perturbation 

solution  which  is  valid  at  low  voltages.  He  assumes  that 

the  solution  to  (4.2)  is  of  the  form  Y = Y + Y,  where  Y is 

o 1 o 

given  by  Eq.  (4.1)  and  Yj^  is  a small  nonequilibrium  correction, 
Ihen  one  only  needs  to  solve 


2 

d Y 

.2  ^1 


X 

- L 


Y dcp 
o ^o 

„ I 2 dt 


sin  cp 


(4.10) 


with  boundary  conditions  Yj^(O)  = Yj^  (L)  = 0.  In  Fig.  4.7(a) 
we  show  the  real  versus  the  imaginary  part  of  the  order 
parameter  at  equal  time  intervals  for  ~ 
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Fia.  4.6  Variation  of  I with 

V 
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Figure  4.7.  Plot  of  the  complex  parameter  at  equal  time 
intervals;  each  line  corresponds  to  the 
values  of  the  order  parameter  along  the 
length  of  the  bridge.  a)  Equilibrium  case 
= 0;  b)  Nonequilibrium  case  for 

u^  = 0,  dcp  /dt  =5,  where  t is  measured 

in  units  of  t , the  current  relaxation  time. 

o' 
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assuming  that  '‘'^(L)  = e , 'i'^(O)  = 1 and  = constant 

(voltage-bias).  Each  line  corresponds  to  a projection  of 
the  spatial  variation  of  the  order  parameter  onto  the  com- 
plex plane  with  the  ends  of  each  line  corresponding  to  the 
values  at  the  ends  of  the  bridge.  In  this  case  the  super- 
current I = I sino/  is  forward  during  the  first  half  of 
the  cycle  and  reversed  during  the  second  half. 

The  magnitude  of  the  order  parameter  as  a function  of 
time  and  position  in  the  bridge  for  the  same  conditions  is 
shown  in  Fig.  4.8(a).  Each  curve  in  Fig.  4.8  corresponds  to 
those  in  Fig.  4.7,  but  now  the  solid  lines  correspond  to  the 
first  part  of  the  cycle  and  the  dashed  lines  to  the  second 
half  of  the  cycle.  In  Fig.  4.7(b)  and  4.8(b)  we  show  similar 
plots  but  now  u^  = 0 and  (dcp  /dt)u^  = 5 so  that  Y = Y + 

Now,  because  of  the  disequilibrium,  during  the  first 
part  of  the  cycle  the  order  parameter  (which  is  decreasing 
in  magnitude)  is  larger  in  the  middle  of  the  bridge  than  the 
u = 0 (equilibrium)  value.  At  the  middle  of  the  cycle  there 
is  little  difference  from  the  equilibrium  case  because  the 
magnitude  of  Y is  stationary  at  cp^  = tt.  But  in  the  second 
half  of  the  cycle  the  order  parameter  in  the  middle  of  the 
bridge  (which  is  now  increasing  in  value)  again  lags  and  is 
now  smaller  than  the  equilibrium  value.  Thus,  as  discussed 
qualitatively  earlier,  the  magnitude  of  the  supercurrent 
is  increased  during  the  forward  part  of  the  cycle  and 
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Figure  4.8.  Plot  of  the  magnitude  of  the  order  parameter 
as  a function  of  space  and  time.  The  solid 
lines  correspond  to  the  first  half  of  the 
cycle  with  the  order  parameter  decreasing 
and  the  dashed  lines  to  the  second  half  of  the 
cycle  when  the  order  parameter  increases. 

There  is  a one-to-one  correspondence  between 
these  curves  and  those  shown  in  Fig.  4.7. 
a)  Ut  = Ul  = 0;  b)  Ut  = 0»  ^cp^/dt)  u^^  = 5. 
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decreased  in  the  reverse  part  of  the  cycle  leading  to  a net 

enhancement  of  the  supercurrent  through  the  bridge.  The 

/dcp  \ 

correction  to  the  supercurrent  for  u^  = 0 and!  °)ut  = 5 

J-  \d  t / L 

is  shown  in  Fig,  4.9.  Golub's  analytic  expression  for  the 
supercurrent  is  only  approximate  in  that  at  a given  time  the 
actual  supercurrent  is  not  spatially  uniform  along  the  bridge 
as  required  by  Eq.  (4.8),  so  that  Golub's  solution  (like 
that  of  Likharev  and  Yakobson  for  u = u , u^  =0)  is  based 
on  the  spatial  average  of  the  supercurrent.  In  general,  the 
deviations  of  J from  spatial  uniformity  are  small  compared 
to  the  magnitude  of  the  enhancement.  For  completeness  the 
phase  as  a function  of  time  and  position  is  shown  in  Fig. 
4.10,  for  u^  = 0,  Uj^  = 0 and  (dcp^/dt)Uj^  = 5. 

If  one  uses  Golub's  expression  for  the  supercurrent,  one 
then  obtains  an  effective  conductance  at  lov/  voltages  which 
is  due  to  the  extra  conductance  provided  by  the  increasing 
supercurrent  as  the  voltage  is  increased.  This  effective 
resistance  is  given  by 


^ 

^EFF 


t2/b2 

A 

Is 


24 


(4.11) 


where  is  the  normal  resistance  of  the  link.  Although  this 
was  derived  under  the  voltage-bias  assumption,  we  assume  that 
in  the  current-biased  case  it  will  describe  the  differential 
resistance  of  the  foot  above  I . As  will  be  shown  later, 

CO 

for  the  materials  and  temperatures  of  interest  here  y » P 
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Figure  4.10.  Spatial  and  temporal  variation  of  the  phase 
for  a)  equilibrium  case  u = u = 0; 

Ij  1 

b)  u„  = 0,  ^cp  /dt)  u = 5.  In  the  second 

X O Xj 

plot  one  CL :ve  has  been  removed  in  the 
regime  where  the  approximation  » Y^^ 


breaks  down 
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so  that  at  low  voltages  characteristic  of  the  foot  the 
effective  resistance  is  given  by 


^FF 


1 + 


1 

(1  + y)u^lV?^ 

24 


(4.12) 


However,  this  resistance  is  that  expected  after  an  initial 
RSJ-like  rise  of  the  voltage.  Because  of  the  initial  assump- 
tion of  steady  state  Boltzmann  equations  for  the  quasi- 
particles at  each  point  of  the  Josephson  cycle,  and  because 
of  the  assumption  of  small  corrections  necessary  for  the 
perturbation  method  of  solution,  Golub's  results  should  only 
be  valid  at  small  voltages. 

4.3.3  Ihe  Theory  of  Aslamasov  and  Larkin  (tj  « t^) 

At  higher  voltages,  the  approach  taken  by  Aslamasov  and 
Larkin  is  more  appropriate.  They  write  the  nonequilibrium 
Ginzburg-Landau  equation  in  the  form 


ttD 

8T 


7a3) 

8n^T^ 

c 


f ■ (E)dE 


0 


(4.13) 


where  f = f-f  is  the  difference  between  the  nonequilibrium 
o 

distribution  function  and  the  Fermi  function.  This  assumes 
that  the  occupation  of  electron-like  and  hole-like  excita- 
tions is  equal,  i.e.. 


that  the  transverse  mode  is  irrelevant 


I 
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to  the  problem.  Note  that  with  neglect  of  certain  terms, 
this  equation  is  completely  equivalent  to  Eq,  (4,8).  In 
particular  A(E  - a )“  is  just  the  limiting  form  of  R2(E) 
for  Tj,A  -♦ 

The  real  difference  between  the  two  approaches  lies  in 
the  calculation  of  the  nonequilibrium  distribution  function 
of  the  quasi-particles,  Golub  retains  an  explicit  time 
dependence  of  the  order  parameter  and  integrates  over  energy, 
so  that  the  magnitude  (and  sign)  of  the  correction  to  the 
equilibrium  distribution  function  of  the  quasi-particles 
depend  explicitly  on  the  time  evolution  of  the  gap,  Aslamasov 
and  Larkin,  on  the  other  hand,  calculate  the  distribution 
function  using  an  approach  in  which  the  quasi-particle  dis- 
tribution in  the  weak  link  is  averaged  over  space  and  time 
very  early  in  the  calculation,  and  only  the  energy  depen- 
dence is  explicitly  retained.  This  leads  to  rather  different 
physical  and  mathematical  approximations,  and  is  justified 
only  when  the  Josephson  period  is  much  smaller  than  the 
relaxation  times.  At  low  voltages  where  the  linear  foot 

is  observed,  the  inelastic  scattering  time  Tg  (of  the  order 

—9  —10 

of  10  to  10  seconds  for  tin)  is  comparable  to  the 
Josephson  period.  Thus,  the  low  voltage  predictions  of 
Aslamazov  and  Larkin  should  not  be  applicable.  At  higher 
voltages  they  calculate  a time  averaged  distribution  func- 
tion for  the  excitations  of  the  form 
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f (E)  = ijl  - j (4.14) 

for  E < A . At  energies  above  A » they  assume  that  the  dis- 
o o 

tribution  is  thermal,  because  of  the  rapid  diffusion.  This 
function  corresponds  to  a smaller  than  equilibrium  occupa- 
tion number  for  those  excitations  with  energies  less  than 
A^,  just  as  was  the  case  for  Golub's  theory  during  the  most 
important  (forward  current)  part  of  the  cycle.  Using  Eq. 
(4.14)  for  f,  they  obtain  that  the  last  term  in  Eq,  (4.13), 
which  they  denote  f (A ) , is  given  by 

/T  A^  A 

$(A)  = ^ ^ (1-  f-)  (4.15) 

When  this  and  the  gradient  terms  are  assumed  to  be  the 

most  important  ones,  the  order  parameter  variations  are 

T^-  T k 

characterized  by  a length  scale  v = E (T)  ( — — ) . if  the 

c 

length  of  the  bridge  L is  » r,  they  obtain  that  the  super- 
current is  enhanced  and  its  maximum  is  given  by 

= K - (4.16) 

^co  ^ 

where  K is  a numerical  coefficient  not  calculated  but  of  order 
one,  Ihis  result,  unfortunately,  is  derived  under  the 
rather  stringent  assumptions  that  5 » ^ » 11/  which  can  be 
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satisfied  only  extremely  near  to  since  % and  r differ 

only  by  a factor  of  three  even  at  0.99  T^;  nevertheless,  it 

may  be  useful  enough  for  comparison  with  experimental  results. 

For  currents  larger  than  I , Aslamazov  and  Larkin  calculate 

cl 

the  shape  of  the  I-V  curve  based  on  certain  assumptions 
about  the  current-biased  nature  of  the  weak  link  near  the 
forward  current  portion  of  tlie  Josephson  cycle.  The  result 
they  obtain  is  that  the  I-V  relation  should  be  given  by 

Since  this  depends  on  the  pulse-like  nature  of  the  response 
of  the  current-biased  junction,  it  is  presumably  valid  only 
for  currents  which  do  not  greatly  exceed 

4.4  Discussion  of  Experimental  Results 

Together,  the  Golub  and  Aslamazov  and  Larkin  models, 
provide  a semi-quantitatively  satisfactory  explanation  for 
the  I-V  curves  that  we  observe. 

In  Fig.  4.11  we  show  Golub's  prediction  Eq.  (4.11)  of 
the  normalized  resistance  of  the  foot  region  as  a function 
of  temperature  for  several  values  of  the  bridge  length  para- 
meter and  two  values  of  the  inelastic-scattering  time  t^. 

In  generating  these  theoretical  curves,  we  have  assumed  an 
average  value  of  A throughout  the  cycle  of  0.8  a value 

Tj,  = 8 X lo”^^  and  2 x lo”^*^  seconds,  and  a coherence  length 


J 
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Fip;.  4.11  Predictions  of  the  model  of  Golub  for  the  norma- 
lized resistancf3  of  the  foot  as  a function  of 
temnerature  for  different  values  of  the  bridtJe^Q 
length.  The  solid  linos  corres’^ond  to  i^=8xl0~ 
seconds  and  the  dashed  linos  to  2xl0~10  seconds. 
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? (0)  = 0,13  am  appropriate  to  our  moderately  clean  films 

(mean  free  path  0,1  um).  The  calculations  are  relatively 

insensitive  to  the  value  of  because  for  temperatures  less 

tiian  0.99  T , is  at  least  100,  so  that  I does  not  de- 

pend  strongly  on  it.  For  bridge  lengths  0,5-0. 8 am  and 

temperatures  of  0.9-0.98  T^,  the  parameter  v in  Eq.  (4.12) 

is  typically  10-20,  while  u^  = 5.79  so  that  the  corrections 

to  can  be  substantial, 

N 

In  Fig.  4.12  we  show  a plot  of  the  normalized  resistance 
of  the  foot  region  for  six  different  bridges  of  similar 
geometries  but  varying  resistances  compared  to  the  theoreti- 
cal curves  for  L = 0.7-0. 9 am  and  7„  = 8 x 10~^*^  seconds 
corresponding  to  that  determined  experimentally  by  Skocpol 
et  al.  Although  SEM  pictures  (shown  in  Chapter  Two)  of 
our  samples  show  that  the  distance  between  the  edges  of  the 
banks  is  usually  0.2  am  smaller  than  that,  the  appropriate 
length  for  comparison  with  the  theory  is  not  well  defined 
and  may  wall  be  different  than  that  distance.  Near  T^,  the 
experimental  effective  resistances  are  somewhat  higher  than 
predicted  by  theory,  Tliis  may  occur  because  the  temperature 
dependence  of  the  theory  is  not  quite  right,  or  it  may  result 
from  the  difference  in  T^,  discussed  in  Chapter  Tnree,  between 
the  bridge  and  the  lianks,  so  that  additional  contributions 
of  the  proximity  effect  affect  the  results.  For  T^  we  have 
used  the  temperature  at  which  the  linear  temperature  depen- 
dence of  the  critical  current  extrapolates  to  zero,  which 
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Fi'’.  4.12  Tho  norMally.ef.  rosintancc  of  tho  foot  vorsiis  temp- 
erature. The  data  points  are  for  six  bridr;es  and 
the  dashed  curves  are  Golub's  predictions  for  var- 
ious L values. 
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may  differ  by  10-20  millidegrees  from  the  temperature  at 
wbich  the  critical  current  actually  appears  to  go  to  zero 
and  by  10-20  millidegrees  from  the  temperature  at  which 
the  additional  resistance  of  the  entire  banks  is  observed „ 

One  feature  of  the  prediction  which  is  difficult  to 

4 

verify  but  appears  somewhat  suspect  is  the  strong  L length 

2 

dependence  (arising  from  tJie  L /D  in  y togetlier  with  the 
2 2 

additional  L arising  from  the  calculation  of  the  enhanced 

order  parameter).  Our  fabrication  method  is  unable  to  make 
bridges  with  substantially  shorter  geometrical  lengths,  and 
attempts  to  make  longer  bridges  by  moving  the  knife  side- 
ways have  been  inconclusive  because  the  first  cut  appears  to 
always  leave  a noticeably  deeper  cut  than  the  scraped  elonga- 
tion. Such  bridges  have  normalized  differential  resistance 
values  comparable  to  those  of  unelongated  bridges.  We  have 
also  varied  the  mean  free  path  of  our  films,  and  found  that 
the  normalized  differential  resistance  of  moderately  dirty 
films  also  does  not  differ  decisively  from  the  range  observed 
in  our  moderately  clean  films.  In  very  dirty  films  hysteresis 
was  present  very  close  to  and  the  foot  could  not  be 
observed  over  a wide  enough  temperature  range.  Since  the 
theoretical  assumption  of  a well-defined  length  witlu  rigid 
banks  is  obviously  an  idealization,  the  extent  to  which 
quantitative  agreement  with  the  normalized  differential 
resistance  and  its  length  dependence  should  be  expected  is 
not  clear. 
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Golub's  approach  does  mucli  better  than  the  expression 
given  by  Aslamazov  and  Larkin  for  the  slope  of  the  foot, 
calculated  assurring  » t Tlie  latter  is  several  orders 

ti  J 

of  magnitude  smaller  than  the  datao  Since  t is  comparable 

ij 

to  Tg  throughout  the  foot  region  in  our  bridges,  that  expres- 
sion should  not  be  expected  to  apply,  although  the  extent  of 
the  disagreement  is  somewhat  surprising  and  may  indicate 
some  difficulty  in  that  particular  calculation  by  A-L. 

At  voltages  above  this  initial  region,  the  Aslamazov 
and  Larkin  solution  is  more  appropriate.  According  to  their 

model  the  super -current  is  enhanced  up  to  a current  I , which 

cl 

is  of  order  (L/-^)  In  Fig.  4,13,  we  plot  the  ratio 

as  a function  of  temperature  for  the  same  six  samples 
used  in  Fig.  4.12.  Near  T^,  the  ratio  grows  very  rapidly  in 
a ■nanncr  consistent  with  the  temperature  dependence 
^r-l  H 

« (1  - T/T^)  . At  lower  temperatures  the  ratio  seems  to 

CO 

saturate  and  remain  constant  rather  than  continuing  to  grow 
as  predicted  by  theory.  At  much  lower  temperatures,  heating 
effects  start  playing  an  important  role  and  this  ratio 
decreases  until  the  two  currents  are  the  same.  Tnis  is  shown 
in  Fig.  4.14  where  the  currents  I , and  I are  shown  for  one 
bridge  over  the  whole  temperature  range  where  a foot  can  be 
observed.  Tne  dashed  lines  in  Fig.  4.13  correspond  to  Eq. 
(4.16).  TJie  appropriate  values  for  KL  vary  from  0.5  qm  to 
0.65  urn.  For  K=0.8  the  range  and  systematic  variation  of  the 
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4.13  The  ratio  of  the  cnlianceci  suoercurrent  I to 
the  equilibrium  critical  current  I as  S func 
tion  of  teinnorature  for  the  same  b??dp:es  as  in 
Fir,.  4.12.  The  dashed  curves  corresnond  to  the 
Aslarnazov  and  T.nrf-.in  nrediction  for  different 
values  of  KL  where  I'.'  1. 
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effective  length  from  sample  to  sample  is  in  agreement  with 
the  values  determined  in  Fig.  4.12.  Tlie  saturation  observed 
at  low  temperatures  for  this  ratio  is  not  contained  in  the 
theory.  As  previously  mentioned,  the  theory  was  derived 
under  tlie  conditions  r)  « L « 5 (T)  '<  both  of  these  conditions 
are  only  marginally  satisfied  in  the  range  shown  in  Fig.  4.13. 
At  most  the  lengtli  of  the  bridges  is  two  times  the  nonequi- 
librium  length  n,  and  at  low  temperatures,  the  coherence 
lengtli  becomes  shorter  than  the  length  of  the  bridge.  In 
Fig.  4.14  the  regime  r,  < L < is  indicated  by  the  dashed 

lines.  The  order  of  magnitude  of  the  effect  does  agree  well 
with  that  predicted  by  the  A-L  theory. 

For  a limited  range  of  currents  above  AL's  pre- 

diction Eq.  (4.17)  for  the  shape  of  the  I-V  characteristics 
can  be  fit  reasonably  well,  although  the  test  is  far  from 
stringent.  Fig.  4.1  includes  three  theoretical  curves  of 
this  form  fit  with  the  value  of  and  adjusting  and 

to  an  appropriate  value  for  each  curve. 

For  still  higher  voltages  the  theoretical  situation  is 

not  at  all  clear.  This  is  tiie  region  of  apparent  excess 

— 55-57 

current  which  is  a feature  of  both  the  phase-slip  model 

58  59  — 

and  the  TDGL  theory,  ' Both  predict  that  should  be  a 
definite  fraction  of  I^.  Presumably  in  the  nonequilibrium 
situation  described  here  t}ie  relevant  critical  current  is 

the  effective  critical  current  for  the  a'-tual  operating 
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conditions.  In  Fig.  4.15,  wo  plot  I , I ,,  and  I for  a 

CO  cl  s 

typical  bridge.  Near  T , I is  about  0.6  I while  at  lower 

c S CO 

temperatures,  it  is  about  0.7  qualitative  change- 

over to  tlie  dependence  on  I , is  clear-cut,  and  within  this 

cl 

interpretation  tends  to  support  I , as  the  effective  criti- 

cl 

cal  current  under  the  operating  conditions.  Note  that  I , 

CO 

not  the  linear  temperature  dependence  extrapolating 

to  T expected  for  the  zero  dissipation  critical  current, 
c 

Tne  behaviour  described  is  typical  of  all  the  microbridges 

and  in  Ibble  4.1,  we  include  the  values  of  at  low  tempora- 

i tures  as  well  as  the  maximum  value  of  I ,/I  observed  in 

' cl  CO 

each  bridge.  Note  that  the  value  of  vary 

much  from  bridge  to  bridge.  The  same  is  not  true  of  I ,/I 

cl  CO 

as  bridges  with  equivalent  resistances  do  show  substantial 
differences  from  a low  of  1.2  to  a high  of  2.0. 

Some  comments  should  be  made  about  the  material  depen- 
dence of  the  effects  discussed.  The  characteristic  feature 
of  interest  has  been  observed  only  in  microbridges  made  of 
tin  and  indium.  Since  these  materials  are  very  similar 
the  theory  should  describe  both  equally  well.  On  the  other 

hand  the  I-V  characteristics  of  short  aluminum  microbridges 

28 

of  lengths  simila  co  ours  do  not  show  similar  effects. 

We  have  calculatjd  the  effective  resistance  of  the  low  volt- 
age region  for  aluminum  microbridgos  of  similar  lengths  and 
find  that  the  change  from  the  normal  resistance  would  be 
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Flp,.  4.15  Plots  of  the  equilibrium  critical  current  Icot  the 

excess  sunercurrent  1..,  the  enhanced  supercurrent 

I , as  a function  of  tomoerature  for  bridge  43. 
cl 
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Ttible  4.1 

Foot  Parameters  for  Some  Tin  Microbridges 


Number 

I ,/I 

—cl'  —CO 

15 

0.06 

2.00 

0.77 

11 

0.33 

1.84 

0.72 

13 

0.09 

1.67 

0.78 

9 

0.36 

1.70 

0.72 

14 

0.20 

1.62 

0.79 

18 

0.13 

1.62 

0.72 

19 

0.06 

1.48 

0.75 

8 

0.44 

1.38 

0.70 

21 

0.07 

1.34 

0.74 

7 

0.43 

1.28 

0.78 

10 

0.85 

1.19 

0.65 

22 

0.62 

1.20 

0.74 
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much  smaller  than  that  predicted  for  tin.  This  is  due  to 

<j 

the  very  long  coherence  length  of  aluminum  which  v/ill  tend  to 

oppose  changes  in  the  order  parameter  different  than  the 

equilibrium  case.  In  order  to  observe  the  supercurrent 

enhancement  in  aluminum,  one  would  require  microbridges  with 

L/ij  ratios  comparable  to  those  of  our  tin  microbridges. 

However,  because  of  the  longer  inelastic  scattering  time  of 
66 

aluminum  the  enhancement  would  appear  at  very  low  voltages 
where  it  might  not  be  practical  to  observe  it. 

In  the  case  of  lead,  which  has  a shorter  scattering 
time,  the  transverse  term  might  become  comparable  to  the 
longitudinal  one,  but  it  is  not  presently  feasible  to  fabri- 
cate bridges  satisfying  the  condition  L < g (T)  because  of 
the  shorter  coherence  length  in  lead.  The  I-V  character- 
istics of  a lead  microbridge  comparable  in  size  to  our  tin 
microbridges  is  shown  in  Fig.  4.16.  Tlie  I-V  characteristics 
are  quite  different  than  those  of  tin.  Near  T^,  tiie  curves 
show  a change  to  the  resistive  state  above  I^^  marked  by  a 
steep  rise  and  a curve  without  inflection  points,  very 
similar  to  that  observed  in  tin.  However,  as  the  temperature 
is  lovered  the  characteristics  are  quite  different.  Tiie 
features  of  the  gap  and  subharmonic  gap  structure  are  not 
visible  on  the  I-V,  though  a feature  at  the  gap  is  seen  if 
dl/dV  is  measured.  At  low  voltages,  where  one  would  expect 
the  foot  to  be  observed  {j ^ is  of  order  10~^^  seconds,  so 
nonequilibrium  effects  should  become  important  at  voltages 
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of  order  50  uV) , instead  of  a foot  one  observes  a sharp 
rise  up  to  a voltage  of  about  20  uV,  followed  by  a region 
witli  an  excess  supercurrent  of  order  0.8,  Tne  region 
below  20  luV  becomes  hysteretic  very  rapidly,  and  at  lower 
temperatures  the  hysteresis  becomes  v/ider.  These  character- 
istics are  not  understood  at  present  and  it  is  clear  that 
further  work  with  high  materials  will  bo  required  before 
a coherent  understanding  of  the  properties  of  these  micro- 
bridges can  be  achieved.  It  might  be  that  the  I-V  character- 
istics observed  are  due  to  flux-flow.  Waile  we  believe  our 
tin  microbridges  are  too  narrow  to  contain  a single  vortex, 
the  same  is  not  true  of  our  lead  bridges.  Ho'vever,  these 

I-V  characteristics  are  different  than  those  predicted  by 
67 

Likharev  for  flux-flow  through  a microbridge. 


aiAPTER  FIVE 


HIGH-FREQUENCY  PROPERTIES  OF  VARIABLE-THICKNESS  MICROBRIDGES 

The  extreme  sensitivity  of  Josephson  devices  to  electro- 
magnetic radiation  suggests  the  possible  at-plication  of  these 
devices  as  detectors  in  the  microwave  and  far-infrared 
regions  of  the  spectrum^  Among  Josephson  devices,  thin  film 
microbridges  are  especially  attractive  for  practical  applica- 
tions because  of  their  compatibility  with  thin-film  technol- 
ogy for  the  fabrication  of  complex  circuits  as  well  as  their 
small  capacitance.  However,  the  performance  at  high  fre- 
quencies in  two-dimensional  thin-film  microbridges  is  con- 
siderably inferior  to  that  of  point  contact  devices.  Skoepol 

4 

et  al  had  shown,  when  this  work  was  initiated,  that  heating 
effects  impose  a limit  on  the  voltage  (and  thus  frequency) 
at  which  Josephson  effects  could  be  observed.  In  that 
analysis,  the  superiority  of  the  point  contact  was  suggested 
to  be  due  to  its  three-dimensional  geometry  which  more 
efficiently  removes  the  heat  generated  in  the  weak  link 
region.  Therefore,  the  variable-thickness  geometry  should 
provide  a system  whose  performance  is  closer  to  that  of  point 
contacts,  since  it  should  remove  the  heat,  if  not  as  effic- 
iently as  the  point  contact,  at  least  more  efficiently  than 
in  two-dimensional  uniform  thickness  microbridges , In  this 
chapter  we  describe  the  response  of  our  variable-thickness 
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microbridges  to  microwave  radiation.  It  is  shown  that  the 
dramatic  improvements  in  performance  observed  are  due  to 
their  reduced  self-heating.  Furthermore,  it  is  found  that 
heating  modifies  the  ac  response  of  these  devices,  as  well 
as  determining  their  maximum  frequency  of  operation.  Heat- 
ing effects  are  found  to  be  sufficiently  reduced  in  the 
variable-thickness  geometry  that  they  represent  a promising 
configuration  for  microwave  and  far-infrared  applications. 

In  Section  5.1  the  experimental  ac  Josephson  properties 
of  our  microbridges  are  presented.  At  low  temperatures  we 
observe  ac  Josephson  steps  to  voltages  as  high  as  3.7  mV  in 
our  best  samples.  This  corresponds  to  frequencies  of  2 THz 
and  is  a dramatic  improvement  over  other  thin  film  micro- 
bridges. Our  experimental  results  led  to  a detailed  heat- 
ing model  for  high-frequency  three-dimensional  Josephson 
devices,  which  is  presented  in  Section  5.2.  The  model  pre- 
dicts an  exponential  decrease  of  the  critical  current  I ^ 

CO 

-P/P 

e where  P is  the  power  dissipated  in  the  weak  link 

region  and  P^  is  a characteristic  power  level  which  is  tem- 
perature and  material  dependent.  We  find  these  predictions 
of  the  model  to  be  in  excellent  agreement  with  our  experi- 
mental observations  in  tin  and  lead  variable-thickness  micro- 
bridges, as  well  as  with  those  reported  by  other  workers 

using  different  materials.  At  low  voltages  and  near  T 

c 

the  predictions  of  the  model  are  found  to  be  in  good  quali- 
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tativc  agreement  with  tlic  observed  cut-off  of  the  oscilla- 
tory Bessel  function  power  dependence  of  the  ac  steps  as  well 
as  with  the  temperature  dependence  of  their  amplitude.  This 
is  discussed  in  section  5.2,5.  From  the  model  it  is  also 
possible  to  establish  optimum  design  criteria  for  variable-thick- 
ness microbridges  to  be  useful  at  high  frequencies.  This 
is  discussed  in  section  5.2,6,  where  some  preliminary  results 
on  the  application  of  these  criteria  are  presented.  Finally, 
experimental  evidence  is  presented  in  section  5.3  which 
clearly  demonstrates  the  connection  between  performance, 
hysteresis,  and  heating.  Here  gap  and  subharmonic  gap  struc- 
ture is  used  as  a local  thermometer  of  the  temperature  rise 
induced  in  the  weak  link  by  dissipation.  The  values  of  the 
gap  obtained  from  the  subharmonic  gap  structure  extrapolated 
to  zero  dissipation  are  found  to  be  in  excellent  agreement 
with  the  BCS  theory. 


5.1  Ihe  ac  Josephson  Effect  and  Variable-Thickness  Microbridges 
For  the  purposes  of  describing  the  ac  response  of  our 
microbridges,  it  is  useful  to  introduce  the  resistively- 
shunted  junction  (RSJ)  model  in  which  the  total  current  flow- 
ing through  the  device  is  assumed  to  flow  through  an  ideal 
Josephson  element  sin  cp^  shunted  by  a normal  resistance 

and  a capacitance  C: 


I 


I 

CO 


sin  CO 
^o 


C 


dV 

dt 


(5.1) 
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whore  co  is  the  total  phase  difference  across  the  weak  link, 
^o 

Ihe  phase  difference  grows  in  time  with  voltage  according  to 
the  Josephson  relation 


dcD 

dt 


2eV 


(5.2) 


Using  Eq.  (5.1)  in  (5.2)  yields  an  oscillatory  supercurrent 
component  flowing  through  the  weak  link.  While  it  is  pos- 
sible to  detect  directly  the  radiation  emitted  by  the  junc- 
tion at  uy^  = 2eV/h,^® ' is  simpler  to  mix  the  internal 

oscillation  of  the  junction  with  an  external  radiation  field. 
Then  the  total  voltage  induced  across  the  weak  link  will  be 
V = + V^^cos(x^^  t) (assuming  an  ideal  voltage  bias). 

When  this  is  combined  with  Eq.  (5.2),  it  yields  a supercurrent 

I 

component  of  the  form 


I = I 


CO 


"n(^)  (=^0  " Ve)  " 


(5.3) 


where  is  an  arbitrary  phase  factor,  and  is  the  Bessel 
function  of  order  n.  Equation  (5.3)  implies  that  whenever 
the  Josephson  oscillation  frequency  equals  a harmonic  of  the 
frequency  of  the  applied  radiation,  a dc  beat  will  be  observed 
in  the  I-V  characteristics.  The  amplitude  of  this  dc  current 
beat  is  given  by 


J 
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2eV 


I = I J 
n CO  n 


rf 


'tTx 


(5.4) 


rf 


In  practice,  however,  the  junction  is  usually  current- 
biased  so  that  the  observed  behaviour  differs  from  Eq.  (5.4), 
and  the  ac-induced  stops  in  the  I-V  characteristics  have 
amplitudes  proportional  to  distorted  Bessel-f unctions  as 
first  determined  by  Russer.^*^  The  boundary  between  these 
two  different  behaviours  is  determined  by  the  product 

which  allows  the  definition  of  a normalized  frequency 
^•‘^r  f 

C = •^■"7 — ^ . For  n > 1 the  resistive  shunting  will  domi- 
^■*■00  N 

nate  and  the  step  amplitudes  will  be  given  quite  accurately 
by  Eq.  (5.4) „ For  0 < 1 most  of  the  ac  current  will  flow 
through  the  Josephson  element  and  the  step  amplitudes  will 
correspond  to  distorted  Bessel  functions.  In  a similar 
manner  Eq.  (5.1)  defines  a characteristic  frequency  l/R^^C 
above  which  the  capacitive  element  dominates.  However,  this 
should  be  of  little  practical  importance  in  variable-thickness 
microbridges  because  of  their  low  capacitance  estimated  to 
be  of  the  order  10”^-10~^  pF. 

Ihe  behaviour  described  by  this  model  is  readily  observed 
in  our  variable-thickness  microbridges.  In  Fig.  5.1  the 
response  to  32  GHz  microwave  radiation  of  one  of  our  tin 
microbridges  is  shown.  At  small  microwave  powers  a step 
appears  at  66  pV  corresponding  to  the  first  harmonic  of  the 
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Fig.  5.1  Response  of  bridge  #3  to  32  GHz  radiation.  The 
responds  to  66  mV. 
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applied  radiation.  As  the  microwave  power  is  increased 
further,  tlie  critical  current  I , and  the  first  and  the 

CO 

second  step  seem  to  behave  in  a manner  qualitatively  con- 
sistent with  Eq„  (5„4).  However,  note  the  presence  of  steps 
at  voltages  corresponding  to  subharmonics  of  the  fundamen- 
tal frequency  ^ x^^,  where  m = 1,2,...  . These  steps 

are  not  predicted  by  either  Eq,  (5,4)  or  the  current-biased 
RSJ  model,  but  are  commonly  observed  in  microbridges 

at  a variety  of  temperature  and  power  levels.  Their  obser- 
vation is  probably  related  to  the  current-phase  relation  not 
being  perfectly  sinusoidal  but  rather  being  slightly  anhar- 
monic.  While  measurements  of  the  current-phase  relation  of 
variable-thickness  microbridges  do  show  an  almost  ideal 
sinusoidal  relation, these  experiments  were  all  per- 
formed very  close  to  T^  and  at  zero  voltages.  At  finite 
voltages,  the  supercurrent  enhancement  discussed  in  Chapter 
Four  loads  to  an  anharmonic  current-phase  relationship  which 
would  explain  the  observed  subharmonic  steps. 

As  the  microwave  power  is  increased  further,  steps 
appear  at  higher  and  higher  voltages  while  the  amplitudes  of 
the  ac  steps  at  low  voltages  change  periodically  in  a way 
qualitatively  similar  to  Eq.  (5.4).  Figure  5.2  shows  I /I 

n CO 

for  the  first  ten  stops  induced  by  10  GHz  radiation  for  a 
typical  microbridge.  While  qualitatively  this  is  similar 
to  a Bessel  function  dependence,  no  amount  of  fitting  of 
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Fitr.  5.2  Full  step  widths  normalized  to  the  critical 

current  I for  the  first  10  steps  at  10  GHz. 
CO 
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either  the  current-biased  or  voltage-biased  models  yields 

detailed  agreement  with  tlie  data.  Typically,  the  step 

heights  observed  are  smaller  than  predicted  by  either  of 

tliese  models  and  the  periodicity  is  hard  to  observe  for 

high  powers  or  the  high  order  steps.  Near  the  width  of 

the  low-order  ac  steps  is  found  to  decrease  as  one  approaches 

in  disagreement  with  what  the  RSJ  model  predicts  for  fl  > 1. 

Tliese  discrepancies  can  be  explained  qualitatively  in  terms 

of  heating  effects  as  discussed  in  Section  5.2.5. 

At  very  high  pov/er  levels,  steps  are  observed  to  very 

high  voltages,  until  a maximum  voltage  is  reached.  Any 

further  increase  of  the  power  yields  no  steps  at  higher 

voltages  and  the  amplitudes  of  tliose  observed  below  V 

max 

begin  to  decrease.  Such  behaviour  is  not  expected  from 
Eq.  (5.4)  which  predicts  only  a smooth  decrease  in  the  ampli- 
tudes of  the  steps  with  no  sharp  cut-off.  This  V is  a 

max 

useful  index  of  the  performance  at  microwave  frequencies,  as 

it  corresponds  to  a maximum  observable  internal  frequency 

under  the  operating  conditions.  We  interpret  this  V as 

max 

the  voltage  at  which  phase  coherence  across  the  microbridge 
is  lost  due  to  the  growth  of  a normal  hotspot  because  of  t]ie 
increased  power  dissipation.  Tnis  interpretation  will  be 
appropriate  only  as  long  as  the  frequency  of  the  applied 
radiation  is  well  below  the  gap  frequency  at  the  operating 
temperature.  At  higher  fundamental  frequencies,  effects 


Ill 


due  to  the  Riedel  peak  and  subsequent  roll-off  in  the 

intrinsic  fre^iucncy  dependence  of  the  supercurrent  might  also 

71  7 3 

be  of  importance.  ' At  microwave  frequencies  and  away 

from  these  roll-off  and  peak  effects  can  be  neglected 

in  our  microbridges  because  fix  « 2 A so  that  V should 

max 

only  be  limited  by  heating  effects. 

Tile  interpretation  of  reflecting  a heating  limit 

in  superconducting  microbridges  has  been  confirmed  experi- 
mentally. Improvements  of  V as  a result  of  using  the 

variable-thickness  geometry  have  been  reported  by  Klapwijk 

. .74  49 

and  Mooij  and  Klapwijk  and  Veenstra.  Using  a mechanical 

scribing  machine  they  produced  aluminum  variable-thickness 

microbridges  with  V = 0.75  mV  at  T/T  = 0.93.  Tnis  rep- 

max  c ^ 

resented  a considerable  improvement  over  their  uniform 
thickness  microbridges  where  the  ac  steps  were  not  observed 
even  at  low  voltages.  Similarly  their  tin  variable-thickness 
microbridges  showed  stops  up  to  = 1.2  mV  at  very  low 

temperatures. 

Even  greater  improvements  are  observed  in  our  micro- 
bridges.  In  Table  5.1  we  present  data  from  a representative 
selection  of  our  microbridges.  Tlie  first  group  in  the  table 
consists  of  microbridges  with  low  impedances  and  fairly 
large  critical  currents.  Within  this  group  the  bank  thick- 
ness varies  from  0.1  qm  (a  uniform-thickness  bridge)  to  3,0 
qm  in  the  variable-thickness  geometry.  Hiere  is  a clearcut 
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l^blo  5.1 

Some  Cliaracter istic  Properties  of 
Representative  Tin  Microbridges 


R 

I (2.1  k) 

Banks 

Width 

V 

max 

(Q) 

''°(ma) 

(um) 

(urn) 

(mV) 

1 

0.15 

22.0 

0.1 

1.5 

0.30 

2 

0.08 

10.0 

1.0 

0. 5-1.0 

0.63 

3 

0.13 

12.5 

1.0 

0.5 

0.98"^ 

4 

0.24 

11.3 

3,0 

~ 1.0 

1.47'^ 

5 

0.27 

17.3 

3.0 

~ 1.0 

1.37'^ 

6 

0.40 

3.5 

0.1 

0.5 

0.70 

7 

0.43 

3.8 

1.0 

~ 0.5 

1.67 

8 

0.44 

7.5 

1.0 

< 0.5 

9 

0.36 

2.5 

1.0 

< 0.5 

2.27 

10 

0.85 

1.9 

1.0 

< 0.5 

2.64 

+ 

11 

0.33 

4.0 

1.7 

< 0.5 

3.70 

+ Indicates  that  observation  corresponded  to  maximum  power 
available  but  the  limit  for  the  bridge  had  not  been 
reached. 
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improvement  in  V for  the  thickest  banks,  of  almost  a 
max 

factor  of  five.  The  second  group  is  more  representative  of 

the  bridges  made  in  the  later  stages  of  our  ivork  during  which 

t]ic  fabrication  technique  and  the  burnout  problem  were  better 

controlled.  All  of  the  bridges  in  this  second  group  have 

higher  resistances,  smaller  critical  currents  and  narrower 

widths  (w  s 0.5  um) . Once  again  an  improvement  of  almost  a 

factor  of  five  from  the  uniform-thickness  bridge  with  0.1 

um  thick  banks  to  the  variable-thickness  bridge  with  thickest 

bank  and  highest  resistance  is  observed.  Within  each  group 

a considerable  difference  in  performance  between  high  and 

low  resistance  microbridges  is  observed.  Tliis  is  easily 

understood  since  the  observation  of  ac  steps  at  V requires 

max 

that  the  bias  voltage  and  the  microwave  voltage  be  of  order 
so  that  the  power  dissipated  in  the  weak  linx  scales 
as  giving  more  dissipation  in  the  low  impedance 

bridges , 

In  Table  5,2  v/e  show  the  values  of  V for  our  lead 

max 

var iable-tliickness  microbridges,  all  of  which  had  the  same 

bank  thickness  d = 1,5  pm.  The  voltage  limits  found  are 

not  too  different  from  those  in  tin;  although  the  highest 

value  of  V in  lead  is  belov/  that  in  tin,  the  number  of 
max 

lead  samples  studied  was  much  smaller.  Note  however,  that 
except  for  the  microbridge  with  = 10  ohms,  which  showed 
no  ac  steps,  the  bridge  with  the  best  performance  had  the 
highest  resistance. 
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Table  5.2 

Some  Characteristic  Properties  of  Lead  ^^icrobridges 


Number 

(0) 

I (4.2K) 
^°(mA) 

V 

max 

(mV^ 

L-6 

0.07 

16.8 

1.6 

L-4 

0.08 

37.5 

1.8 

L-5 

0.08 

20.5 

1.8 

L-8 

0.10 

4.1 

1.1 

L-2 

0.16 

33.6 

1.9 

L-3 

0.60 

3.5 

1.8 

L-1 

0.80 

1.9 

3.0 

L-7 

10.00 

0.042 

No  ste] 
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I 

I 

From  the  above  results,  one  obvious  possibility  is  to 

[ 

I use  thicker  and  thicker  banks  and  smaller  microbridges  in 

order  to  reduce  heating  effects  further.  However,  the  cut- 
ting of  ever  smaller  bridges  with  ever  thicker  banks  is 
increasingly  difficult.  It  was  found  that  with  a bank 
thickness  > 1.5  lum  the  yield  of  good  microbridges  was  sig- 
nificantly reduced  since  the  edge  of  the  diamond  knife  did 
not  penetrate  uniformly.  T:ius  a bank  thickness  of  1.5  ^j.m 
was  selected  as  a compromise  between  good  yield  and  good 
cooling.  I^ie  problem  of  increasing  bank  thickness  is  not 
restricted  to  our  crude  fabrication  technique.  In  more 

^ sophisticated  techniques  using  electron-resists,  the  bank 

thickness  is  limited  by  the  resist  tliickness  (~  1.0  um),  so 
that  any  further  increase  in  the  bank  tiiickness  requires 
the  use  of  many-stage  deposition. 

Thus,  the  observed  experimental  behaviour  is  consistent 
with  heating  determining  the  high  frequency  perfOx-mance  of 
thin  film  microbridges  in  general.  In  the  next  section 
further  implications  of  heating  are  explored  within  the  frame- 
work of  a simple  heating  model  for  the  case  of  three- 
dimensional  metallic  cooling  appropriate  to  the  variable- 
thickness geometry. 
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5.2  Tno  Heating  Model 
5.2.1  Intxoluction 

T:ie  increased  dissipation  as  the  voltage  and  the  rf 
pover  are  increased  produces  a gross  disequilibrium  inside 
tJia  small  weak  link.  One  can  approximate  tliis  situation  by 
defining  a local  temperature  distribution  inside  the  weak 

4 

link.  Earlier  analysis  of  this  problem  showed  that  in  a two- 

dinensional  system  should  scale  roughly  as  the  energy 

gap  2i/e  and  in  three  dimensional  systems  it  should  scale 

— t — 

approximately  as  (24/e) (^/a)  ^ where  § is  the  coherence 
length  averaged  over  the  temperature  distribution  inside  the 
link,  and  a is  the  neck  radius.  Tne  analysis  presented  in 
tliis  section  extends  the  previous  analysis  for  the  three- 
dimensional  geometry  by  using  the  collinear  cone  gecraetry 
intrO'duced  in  Cnapter  Hiree. 

In  the  presence  of  a current  I the  total  amount  of  power 

2 

dissipated  inside  the  weak  link  is  given  by  I R.  Tlie  power 
dissipation  will  not  all  occur  at  tlie  center  of  the  weak 
lin?<  but  rather  will  be  a function  of  radius  r.  Tlie  tempera- 
ture gradient  required  to  carry  tlie  power  away  from  the 
center  of  the  weak  link  is  given  by 


dT  P(_i^ 


dr 


« ^eff  ^ 


(5.5) 


where  is  the  thermal  conductivity  as  modified  by 


I 


j 
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boundary  scattering.  Tiiis  corresponds  to  a position  depen- 
dent resistivity  as  discussed  in  section  3.2.1, 

Tno  power  dissipated  inside  a radius  r is  given  by 


^ J 


eff  , 


(5.6) 


whore  we  have  assumed  the  validity  of  the  Wiedemann-Franz 
law  and  ~ ^ Jl/2r  sin6)  as  obtained  in  Chapter 

'Diroo.  Integrating  Eq.  (5,5)  from  r = <=  to  r we  obtain 


T^(r)- 


+ !eff 

2rT^<^nr  ' *^o 


(5.7) 

From  Eq,  (5.7)  the  maximum  temperature  in  the  middle  of  the 
weak  link  is  for  I « r^^sinG 


T = T(r-  ) 
m 1 


(5.8) 


Tliis  is  a very  general  result  valid  for  all  three-dimensional 
7 5 7 6 

geometries.  ' Furthermore,  all  material  properties  have 
dropped  out. 


Frotn  Eq,  (5.8)  one  can  estimate  the  temperature  rise 
produced  by  the  dc  voltage  and  the  microwave  radiation  in  our 
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variablc-tJiicknoss  microbr idgos , For  useful  response 

should  bo  of  order  so  that  the  total  power  dissipated 

2 ^ % 

is  V = (V  + 0.5  V c).  For  the  microbridge  with  the 
rms  dc  rf'  ^ 

best  performance  in  Table  5,1,  V = 3.7  mV  and  T,  = 2,1  K 

max  D 

which  yields  a maximum  temperature  in  the  weak  link  of  15  K. 
Thus  the  departures  from  equilibrium  are  fairly  large  at 
the  operating  voltages  of  interest.  From  Eq.  (5.8)  the 
voltage  necessary  to  raise  the  temperature  in  the  middle  of 


the  bridge  above  T^  is 


V - 0.312  T (1  ^ ) (mV) 

c 


(5.9) 


At  1^  = 0 a voltage  V ~ 1.2  mV  is  required  to  raise  the 

temperature  of  a tin  bridge  above  T^.  This  voltage  will  be 

larger  in  materials  with  high  T^  and  smaller  at  temperatures 

closer  to  T , 
c 

It  is  important  to  note  that  one  effect  of  these  large 
temperature  rises  in  the  weak  link  is  to  add  additional 
noise  due  to  the  hot  electrons  present  in  the  weak  link.  If 
one  assumes  the  weak  link  is  a series  of  Johnson  noise 
sources  it  can  be  shown  that  the  noise  temperature  is  given 
approximately  by 


T 

n 


1 

2 


(T  + T ) . 
b m 


(5.10) 
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5,2,2  Reduction  of  by  Heating 

n:e  Aslamazov  and  Larkin  model  introduced  in  section 
3,2,1  can  be  used  to  take  heating  into  account  for  the  cal- 
culation of  the  critical  current  in  the  presence  of  dissipa- 
tion, In  Eq,  (3,1)  the  order  parameter  in  the  weak  link  was 
found  to  be  of  the  form 


o 


icp 

fC?)  + (1  - f(?))e  ° 


(5,11) 


In  tlie  presence  of  a voltage  this  corresponds  to  the 
oscillating  order  parameter  shown  in  Chapter  Four,  In 
the  presence  of  dissipation  Y will  still  oscillate,  but, 
because  of  heating,  the  order  parameter  will  have  a maximum 
which  is  below  that  given  by  (5,11),  One  can  write  Eq, 
(4,11)  in  a more  general  form  to  include  heating 


Y 


f(7)  + f(-r) 


(5,12) 


Ihe  total  supercurrent  flowing  through  the  weak  link  is 

now 


J 

s 


f(-7)vf  (?) 


f (r)vf  (-?) 


sincp^ 


(5,13) 


v^ere  in  the  absence  of  heating  the  factor  in  parenthesis 
reduces  to  Vf,  In  the  middle  of  the  bridge  this  factor 
reduces  to  2f(0)  7f(0),  In  the  presence  of  heating  f(0) 
will  be  depressed  below  1/2  so  that  (y|  will  be  depressed 
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2 

to  (2f(0)'i^l  . Since  vf(0)/f(0)  remains  approximately  the 
same,  of  order  ~ (4/'Ta ) tan  (G/2 ) , in  the  collinear  cone 
geometry,  in  the  presence  of  heating  will  be  reduced 

in  proportion  to  (f(O)j  . 

In  tlie  case  of  large  dissipation  wlaich  is  of  interest 
here,  the  middle  of  the  weak  link  is  above  its  transition 
temperature  and  one  has  to  estinuite  the  ability  of  the  pairs 
to  penetrate  across  this  normal  barrier.  One  can  use  the 
Ginzburg-Landau  equation 

2 

- + ai'  + = 0 (5.14) 

2m 

In  the  normal  region  the  order  parameter  is  very  small  and 
the  nonlinear  term  in  Eq.  (5.14)  can  be  dropped.  The  radial 
equation  for  (5.14)  is 

(j-y)  = —J: (rY)  (5.15) 

dr^  r(r) 

where  now  the  coherence  length  Z (T)  is  a position  dependent 
quantity  because  of  the  large  temperature  gradients  at  either 
side  of  the  weak  link.  Then 

?‘'^(r)  = Z~^(0)  [t(r)  - 1]  (5.16) 

inside  of  the  normal  region.  In  order  to  estimate  how  Y is 
attenuated  in  this  normal  region  we  use  the  WKB  approximation 
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which,  from  the  loading  term,  gives  an  attenuation  factor 
for  rY 

2f(0)  = exp(-W)  (5.17) 

wliere 

r(Tc)  ^ 

W = J [t(r)-  1]^  dr  . (5.18) 

^1 

This  integral  can  be  written  in  terms  of  the  average  inte- 
grand as 

r(T^)-ri  / 

W = < [t(r)-  1]  ? (5.19) 

5(0)  \ y/AVE, 

Notice  that  Eq.  (5.19)  involves  simply  1/5  averaged  over  the 

temperature  profile  in  the  normal  region.  For  high  powers, 

r(T  ) is  given  by  Eq.  (5.7)  and  is  much  larger  than  r , a, 

C J. 

£ and  P so  that  all  terms  in  (5.7)  of  order  £/r (T  ) and 
o c 

rj^/r(T^)  can  be  neglected.  Furthermore,  the  average  in  Eq. 
(5.19)  can  be  replaced  to  a good  approximation  by  its 
value  in  the  middle  of  the  range  of  r 
r(T  ) , , A 

Prom  Eq . (5.7) 

9 . 3 e^p  P 

T^(r)  - = -5 — (5.21) 

^ TT^k  n r 
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from  which  tlie  value  of  r (T  ) is 

c 

. 2 

3 e Q p 

r(T)  = v-y  n ■ ■ y- 

" a T^d  -t^) 


(5.22) 


Since  t[j  r(T^)3  = (2  - t^)  so  that  r(T^)]-  1}  ^ = 


[■jd-  , then 


1 3 °o^ 

[2d-  ?(0)u 

D ^ 


(5.23) 


Since  is  expected  to  decrease  as  |2f(0)l  = exp(-2W),  then 

the  effect  of  heating  will  be  to  decrease  the  critical  cur- 
rent as 


I (P)  = I (P  = 0)  exp(-P/P  ) 
CO  CO  o 


(5.24) 


where 


P„(T)  = — (1-  t?)^"  K(T  )T  P(0)n 
o rr  b c c 


(5.25) 


P^(T)  is  a temperature-dependent  characteristic  pou-er  level 
which  is  geometry  and  material  dependent.  While  this  form 
was  obtainec  for  high  power  levels  it  can  be  shown^*^  that 
the  form  for  the  decrease  of  the  critical  current  at  low 
power  levels  is  essentially  the  same.  Hie  temperature  depen- 
dence in  Eq.  (5.25)  reflects  the  higher  power  required  to 
raise  the  temperature  of  the  weak  link  above  as  the 
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* 


temperature  is  reduced  below  T . The  material  dependence, 
which  scales  as  ? (0)/p^,  reflects  the  ability  of  the 
. material  to  remove  the  heat  as  well  as  its  ability  to  maintain 

phase  coherence  through  its  coherence  length.  Materials 
witli  high  resistivity  tend  to  have  short  coherence  lengths. 

This  can  be  balanced  by  a higher  so  that  at  the  same 
bath  temperature  more  power  will  be  required  for  a higher  T 

c 

bridge  to  be  driven  normal.  At  low  temperatures  assuming  a 

2 

residual  resistance  ratio  of  10,  F (0)/p^  is  approximately 
the  same  for  Al,  In,  Sn,  and  Nb.  It  is  much  smaller  in  higher 
[ T materials  because  their  short  coherence  lengths  and  high 

t 

1 

resistivity  more  than  offset  the  advantage  of  their  higher  T^'s. 

Assuming  t = 0 and  a cooling  solid  angle  of  tt  for  these  mater- 
i ^ 

I 

ials,  10  uW.  Thus,  from  this  characteristic  power  level, 

it  would  require  approximately  50  pW  in  order  to  reduce  the 

critical  current  to  1%  of  its  equilibrium  value. 

One  must  be  careful  in  assessing  these  results.  From 

our  assumptions  Nb  would  appear  to  be  favored  because  its 

^ higher  resistivity  would  yield  higher  resistance  microbridges 

with  equivalent  P . However,  while  this  would  be  true  of  bulk 

o 

‘ niobium,  it  is  apparently  not  possible  to  make  clean  thin 

films  of  niobium,  so  the  actual  values  of  P would  be  con- 

o 

siderably  below  10  pW. 


5.2.3  The  Variable-Thickness  Geometry 

In  the  preceeding  discussions  we  have  assumed  that  the 
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col linear  cone  geometry  is  a good  approximation  to  the 
variable-thickness  geometry,,  While  this  is  true,  it  is 
also  important  to  consider  the  case  of  a finite  bridge 
length  L which  is  attached  to  very  thick  banks  which  are 
approximated  by  the  collincar  cone  geometry.  In  this  manner 
it  is  possible  to  assess  the  importance  of  the  finite  bridge 
length  as  well  as  to  evaluate  the  performance  of  long  variable- 
tliickness  microbridges  where  the  simple  radial  cooling 
approximation  breaks  dO'wa.  For  very  large  power  levels  it 
is  possible  to  approximate  the  bridge  contribution  by  simply 
assuming  the  whole  bridge  is  at  the  maximum  temperature  T^. 

Tliis  establishes  an  upper  bound  to  the  additional  contribu- 
tion to  tlie  WKB  exponent  from  the  finite  length  of  tlie 
bridge.  Tiien  from  Eq.  (5.1?)  this  additional  contribution  is 


W . , 
oridge 


L__ 

2§  (0) 


(t  - 1)*^^ 
m 


(5.26) 


whore  t = T /T  . Tlie  critical  curr*.  nt  is  reduced  by 
m -ir  c ^ 


-2W  T-2W, 

I (P)  = I (P  .=  o)e  oridge 

CO  CO 


(5.27) 


wliere  given  by  Eq.  (5.23).  Note  that  the  bridge 

contribution  is  essentially  what  one  expects  for  the  reduc- 
tion of  the  critical  current  in  an  S-N-S  junction  with  co- 
herence length  ? (t  ) . 

m 
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At  high  voltages  eV/nkT^  ~ 5 for  our  tin  microbridges  and  2.5 
for  lead  so  the  bridge  contribution  will  be  important  if 
L ~ 12  rj^  in  tin  and  L > 5 r^^  in  lead.  While  in  tin  this 
contribution  should  not  be  important  (recall  r^^  ~ 0.04  um)  , 
in  lead  microbridges  and  other  high  materials  it  can  be 
important.  At  low  voltages  and  near  T^,  is  more 

important  than  given  by  Eq.  (5.30);  however,  the  nature  of 
the  approximations  made  and  the  uncertainty  of  the  length  L 
make  it  rather  difficult  to  make  quantitative  comparisons 


with  the  model 
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5.2.4  Voltage  Limit  for  Microwave  Induced  Steps:  Temperature 
and  Material  Dependence 

The  effect  of  dissipation  is  to  reduce  the  critical 

current;  this  in  turn  reduces  the  widths  of  the  ac  steps 

observed.  At  high  voltages  this  sets  a limit  for  the 

observability  of  the  ac  steps,  which  can  be  compared  to  our 

data.  In  Section  5.1  we  found  the  best  performance  in  a 

bridge  with  R,  = 0.33  ohm,  for  which  V = 3.7  mV.  Using 
N max  ^ 

Eq.  (5.8)  this  corresponds  to  T =15  K.  The  total  power 

m 

dissipated  in  the  bridge  is  62  qW,  so  that  if  ~ 10  gW 
the  critical  current  is  reduced  to  less  than  1%  of  its 
original  value.  This  is  sufficient  to  reduce  the  height  of 
the  step  below  detectability. 

A better  way  to  compare  the  model  with  the  data  is  to 
examine  the  power  dependence  of  the  ac  steps  in  detail.  In 
Fig.  5.3(c)  we  show  the  experimental  data  for  dl/dV  vs.  V 
for  the  same  bridge  with  V =3.7  mV.  Each  of  the  equally 
spaced  lines  corresponds  to  an  individual  step,  and  the  deriva- 
tive is  proportional  to  the  step  amplitude.  At  voltages  of 
order  2.5  mV  a striking  periodicity  develops  in  the  step 
amplitudes.  In  Fig.  5,3  (a)  we  plot  the  amplitude  of  the 
Bessel  function  for  the  case  = a>  (perfect  cooling).  The 
same  periodicity  observed  in  the  experimental  data  occurs 
at  high  n for  J^(x).  In  Figs.  5,3  (b),  (d)  and  (e)  we  plot 

the  same  Bessel  function  but  now  its  amplitude  is  affected 
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Figure  5.3.  Step  amplitudes  at  10  GHz  for  different  values 
of  compared  with  the  experimental  ampli- 
tudes for  tin  bridge  #11, 
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by  heating  as  is  decreased.  As  the  voltage  is  increased  the 
step  width  amplitudes  decrease  in  a manner  consistent  with  the 
data.  Hie  best  fit  for  appears  to  be  somewhere  between 
10-15  uW.  This  value  of  P^  is  in  very  good  agreement  with 
the  heating  model  presented  above,  given  the  many  approxi- 
mations and  assumptions  made. 

In  Section  5.1  we  discussed  the  power  dependence  of  the 
ac  steps  and  established  that  for  Q = (h!D/2eI^^R^)  > 1 the 
step  amplitudes  should  behave  like  Bessel  functions.  Since 
at  high  powers  is  considerably  reduced  below  its  P = 0 

value  it  is  valid  to  compare  the  amplitudes  of  the  ac  steps 
to  Bessel  functions.  At  lower  voltages  (hence  powers)  this 
breaks  down  as  the  normalized  frequency  Q < 1,  and  compari- 
son should  be  made  to  the  distorted  Bessel  functions  of 
Russer.^^  All  in  all,  the  agreement  at  high  voltages  is  a 
remarkable  confirmation  of  the  Bessel  function  dependence 
of  the  amplitude  of  the  ac  steps  to  the  highest  order 
reported  to  date. 

While  the  observed  data  is  in  very  good  agreement  with 
the  predictions  of  the  model,  it  is  important  to  consider  a 
criterion  for  the  observation  of  the  ac  steps.  Our  low 
impedance  microbridges  have  large  critical  currents.  However, 
as  the  dissipation  scales  as  V /R^  these  critical  currents 
decrease  significantly  with  power  and  the  ac  steps  become 
very  small.  A reasonable  assumption  based  on  the  work  of 
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Ambcgaokar  and  Halpcrin^^  and  Stephen^®  is  that  the  n^^  step 

in  the  absence  of  noise  should  exceed  ekT  /h,  where  T is 

N'  N 

tJie  noise  temperature.  Since  J (x)  oscillates  with  an 

n 

amplitude  ~ x until  it  reaches  a maximum  at  n x,  as  shown 

in  Fig.  5.3(a),  the  most  favorable  amplitude  of  tlie  rf  power 

is  that  for  which  x = n,  in  which  case  J has  a maximum 

n 

'f/2n  = — -T  (haj^^/2eV^^)  The  increased  dissipation  in  the 

weak  link  raises  the  effective  noise  temperature  = 

. All  these  requirements  give 


(5.31) 

Since  P 1,5  for  a given  value  of  R„,  P , T and 

dc  N N o c 

X _ . Eq,  (5.31)  gives  the  best  value  of  V for  these  para- 
rf  ^ max 

meters.  While  in  Eq.  (5,31)  one  could  use  tlie  value  of 
I^^(P  =0)  predicted  by  theory,  experimental  values  are 
usually  higher  and  it  is  better  to  use  the  actual  experimen- 
tal value.  Note  the  weak  frequency  dependence  of  Eq.  (5.31), 

Experimentally  we  find  the  same  value  of  V (T)  for  10  and 

max 

32  GHz,  except  in  a very  narrow  range  near  T , 

c 

One  can  use  Eq.  (5.31)  to  calculate  experimental  values 

of  P at  each  temperature,  from  the  measured  V values, 
o max 
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Tills  is  shown  in  Fig,  5„4  for  one  of  our  good  performance 
microbridges  in  which  at  low  temperatures  is  expected  to 
be  determined  exclusively  by  the  radial  contribution  to  the 
decrease  of  the  critical  current.  The  overall  agreement 
between  experimental  data  and  the  temperature  dependence  of 
the  heating  model  fitted  at  the  lowest  temperature  is  excel- 
lent. Very  near  T^,  the  experimentally  determined 
deviates  below  the  model  in  a consistent  way.  This  behaviour 
is  observed  in  all  of  our  high  performance  tin  microbridges 
and  we  attribute  it  to  the  bridge  contribution  Eq.  (5.26) 
becoming  more  important  near  T^  as  predicted  by  the  model. 
While  the  calculation  for  is  only  approximate  and 

represents  an  upper  limit  and  the  deviations  from  the  model 
are  small,  it  is  nevertheless  interesting  to  use  the  experi- 
mental data  to  attempt  to  obtain  an  estimate  of  the  length 
required  for  yield  the  observed  results.  If  the 

bridge  contribution  becomes  important,  rather  than  measuring 

P from  Eq.  (5.31)  we  measure  the  decrease  of  the  critical 
-P/P  ff 

current  as  e where 


P 


o,  ef  f 


P 

o 

1+  bridge 
radial 


(5.32) 


From  the  experimental  data  P /P 

o'  o,eff 


, /W  , 

bridge  radial 


0,5.  A typical  V 


max 


1.5  near  T 
at  T/T 


c 

0. 


which  yields  t^  = 1,43  and  an  effective  length  L 


so 

97  is  1 mV 
0,25  qm. 


32 


Tc 


P'ip;.  5.4  Temperature  dependence  of  P^(T)  for  a hitrh  per- 
formance microbridoe  (bridge  //16). 
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Tlius,  the  collinear  cone  approximation  is  indeed  a good 
approximation  as  this  length  is  smaller  than  that  measured 
from  SEM  pictures  of  our  microbridges.  One  should  be  care- 
ful however  as  this  is  only  an  approximate  value  for  L 
because  of  the  approximations  involved  in  the  calculation. 

Tliis  length  is  a factor  of  3 smaller  than  that  obtained  in 
Chapter  Four  for  a good  fit  to  the  effective  resistance  of 
the  foot  region.  This  discrepancy  is  of  no  concern  as  the 
nature  of  the  two  lengths  is  quite  different.  The  length 
measured  in  Chapter  Four  is  the  length  required  for  the  order 
parameter  in  the  bridge  to  reach  the  value  of  the  order  para- 
meter in  the  banks.  In  the  present  case,  the  effective 
length  measures  the  contribution  from  the  length  of  the 
bridge  which  is  not  contained  in  the  collinear  cone  approxi- 
mation. From  the  SEM  pictures  shown  in  Chapter  Two,  it  is 
clear  that  the  collinear  cone  geometry  would  approximate  at 
least  part  of  the  bridge. 

One  important  conclusion  of  the  above  analysis  is  that 
for  high  T^  materials  the  contributions  from  the  bridge  will 
be  even  more  important  than  in  tin  for  comparable  size  micro- 
bridges. In  Table  5.3,  the  best  values  of  V and  P (T  = 0) 

max  o 

obtained  from  Eq.  (5.25),  extrapolating  to  T = 0,  are  shown 
for  variable-thickness  microbridges  made  of  differeni  materi- 
als. Ihese  values  can  be  considered  to  be  in  good  agreement 
with  the  model.  The  values  for  the  lead  microbridge  are  below 
that  of  tin;  however,  the  zero  temperature  coherence  length 
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Table  5.3 

Parameters  For  The  Dost  Performance  Obtained  in 
Variable-Thickness  Microbridges  of  Different  Materials 


Material 

V„  (mV) 
Max 

P (T=0)(uW) 
o 

Tin 

0.54 

3.7 

8.1 

Lead 

0.28 

3.0 

5.9 

Niobium 

0.94 

0.45 

1.0 

Aluminum 

0.92 

0.97 

0.62 
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F (0)  is  almost  a factor  of  two  smaller  in  lead  than  in  tin 

so  that  the  contributions  from  the  finite  length  of  the 

bridge  are  more  important.  The  niobium  result  corresponds 

to  a fairly  long  {~  1 pm)  var iable-tliickness  bridge  made  of 
. 78 

very  dirty  niobium  so  that  both  the  thermal  conductivity 

K and  5(0)  are  very  small.  Tlie  low  values  of  P obtained  by 

o 

79 

Klapwijk  and  Mooij  in  aluminum  variable-thickness  micro- 
bridges are  harder  to  interpret.  The  values  are  con- 
siderably below  those  predicted  by  the  model,  but  no  system- 
atic search  for  the  largest  V was  made.  The  value  given 

max 

in  the  table  corresponds  to  a bridge  with  = 0.85  ohms 
at  a bath  temperature  of  1„07  K (T^  = 1.16).  At  the  micro- 
wave  power  at  which  the  measurement  was  made  the  width  of 
the  step  was  5 pA,  which  is  4%  of  the  zero  dissipation  criti- 
cal current.  The  frequency  used  was  35  GHz  so  V = 0.965 

max 

mV  corresponds  to  the  13^^  step.  The  Bessel  function 
has  a maximum  of  order  0.3,  for  which  the  step  width  observed 
is  16%  of  the  maximum  width  expected  in  the  absence  of  dis- 
sipation. "niis  is  a sizable  step  and  it  suggests  that  more 
steps  could  have  been  observed  with  this  sample.  It  should 
be  noted  that  their  tin  microbridges  of  comparable  size  to 
their  aluminum  microbridges  showed  steps  to  relatively  high 
voltages  although  they  had  values  of  below  those  obtained 

in  this  work.  While  these  explanations  seem  plausible,  with- 
out further  experimental  work  it  can  not  be  confirmed. 

The  results  of  the  model  can  also  be  used  to  estimate 
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V for  the  Nb  point  contacts  of  McDonald  ot  al.  For  a 

n.ax 

typical  resistance  of  ~ 10  ohms  is  estimated  to  be 

of  the  order  of  16  mV,  compared  to  17  mV  observed  experi- 
mentally. This  agreement  is  quite  good,  the  superior  per- 
formance of  the  point  contact  arising  from  its  smaller  size 
and  hence  larger  resistance.  Furthermore,  the  niobium  points 
were  made  of  bulk  niobium  wires  with  long  moan  free  paths  and 
thus  good  thermal  conductivity  compared  to  the  dirty  niobium 
films . 

5.2„5  Heating  Effects  at  Low  Voltages 

While  we  have  been  primarily  interested  in  the  high 

voltage  lim.its  for  the  observation  of  ac  Josephson  steps, 

heating  effects  can  also  significantly  modify  the  ac  steps  at 

lovv'  voltages  near  the  critical  temperature.  This  is  simply 

a consequence  of  the  small  power  levels  required  near  T to 

c 

raise  the  temperature  of  the  bridge  above  its  critical  tem- 
perature. In  terms  of  the  heating  model  it  can  be  said  to 

be  a consequence  of  the  small  value  of  P near  T „ 

o c 

Near  T^  quantitative  comparisons  with  the  heating  model 
become  rather  difficult  because  of  the  uncertainty  of  the 
bridge  contribution  to  tlie  decrease  of  the  critical  current. 
Furthermore,  one  has  to  assume  a particular  dependence  of  the 
widths  of  the  ac  steps  on  power  in  the  absence  of  heating 
effects.  In  Fig.  5,5  we  plot  tlie  maximum  full-width  of  the 
first  ac  step  normalized  by  the  zero  power  critical 

current  I^^  against  the  normalized  frequency  Q.  The  data 
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FifJ.  5.5  Experimental  values  of  the  maximum  width  of  the 
first  Josenhson  steo  as  a function  of  the  nor- 
malized frequency.  The  solid  line  corresnonds 
to  the  predictions  of  the  current-biased  RSJ 
model . 
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wore  taken  on  a tin  and  a lead  microbridge  at  10  GHz  and  on 
a second  tin  microbridge  at  both  10  and  32  GHz.  The  solid 
lino  corresponds  to  the  value  of  predicted  by  the 

current-biased  RSJ  model.  Wliile  the  values  for  Q « 1 for 
bridges  L-2  and  21  agree  well  with  this  model,  as  Q increases 
witli  decreasing  near  T the  observed  step  widths  are  below 
the  tlieoretical  prediction.  In  the  case  of  bridge  19  the 
data  at  32  GHz  agrees  well  up  to  0 ~ 1 but  tlie  values  at 
10  GHz  are  considerably  above  the  theoretical  curve.  Itiese 
variations  are  typical  of  our  variable-thickness  microbridges 
and  illustrate  the  difficulty  of  quantitative  comparisons 
with  any  predictions  of  the  heating  model  when  the  behaviour 
with  no  dissipation  is  not  precisely  known. 

In  Fig.  5.6  we  show  the  effects  of  heating  on  the  maxi- 
mum width  of  the  first  ac  step  if  the  behaviour  in  the  absence 
of  dissipation  is  given  by  the  current-biased  RSJ  model.  The 
curve  a = 0 corresponds  to  no  dissipation  where 

a = 1.5  V^c/Rj^P^(T  =0).  Since  corresponds  to  the  volt- 

age at  which  the  first  step  is  observed;  these  effects  are 
strongly  dependent  on  frequency.  The  trends  in  Fig.  5.6  are 
qualitatively  similar  to  those  obtained  in  tlie  experiment. 
Besides  the  difficulties  noted  previously,  quantitative  com- 
parisons are  difficult  because  the  actual  dissipation  is 
not  necessarily  equal  to  extent  which  depends  on 

the  normalized  frequency  Q. 
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Fipc.  5.6  Effects  of  heating  on  the  maximum  width  of  the  first 
Josephson  stgp  for  different  values  of  a =P/P  , 
assuming  P=V‘^/Rj,j  at  the  step.  ° 
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Behaviour  similar  to  that  shown  in  Fig.  5.6  has  been 

rei^ortod  as  a function  of  frequency  in  cross-scratched  indium 

80 

microbridges  by  Clark  and  Lindelof.  They  attributed  this 
behaviour  to  the  intrinsic  relaxation  time  of  the  order 
parameter  limiting  the  frequency  of  the  Josephson  oscillation 
in  tlie  microbridge.  Using  a phenomenological  relaxation  time 
for  indium  t = 10  (1  - T/T^)  ^ seconds  they  generated  the 

widths  of  tlic  ac  steps  by  incorporating  this  relaxation  time 
into  the  TDGL  theory.  Tlie  curves  they  obtained  in  this  man- 
ner are  qualitatively  similar  to  those  in  Fig.  5.6  for  fi  > 1, 
and  agreed  well  with  experiment.  However,  the  use  of  this 
time  has  no  justification  and  it  is  clear  that  heating  can 
play  a significant  role  since  their  microbridges  had  small 
P^'s  even  at  low  temperatures. 

Wliile  it  is  expected  that  an  intrinsic  relaxation  time 
should  play  a role  in  limiting  the  frequency  of  the  Josephson 
oscillation,  none  of  our  experimental  results  suggest  this 
limit  has  been  reached.  Wiile  the  exact  manner  in  which  the 
intrinsic  relaxation  time  will  impose  tliis  limit  is  not  under- 
stood, wo  have  observed  stops  near  T at  frequencies  v = 10^^ 

c 

Hz  at  a temperature  where  the  Ginzburg-Landau  time  corresponds 
to  ^ = 2 X 10'^^  Hz. 

The  effects  of  self-heating  can  manifest  themselves  in 
yet  another  way  by  cutting  off  the  oscillatory  Bossel- 
function  dependence  (on  the  rf  power  level)  of  the  amplitude 


Effects  of  heatin^  on  the  second  maximum  of  the  width  of  the  first  ac 
step  as  a function  of  frequency  for  different  values  of  P R,(mV)^. 
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of  a given  step.  As  the  rf  power  is  increased  above  the  value 
which  gives  the  first  nuiximum  for  each  step,  the  step  widths 
at  subsequent  maxima  are  substantially  reduced.  Our  variable- 
thickness microbridges  show  periodicity  with  increasing  power 
at  both  10  GHz  (as  shown  in  Fig.  5.2)  and  at  32  GHz.  At 
high  powers  the  steps  and  their  periodicity  become  harder  to 

observe.  In  two-dimensional  microbridges  Gregers -Hansen  et 
3 24  45 

al.  and  Jillie  ' have  reported  the  absence  of  periodi- 
city at  frequencies  above  10  GHz.  In  Fig.  5.7  we  show  the 

max 

reduction  of  the  second  peak  for  the  first  step  1^^  (P)/ 

max2  2 

I^  (P=0,  versus  frequency  for  various  values  of  P^R^(mV)  . 

Here  we  have  assumed  that  the  normalized  frequench  Q ^ 1 in 

order  to  use  the  simple  voltage-bias  approximation.  For  our 

high  performance  microbridges  at  low  temperatures  P R ~ 2. 0-4.0 
2 max^  ° 

(mV)  so  that  the  value  of  1^^  will  be  reduced  signifi- 
cantly only  at  frequencies  of  the  order  of  500-1000  GHz. 

2 

In  poorer  performance  bridges  near  T^,  ~ 0.02  (mV) 

so  that  at  the  frequencies  used  in  this  work,  the  periodicity 

would  be  harder  to  observe.  As  an  example,  a good  perform- 

2 n\ax_ 

ance  bridge  near  T has  P (T)R„  ~ 0.2  (mV)'^  and  I,  would 

c o N — i 

be  reduced  to  1%  of  its  zero  dissipation  value  at  a frequency 

of  the  order  of  100  GHz.  At  this  same  temperature  a poor 

performance  bridge  [V  (T=0)  =1.0  mVj,  would  have 
^ ^ max 

max_ 

Ij^  reduced  to  1%  at  38  GHz.  This  qualitatively  confirms 
our  belief  that  for  uniform-thickness  microbridges  the 
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absence  of  the  periodicity  of  the  ac  steps  at  higher  fre- 
quencies is  due  to  self-heating  effects,  since  we  have 
assumed  P = thus  the  cut-off  frequency  obtained 

is  only  an  upper  limit. 

5.2.6  Optimum  Design  of  Thin  Film  Microbridges  for  High- 
Frequency  Applications 

Ihe  remarkable  performance  improvements  at  high  fre- 
quencies observed  in  our  tin  and  lead  variable-thickness 
microbridges  do  not  imply  that  they  form  the  best  possible 
configuration  for  high  frequency  applications.  In  order  to 
design  and  fabricate  a Josephson  device  useful  at  high  fre- 
quencies one  has  to  optimize  the  different  parameters  of  the 
bridge  to  satisfy  the  following  criteria:  (a)  heating  effects 
should  be  minimized  in  order  to  extend  the  usefulness  of  the 
device  to  the  highest  possible  frequency;  (b)  the  impedance 
of  the  device  should  be  as  close  to  that  of  free  space  as 
possible  in  order  to  provide  efficient  coupling  to  the  electro 
magnetic  radiation  one  wishes  to  detect;  (c)  the  prod- 

uct should  be  maximized,  since  above  the  frequency  = 

the  resistive  shunting  degrades  the  performance 
of  the  device;  (d)  the  capacitance  C should  be  small  since 
the  RC  time  constant  defines  a frequency  above  which  this 
capacitive  shunting  dominates,  resulting  in  poorer  performance 

The  limitations  due  to  capacitance  are  of  little  signi- 

_2 

ficance  in  microbridges  since  C ~ 10  pF.  On  the  other 
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hr. nd,  one  requires  an  impedance  as  close  to  that  of  free 
space  as  possible.  While  tiiis  criterion  suggests  using 
materials  with  high  which  typically  have  high  resistivi- 
ties and  high  products,  the  heating  model  indicates 

that  these  materials  are  inferior  because  of  their  short  co- 
herence lengtli  and  low  termal  conductivity.  Iherefore,  the 
optimum  materials  appear  to  be  In,  Sn,  Pb  and  Nb.  We  have 
discarded  A1  because  of  its  low  and  product.  While 

niobium  is  still  impractical  because  of  its  poor  quality  in 
thin  films  it  is  appealing  in  terms  of  its  higher  resistivity. 
:he  optimum  solution  to  this  problem  should  be  the  fabrica- 
tion of  a variable-thickness  bridge  in  which  the  bridge  and 
tlie  banks  are  made  of  different  materials.  The  bridge 
could  be  made  of  a high  resistivity  material  such  as  niobium, 
and  the  banks  of  a high-conductivity  material  such  as  tin. 
While  the  niobium  bridge  will  have  a very  small  5(0),  the 
proximity  to  the  banks,  the  high  resistance,  and  the  good 
cooling  provided  by  the  banks  should  outweigh  this  disad- 
vantage. 

Hie  optimum  design  described  above  requires  sophistica- 
ted fabrication  techniques,  especially  in  guaranteeing  inti- 
mate electrical  and  thermal  contact  of  the  bridge  with  the 
bank  materials.  Two  separate  depostion  processes  would  be  re- 


quired. However,  we  have  tested  the  validity  of  the  concept  with 
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design  shown  in  Fig,  5,8.  This  is  a bridge  consisting  of  a 
very  dirty  0.13  pm  thick  tin  film  on  top  of  which  a clean, 
thick  (1.5  pm)  tin  film  is  evaporated.  This  is  done  sequen- 
tially in  the  same  deposition  to  insure  good  contact  between 
the  two  films.  In  this  manner  when  the  film  is  cut  with  a 
diamond  knife  the  bridge  is  made  out  of  the  very  dirty 
material  and  thus  should  have  a high  while  the  thick 

clean  film  should  provide  good  cooling.  Initial  tests  of 
t}iis  concept  have  yielded  a 2.5  ohm  microbridge,  totally  non 
hysteretic  even  at  1.5  K,  which  showed  ac  steps  up  to  4.2 
mV.  Ihis  performance  is  better  than  any  of  our  single- 
material  tin  or  lead  microbridges  and  tends  to  confirm  the 
validity  of  our  arguments  based  on  the  heating  model, 

5.3  Further  Evidence  of  The  Well-Cooled  Nature  of  The 

Variable  Thickness  Geometry 

In  the  preceding  sections  we  have  made  numerous  referen 
ces  to  the  well-cooled  nature  of  the  variable-thickness 
microbridges  based  on  their  high-frequency  performance. 
Additional  evidence  of  this,  some  of  it  strikingly  quanti- 
tative, can  be  obtained  from  the  dc  I-V  characteristics  of 
our  variable-thickness  microbridges. 

5.3.1  Hysteresis 

4 

Skoepol  et  al,  have  presented  strong  arguments  for  the 
connection  between  hysteresis  and  heating  in  both  long  and 
short  microbridges.  They  argued  that  since  no  dissipation 
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FiR.  5.8  Schematic  of  a bilaver  microbridfre  for  improved  hiph-frequencv  oerformance. 
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occurs  at  zero  voltage,  the  microbridge  would  reach  its  full 

critical  current  I . Wlien  this  current  is  exceeded,  dis- 
co 

sipation  occurs  and  if  tlie  bridge  is  not  well  cooled  a normal 

spot  would  develop.  The  current  required  to  sustain  this 

hot  spot  would  in  many  cases  be  smaller  than  I and  thermal 

hysteresis  would  occur  especially  at  low  temperatures.  This 

explanation  is  directly  confirmed  by  our  experiments.  In 

Fig.  5.9(a)  the  I-V  characteristics  are  shown  for  one  of  our 

best  performance  microbridges  with  V =3.7  mV.  Even  at 

max 

2.1  K the  curves  are  smooth  and  no  hysteresis  or  multi- 

valuedness  occur.  In  contrast  the  I-V  characteristics  of  a 

somewhat  less  optimal  microbridge  are  shown  in  Fig.  5.9(b). 

This  bridge  has  a V =2.3  mV  and  P =2.6  uW.  At  the 

max  o 

lowest  temperature  a multivalued  curve  begins  to  develop  and 
instabilities  are  observed  at  voltages  at  which  the  sub- 
harmonic gap  structure  (SGS)  is  observed.  In  microbridges 
showing  even  poorer  high-frequency  performance,  this  hyster- 
esis becomes  more  pronounced  and  is  first  observed  at  tem- 
peratures much  closer  to  T^.  Moreover,  little  or  no  change 

in  the  I-V  characteristics  or  V occurs  when  the  bath  tern- 

max 

perature  passes  through  the  lambda  point.  Tliis  indicates 
that,  as  was  assumed  in  the  heating  model,  the  predominant 
temperature  rise  is  associated  with  the  thermal  resistances 
of  the  metal  rather  than  that  of  surface  heat  transfer  to 


the  bath 
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5.3.2  Tlie  Subharmonic  Gap  Structure 

Further  quantitative  and  qualitative  evidence  of  the 

amount  of  heating  present  in  our  microbridges  is  provided  by 

the  gap  and  subharmonic  gap  structure  (SGS)  observed  on  the  I-V 

characteristics  at  the  gap  voltage  2A/e  and  submultiples  of 

it.  Tfiis  structure  has  previously  been  observed  in  thin-film 

microbridges  and  various  explanations  for  its  origin  have 

— fi  1 

been  proposed  and  extensively  debated.  ~ 

Tlie  reduced  hysteresis  of  our  microbridges  allows  us  to 
observe  this  structure  over  the  whole  bath  temperature  range 
in  our  well-cooled  microbridges.  In  Fig.  5.9(a)  the  struc- 
ture is  quite  subtle,  especially  at  the  subharmonics  of  the 
gap  = 2a (T)/ne  where  n > 1.  In  Fig.  5.9(b),  this  struc- 
ture is  quite  pronounced  both  for  the  gap  2A(T)/e  and  its 
submultiples.  Typically  the  less  well-cooled  the  bridge,  the 
sharper  the  observed  structure  becomes.  When  we  plot  -dl/dV 
against  voltage  as  shown  in  Fig„  5.10  (bridge  #4)  this 

structure  can  clearly  be  seen  as  peaks.  We  identify  the 

til 

characteristic  voltage  for  the  n peak  as  the  sharpest 
feature  in  the  I-V  characteristics,  which  corresponds  to 
the  inflection  point  to  the  right  of  each  peak  in  tlie  dl/dV 
measurement.  The  results  would  not  vary  significantly  if 
any  other  part  of  the  peak  were  chosen  as  the  index,  as  long 
as  the  microbridge  was  fairly  well-cooled. 


150 


As  shown  in  Fig.  5.10,  as  the  temperature  is  lowered 
the  peaks  become  sharper  and  more  peaks  are  observable. 
Typically  from  four  to  seven  of  these  peaks  are  observed, 
depending  on  the  microbridge.  In  not  so  well-cooled  micro- 
bridges as  in  Fig.  5.9(b)  instabilities  and  jumps  indicating 
the  presence  of  negative  resistance  regions  appear  first  in 
the  gap  peak  (n  = 1)  and  later  in  the  higher  order  peaks. 
Finally,  at  very  low  temperatures  the  whole  region  becomes 
hysteretic. 

When  the  gap  values  are  obtained  from  the  SGS,  v/e  ob- 
[ serve  that  the  higher  order  peaks  (n  > 2)  seem  to  correspond 

I 

to  approximately  the  correct  value  of  the  gap  from  the  BCS 

theory,  but  the  n = 1 and  n = 2 peaks  can  be  significantly 

depressed  in  voltage.  These  differences  become  even  more 

pronounced  at  low  temperatures.  In  Fig.  5.11(a)  and  (b) 

we  have  plotted  the  inferred  gap  voltages  nV  versus  the 

n 

voltage  for  two  microbridges.  If  these  structures  cor- 

respond to  the  same  value  of  the  gap  they  should  yield  a 
horizontal  line  nV^  = 2A/e.  As  shown  in  Fig.  5.11(c)  the 
nV^  values  for  large  n extrapolated  to  zero  dissipation 
(zero  voltage)  agree  well  with  the  theoretical  BCS  gap  func- 
tion for  tin  evaluated  at  that  temperature.  The  depressed 
values  for  large  dissipation  seem  to  correspond  to  definite 
rises  of  the  temperature  above  that  of  the  bath  due  to  the 
increased  dissipation  at  higher  voltages.  This  is  shown  in 
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Fip;.  5.10  Derivative  curves  showini?  the  variation  of  the 

subharmonic  p;ap  structure  with  temperature  ( ) . 


A 
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Fig.  5.9(a)  where  Uie  gap  value  at  the  lowest  temperature 
is  below  the  value  expected  from  the  BCS  theory,  but  above 
that  observed  at  the  same  temperature  in  the  less  well- 
cooled  bridge  shown  in  Fig.  5.9(b).  Thiis  is  further  con- 
firmed in  Tbble  5.4  where  we  show  the  experimental  value  of 

V = 2A/e  at  2.1  K as  well  as  the  value  of  V for  a sel- 

1 max 

ection  of  microbridges.  The  correlation  between  high  V 

^ max 

and  high  is  very  good.  Wo  have  found  that  by  observing 
the  value  of  the  gap  at  low  temperatures  and  comparing  it  to 
other  microbridgos,  we  can  predict  an  approximate  value  of 

V from  the  I-V  characteristics  even  before  any  microwave 
radiation  is  applied  to  it. 

It  is  possible  to  discuss  the  observed  temperature  rises 
in  a more  quantitative  fashion.  In  Fig.  5.12  we  plot  the 
values  of  the  temperature  rise  inferred  from  the  observed  gap 
for  the  n = 1 and  n = 2 peaks  for  one  of  our  microbridges. 

To  check  whether  the  inferred  temperature  rises  depend  on  the 
power  dissipation  in  a reasonable  way,  we  have  adopted  a 
simple  model.  The  relevant  gap  value  is  taken  to  be  deter- 
mined in  a region  at  a distance  of  the  order  of  the  coherence 
length  from  the  bridge  center  since  the  chemical  potential  of 


tlie  pairs  , should  approach  a constant  outside  a region  of 
P 

size  ~ P within  which  Up  is  time-dependent.  Assuming  a 
three-dimensional  radial  heat  flow  which  carries  away  the 


total  power  dissipated,  we  integrate  from  » to  § and  obtain 


POWER  (;xW) 

Inferred  temoerature  rise  from  the  rap  and  n=2  subharmonic  of  brid(?e  ^4 
as  a function  of  the  power  dissioated. 
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INFERRED  TEMPERATURE  RISE  (K) 
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Table  5.4 


Values 

of  V and  V,  at  2.1  K for  a 

max  1 

of  Microbridges 

Selection 

Bridge  # 

V (mV) 

max 

(mV) 

23 

4.2 

0.98 

17 

> 1.5* 

0.90 

11 

3.7" 

0.86 

16 

3.7 

0.95 

10 

2.6 

0.84 

9 

2.3 

0.84 

18 

2.0 

0.73 

3 

0.98 

0.73 

7 

1.7 

+ 

0.55 

4 

1.5 

0.56 

5 

1.4^ 

0.54 

6 

0.7 

0.50 

14 

< 1.0 

0.50 

13 

< 0.8 

0.50 

+ Limit  had  not  been  reached  at  maximum  microwave  power 
available 


* Bridge  was  used  for  far-infrared  experiments  described 

in  Chapter  Six,  and  V could  not  accurately  be  measured. 

max 


* 
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the  temperature  rise 


'^local“  "^bath  ~ 2f:K? 


(5.33) 


where  K is  the  thermal  conductivity  and  f*  is  the  solid 
angle  for  heat  flovv/  in  each  of  trie  banks.  In  Fig.  5.13  we 
plot  the  temperature  rise  from  the  SGS  but  now  against  P/' 
where  we  have  used  the  coherence  length  at  the  bath  tempera- 
ture. The  data  now  includes  the  data  from  the  n = 1,2  peaks 
of  the  bridge  in  Fig.  5.12  as  well  as  the  n = 1 peak  of  our 
best  cooled  microbridge.  The  data  for  both  microbridges 
fall  on  the  same  straight  line  as  expected  from  Eq.  (5.33). 

If  we  use  r.  = 2'*/3  and  the  observed  slope,  we  obtain  a very 
reasonable  value  for  the  therm.al  conductivity  of  tin  of  K = 
0.065  W/cm-K. 

Another  indication  of  the  well-cooled  nature  of  the 

microbridges  is  obtained  when  they  are  irradiated  by  micro- 

81 

waves.  In  contrast  to  other  cross-scratched  bridges,  the 
development  of  the  ac  steps  does  not  cause  the  SGS  to  fade 
away  until  the  power  level  is  increased  sufficiently  that  ac 
steps  are  observed  at  voltages  well  above  the  SGS.  When  tiie 
structure  starts  to  disappear,  it  first  decreases  in  strength 
and  later  shifts  to  lower  voltages  as  the  microwaves  heat 


the  bridge. 

Ihe  effects  of  nonequilibrium  dissipation  make  it 
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u« 


the  power  dissipated. 
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Jifficult  to  interpret  the  shape  of  the  gap  and  subharmonic 
gap  structures.  The  various  explanations  proposed  for  the 
origin  of  these  structures  are  based  on  either  single- 
particle tunneling  processes,  multiparticle  tunneling  pro- 
cesses or  self-coupling  of  the  Josephson  oscillation  to  the 
bridge.  These  models  yield  results  which  are  hard  to  dif- 
ferentiate in  experiments,  except  for  possible  correlations 

84 

with  other  properties  of  the  weak  links.  In  our  micro- 
bridges the  SGS  is  usually  a decrease  in  the  current,  which 
tends  to  support  the  explanation  based  on  the  self-coupling 
of  the  Josephson  radiation  in  the  bridge.  In  this  model, 
whenever  the  voltage  exceeds  V = 2i/me,  where  m is  an  even 
integer,  additional  pair  breaking  by  absorption  of  the 
Josephson  radiation  occurs  in  the  bridge.  This  would  reduce 
the  total  supercurrent  flowing  through  the  microbridge  as 
observed  in  our  experiment.  However,  this  can  only  explain 
the  even  series  of  the  structure,  unless  one  invokes  pro- 
cesses by  which  the  pair-breaking  is  accompanied  by  quasi- 
particle tunneling  across  the  bridge.  This  makes  it  diffi- 
cult to  understand  the  near  equality  of  the  intensities  of 
the  odd  and  even  subharmonics  observed  in  experiments. 

In  Fig.  5.14  the  gap  structure  is  shown  in  detail  for 
seven  different  microbridges  at  approximately  the  same 
reduced  temperature  T/T^  = 0.885.  (Because  of  the  difference 
in  T^  and  variations  in  the  temperature  at  which  the  curves 
were  taken,  T/T^  values  range  from  T/T^  = 0.882  to  0.891, 
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which  would  result  in  only  small  changes  in  the  expected 

value  of  the  gap.)  Tlie  current  scale  in  Fig.  5.14  has  been 

scaled  by  the  resistances  of  the  bridges  so  that  all  the 

curves  have  approximately  the  same  slope,  and  the  curves 

plotted  in  order  of  increasing  V^.  Tlie  shape  of  the  I-V 

characteristic  varies  significantly  from  the  bridges  with  a 

low  V , to  those  with  the  largest  V,  . Wliile  accurate  V 

1 1 max 

data  is  not  available  for  all  the  bridges  shown  in  Fig.  5.14, 

where  it  is  available  the  trend  in  the  shape  correlates  with 

V in  the  same  way  as  witli  V,  . In  the  case  of  those 
max  ^ 1 

bridges  with  the  highest  such  as  #11  and  #17,  V was 

1 max 

very  high  for  one  bridge  and  was  not  studied  in  detail  for 
the  second  (#17).  From  other  indications  such  as  the 
hysteresis,  we  can  determine  that  these  were  very  well- 
cooled  bridges.  In  these  well-cooled  bridges,  the  gap  struc- 
ture is  a small  depression  of  the  current  at  relative  to 
the  trend.  In  the  not  so  well-cooled  bridges,  as  determined 
from  lower  V values,  lower  V,  values  and  large  amounts 
of  hysteresis  at  low  temxxjra turos , the  shape  of  the  gap 
structure  is  quite  different.  In  these  bridges  there  is  a 
broad  region  of  negative  resistance  accompanied  by  a large 
shift  to  lower  currents  above  the  characteristic  voltage 
As  noted  earlier,  tliis  sJiift  to  lower  currents  would  be 
expected  if  the  self-coupling  of  the  Josephson  oscillation 
in  the  bridge  loads  to  a largo  amount  of  pair-breaking  and 
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t_iiOj.of ore  a reduction  of  the  supercurrent  flowing  through 
Uie  bridge. 

We  believe  that  the  voltage  width  of  the  gap  structure 
ooserved  in  our  microbridges  can  be  attributed  to  the  tem- 
perature rise  producing  a gradient  of  the  gap  in  the  micro- 
bridge  and  p^irt  of  the  banks.  Thus,  the  onset  of  the  pair 
braking  wilJ  not  bo  sharply  defined  at  = 2A(T^)/e,  cor- 
responding to  the  l^atli  temperature  but  rather  spread  over  a 
range  of  gap  values  determined  by  the  temperature  rise  near 
the  bridge.  Tfie  actual  origin  of  the  negative  resistance 
regions  is  more  difficult  to  explain.  Tlie  current  jumps  sud- 
denly to  a stable  lower  value  at  higher  voltages.  One  could 
envision  some  sort  of  instability  created  by  tlie  nonoqui librium 
nature  of  the  quasi-particles  created  at  the  gap  edge,  which 
will  not  readily  diffuse  out  of  the  region  of  depressed  gap 
in  a poorly-cooled  microbridge.  Tfie  occupation  numbers  of 
the  quasi-particles  near  the  gap  will  be  strongly  modified, 
leading  to  changes  in  the  gap  and  thus  the  total  super- 
current through  the  bridge.  Any  further  interpretation  of 
these  features  in  poorly-cooled  bridges  is  complicated  by 
the  temperature  distribution  present  as  well  as  by  the 
"feedback  effect"  that  a change  in  the  gap  will  have  in  the 
pair -breaking  process. 

Tlie  overall  trend  observed  in  the  gap  structure  as  the 

34 

cooling  is  improved  is  similar  to  that  reported  by  Yanson  in 
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microbridges  formed  as  shorts  through  dielectric  layers. 

In  Yanson's  bridges,  however,  a change  is  observed  from  a 

reduction  of  the  total  current  above  V^,  to  an  increase  of 

tine  total  current  above  V^.  Tlie  latter  v/as  only  observed 

in  microbridges  with  normal  state  resistances  higher  than 

those  of  our  microbridges.  This  increase  of  the  current 

84 

at  the  gap  is  typically  observed  in  high-resistance  well- 
cooled  point  contacts,  and  suggests  a common  origin  for  the 
SGS  in  these  two  types  of  weak  links  when  they  are  very 
well-cooled o 


CHAPTER  SIX 


OTHER  TOPICS 

» 

6.1  Detection  of  Far-Infrared  Radiation  with  Variable- 

Thickness  Microbridqes 

Tlie  observation  of  microwave-induced  Josephson  steps 
up  to  voltages  corresponding  to  far-infrared  frequencies  in 
our  variable-thickness  microbridges  encouraged  us  to  attempt 
the  observation  of  an  a step  induced  by  far-infrared  (FIR) 

I radiation,  Hiis  unsuccessful  attempt  is  described  in  this 

i 

' section  and  the  experimental  difficulties  encountered  are 

discussed . 

An  optically-pumped  dielectric-waveguide  FIR  laser 
built  in  our  laboratory  by  D.  Weitz  was  used  in  these 
experiments.  We  used  the  496  um  line  of  CH^F  gas  with  a 
maximum  output  power  of  ~ 10  mW.  The  radiation  was  focused 
on  the  microbridge  through  two  crystal  quartz  windows  in  the 
side  of  the  helium  dewar. 

Ibe  samples  used  in  these  experiments  were  similar  to 
, those  used  in  our  earlier  work,  except  tliat  a thin-film  half- 

wave dipole  antenna  which  was  cut  with  a diamond  knife  next 
to  the  bridge  as  shown  schematically  in  the  inset  of  Fig. 
6.1(a).  The  purpose  of  this  antenna  was  to  improve  the 
coupling  of  the  radiation  to  the  microbridge.  Tlie  actual 
length  of  the  antenna  was  150-250  qm  depending  on  the  sample 
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anu  run.  The  length  of  tlie  antenna  was  purposely  made 
shorter  than  half  the  wavelength  of  the  radiation  in  order 
to  compensate  for  the  reduction  factors  due  to  the  finite 
width  of  the  antenna  arms,  the  dielectric  constant  of  the 
substrate,  and  the  thickness  of  the  sapphire  substrate.^^' 

At  4.2  K,  the  absorption  of  the  far-infrared  radiation  by 
tlie  sapphire  sliould  be  small  and  no  significant  heating  of 
the  substrate  should  occur.  The  sample  was  mounted  so  that 
the  antenna  arms  could  be  rotated  about  an  axis  perpendicu- 
lar to  the  polarization  and  the  direction  of  the  incident 
radiation. 

A total  of  six  samples  with  resistances  ranging  from 
0.15-1.25  ohms  were  tested.  In  none  of  these  samples  was 
the  strength  of  the  coupling  sufficient  to  detect  the  effects 
of  the  FIR  directly  on  tlie  I-V  characteristics  or  its  deriva- 
tive dl/dV.  No  ac  Josephson  steps  were  observed  at  1.25  mV. 

A small  response  of  the  sample  to  the  radiation  could  be 
observed  by  chopping  the  laser  beam  and  detecting  the  cur- 
rent change,  I , through  the  bridge  by  measuring  the  vol- 
tage  across  the  0.1  ohm  resistor  in  series  with  the  bridge 
at  the  chopping  frequency  with  a PAR  124  lock-in  amplifier. 

In  Fig.  6,1  (a)  and  (b)  we  show  the  I-V  characteristic  at  3.4  K 

and  we  show  the  observed  response  I as  a function  of  vol- 

ac 

tage  for  tlie  best  sample  used.  This  sample  had  a normal 
state  resistance  of  1.25  ohm  and  therefore  was  essentially 
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voltage-biased  by  the  biasing  circuit  and  current-biased 

by  the  radiation.  The  sample  was  completely  nonhysteretic 

at  2.1  K with  a gap  structure  at  = 0.90  mV.  From  the 

discussions  of  Chapter  Five,  one  would  expect  a V of  at 

least  3.0  mV  for  these  characteristics.  Microwave  steps 

were  observed  up  to  1.5  mV  by  irradiating  the  sample  with  32 

GHz  radiation  with  the  Gunn  diode  placed  three  inches  away 

from  the  crystal  quartz  windows  in  the  side  of  the  dewar. 

This  did  not  constitute  an  upper  limit  as  not  enough  micro- 

waves  could  be  coupled  with  this  crude  arrangement  to 

determine  V . The  FIR  response  shown  in  Fig.  6.1(b)  is 
max 

the  maximum  response  observed  for  any  orientation  of  the 
sample.  The  response  was  found  to  be  quite  sensitive  to  the 
orientation,  suggesting  some  form  of  antenna  response.  How- 
ever, the  small  response  observed  made  it  difficult  to 
accurately  determine  the  antenna  pattern.  If  one  assumes 
that  the  depression  of  the  critical  current  in  the  presence 
of  the  radiation  is  given  by  the  zero  order  Bessel  function, 

then  for  small  signal  currents  this  reduction  I can  be 

o 

written  as 
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rf 


(6.1) 


where  V^^  is  the  rf  voltage.  Then  assuming  the  I-V  charac- 
teristic is  described  by  the  RSJ  model,  the  current  change 

I will  be  proportional  to 
Si  c 
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(6.2) 


as  long  as  I > While  the  observed  response  is  quali- 

tatively similar  to  Eq.  (6.2)  at  least  at  low  voltages,  this 
does  not  necessarily  rule  out  the  possibility  that  the 
observed  response  is  due  to  heating  in  the  metal. 

From  Eq.  (5.24),  the  change  in  in  the  presence  of  small 

dissipation,  also  has  the  form  of  Eq.  (6.1).  In  order  to 
separate  these  two  effects,  one  would  need  a more  extensive 
knowledge  of  the  coupling  or  the  direct  observation  of  an  ac 
step. 

If  the  observed  response  corresponded  to  a true  Joseph- 
son  response,  it  would  not  be  surprising  that  no  ac  step  was 
observed  directly  on  the  I-V  characteristics.  Thermal  noise 

rounding  would  prevent  the  observation  of  a step  until  I^  = 

—8 

Ico^n  ^^®^r f^'^^rf  ^ ~ our  experiment  Ij^  ^ 10~  amps 

_7 

in  the  best  case  and  ekT^/h  <=»  10  so  that  better  coupling 
would  have  been  required  in  order  to  observe  the  first  ac 
step  of  the  496  ijm  radiation. 

These  results  emphasize  the  disadvantage  of  the  low 
impedance  of  our  microbridges.  If  their  impedance  could  be 
increased  to  » 10  ohm  without  significantly  decreasing 
the  thermal  conductivity  of  the  banks  as  suggested  in  section 
5.2.6,  the  coupling  should  be  significantly  improved  and  an 
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ac  step  might  bo  observed » Further  improvements  could  be 
made  in  the  antenna  design  by  reducing  the  width  of  the 
antenna  or  by  attempting  to  determine  its  optimum  length. 
However,  determining  this  length  can  be  very  difficult  as 
it  would  require  a series  of  runs  in  which  the  antenna 
length  was  shortened  each  time  by  small  amounts  and  the 
same  microbridge  used  without  any  change  of  its  properties. 

6.2  Arrays  of  Microbridges  in  Series 

In  the  course  of  our  experiments,  we  became  interested 
in  exploring  the  feasibility  of  studying  arrays  of  micro- 
I bridges  made  with  our  fabrication  technique.  This  effort 

was  not  a comprehensive  study  of  such  arrays,  but  rather  a 
test  of  whether  we  could  produce  with  sufficient  reliability 
many  microbridges  in  series  with  properties  as  nearly  iden- 
tical as  possible.  Arrays  of  microbridges  are  quite  inter- 
esting from  the  point  of  view  of  applications  and  as  useful 
tools  for  the  study  of  the  nonequilibrium  properties  of 
superconductors.  From  a practical  point  of  view,  if  it  were 
' possible  to  fabricate  a coherent  array  of  n identical  micro- 

bridges  then  the  impedance  of  the  array  would  be  nR^^,  where 
R,,  is  the  normal  state  resistance  of  each  individual  micro- 

bridge.  Furthermore,  the  power  of  the  radiation  emitted  by 

2 

such  an  array  would  scale  as  n as  long  as  nR^  « Z^,  the 

8 6 

impedance  of  free  space.  This  super -radiant  property  of 


arrays  has  boon  observed  in  scries  arrays  of  tunnel 
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87  88 

junctions  as  well  as  proximity  effect  microbridges. 

Tlicse  arrays  can  also  be  useful  in  studying  nonequilibrium 
properties  of  superconductors  by  investigating  the  inter- 
actions betv/een  pairs  of  microbridges  which  are  independently 


biased  with  currents  flowing  in  either  the  same  or  opposite 


directions 


24,45,89-91 


Tlic  samples  used  in  our  experiments  were  fabricated 
using  the  same  techniques  outlined  in  Chapter  IWo.  Care  was 
taken  that  the  pressure  applied  in  making  the  cuts  with  the 
diamond  knife  was  similar  for  both  bridges.  Ihe  thickness 
of  the  banks  had  to  be  less  than  0.5-0. 7 qm  in  order  to 
fabricate  closely  spaced  pairs  of  microbridges.  Otherwise 
when  we  attempted  to  fabricate  microbridges  separated  by 
less  than  3 qm  from  center-to-center , making  the  second 
microbridge  would  push  the  thick  bank  material  onto  the  first 
bridge  and  it  would  be  severely  distorted  if  not  destroyed. 

All  of  our  pairs  of  microbridges  were  biased  in  series. 
For  separations  larger  than  6 qm  our  pairs  behaved  as  two 
independent  microbridges  with  two  distinct  critical  currents. 
It  was  possible  to  observe  separately  the  gap  and  subharmonic 
gap  structure  for  each  microbridge.  When  irradiated  with 
microwaves  with  low  amounts  of  power,  the  total  voltage 
across  the  two  bridges  showed  ac  steps  with  nonzero  differ- 
ential resistance,  corresponding  to  the  sum  of  a flat  step 
in  one  bridge,  and  the  differential  resistance  of  the  other. 
At  fairly  high  microwave  powers  these  ac  steps  would  overlap 
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in  current  and  a zero  differential  resistance  would  bo 

observed  in  this  region  of  overlap  at  twice  the  Josephson 

voltage  V for  a single  bridge.  This  so-called  "synchroni- 

zation"  has  been  observed  in  arrays  of  flat  microbr idges 

24  45 

separated  by  as  much  as  50  pni.  ' 

For  separations  smaller  than  6 pm,  a different  behaviour 
was  observed  both  in  the  absence  and  presence  of  microwave 
radiation.  A single  critical  current  I could  be  identified 

CO 

for  a wide  range  of  temperatures.  This  form  of  critical 

current  locking  was  observed  in  most  raicrobridges  separated 

by  less  than  6 pm.  Similar  results  have  previously  been 

reported  in  pairs  of  flat  indium  microbridges  by  Jillio  et 
24  89 

al.  ' They  interpreted  this  behaviour  as  being  due  to 
the  mutual  reinforcement  of  the  bridges,  when  they  are 
separated  by  a distance  smaller  than  the  quasi-particle  dif- 
fusion length.  In  our  samples  we  observe  this  critical 
current  locking  over  a much  wider  temperature  range  (down  to 
T/T^  of  0.3  or  less).  Tliis  was  possible  because  heating 
effects  in  our  samples  do  not  become  important  until  much 
lower  temperatures.  For  temperatures  below  T/T^  ~ 0.75 
a very  small  feature  suggestive  of  a second  critical  cur- 
rent was  observed  when  dl/dV  was  plotted  against  voltage. 

Ihis  second  feature  was,  however,  quite  different  than  that 
observed  when  the  first  bridge  began  to  develop  a voltage  and 
it  was  not  clear  whether  it  truly  corresponded  to  a second 
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critical  current.  Whenever  a single  critical  current  was 
observed,  a single  gap  structure  was  observed  at  twice  the 
gap  voltage  for  a single  bridge,  at  all  temperatures, 
suggesting  the  I-V  curves  of  the  two  bridges  were  indeed 
quite  similar  to  each  other. 

The  most  interesting  effect  observed  in  our  pairs  of 
microbridges  occurred  for  our  highest  impedance  pair  ~ 2 
ohm  where  the  bridges  were  separated  by  6 pm.  This  pair 
exhibited  critical  current  locking  and  when  irradiated  with 
microwaves  it  showed  the  behaviour  seen  in  Fig.  6.2.  At 
very  small  microwave  powers  a small  ac  step  with  nonzero 
differential  resistance  was  observed  at  2Vj  (~  132  pV) . 

Any  increase  of  the  power  yielded  zero  differential  resis- 
tance steps  at  2V  Note,  however,  the  large  width  of  the 

U 

step  at  2V  . If  one  compares  the  maximum  width  of  this  first 

U 

step  to  the  zero  power  critical  current,  a value 
~ 1.5  - 2.0  is  obtained,  which  is  twice  as  large  as  any 
observed  in  our  single  microbridges  as  described  in  section 
5.2.5  and  is  almost  twice  as  large  as  the  maximum  width  pre- 
dicted by  the  RSJ  model  for  any  bias  condition.  Note  also 
that  there  is  also  a large  Wyatt-Dayem  effect,  which  is  quite 
large  considering  the  high  resistance  of  the  microbridges. 

At  higher  voltages  the  next  2-3  steps  also  show  a zero  dif- 
ferential resistance  and  large  widths.  For  high  microwave 
powers  the  ac  steps  begin  exhibiting  a finite  slope  as  seen 


As  the  temperature  is 


17  3 

in  the  last  curve  of  Fig.  6,2. 
lowered  the  observed  behaviour  is  similar  to  that  shown  in 
this  figure  down  to  T/T^  ~ 0.74.  Below  this  temperature  a 
zero  differential  resistance  step  is  only  observed  at  high 
microwave  powers  when  the  individual  steps  are  large  enough 
to  partially  overlap.  Because  of  the  large  separation  of  the 
bridges  (~  6 pm),  this  suggests  the  interaction  occurs 
through  the  branch-mixing  length  At  these  temperatures 

the  width  of  the  ac  step  at  V becomes  comparable  to  that  of 

U 

the  step  at  2V  The  interpretation  of  these  results  is  quite 

U 

difficult  without  knowing  the  behaviour  of  the  ac  steps  for 
each  independent  microbridge.  From  our  experiments  on  single 
microbridges  it  appears  quite  unlikely  that  the  individual 
steps  would  be  so  large  for  each  microbridge.  It  could  be 
that  the  large  steps  are  due  to  a voltage-locking  of  both 
microbridges  whenever  one  of  the  microbridges  is  biased  to  a 
constant  voltage  step. 

Tliese  observations  in  arrays  of  microbridges  indicate 
that  our  fabrication  technique  could  be  quite  useful  in 
studying  the  interactions  of  microbridges.  Because  of  reduced 
heating,  the  variable-thickness  geometry  would  have  the 
advantage  of  providing  a wider  temperature  range  to  study 
the  interactions  of  the  bridges.  Obviously  in  order  to  fur- 
ther clarify  the  nacurr  of  the  effects  observed,  independent 
bias  of  the  microbridges  would  be  necessary. 


CHAPTER  SEVEN 


CONCLUSIONS 

When  our  experimental  investigation  of  the  properties 
of  variable-thickness  microbridges  began  it  was  clear  that 
heating  effects  severely  limited  the  high  frequency  per- 
formance of  two-dimensional  microbridges.  It  was  our  initial 
objective  to  fabricate  microbridges  in  which  heating  effects 
would  be  significantly  reduced  by  the  thicker  banks  so  that 
their  performance  at  high  frequencies  could  be  extended  to 
the  far-infrared  region  of  the  spectrum.  We  found  that 
besides  obtaining  the  improved  high  frequency  performance, 
our  variable-thickness  microbridges  were  quite  useful  in 
studying  the  properties  of  microbridges,  as  tlie  improved 
cooling  allowed  us  to  study  them  over  a wider  temperature 
range  than  had  previously  been  possible.  Furthermore  by 
fabricating  brid  es  with  varying  degrees  of  cooling,  it  was 
possible  to  study  and  explain  how  the  various  regimes  of  the 
I-V  characteristics  and  the  ac  response  of  the  microbridges 
were  affected  by  heating.  Our  study  of  the  wide  variety  of 
phencmiena  observed  in  our  microbridges  in  different  regimes 
has  helped  identify  and  explain  the  extent  to  which  their 
Josephson  behaviour  is  modified  by  nonequilibrium  dissipa- 
tion and  other  nonequilibrium  effects.  In  what  follows  we 
will  consider  regions  of  increasing  voltage  along  the  I-V 
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curves  of  the  bridges  summarizing  our  current  understanding 
of  the  phenomena,  and  identifying  areas  where  further  in- 
vestigation is  required. 

The  zero  voltage  critical  currents  of  both  our  tin  and 
lead  microbridges  were  found  to  have  a Josephson-like  tem- 
perature dependence,  as  expected  from  their  small  size  and 
tile  presence  of  the  thick  banks.  This  is  to  be  contrasted 
with  the  wide  range  of  behaviour  observed  in  two-dimensional 
microbridges.  The  magnitudes  of  the  products  of  our 

microbridges  at  all  temperatures  were  found  to  be  larger  than 

those  expected  for  weak  link  devices.  While  these  larger 

32 

values  are  consistent  with  recent  microscopic  calculations 
at  low  temperatures,  no  explanation  is  yet  available  for  the 
magnitudes  of  the  critical  currents  near  T^.  In  this  regime 
it  is  clear  that  further  theoretical  analysis  is  required 
for  the  case  of  the  neck  radius  of  the  weak  link  being  com- 
parable to  the  mean  free  path  or  tlie  coherence  length  of  the 
material.  Even  at  the  lowest  temperatures  the  onset  of 
resistance  at  I in  our  best-cooled  microbridges  was  found 
to  be  completely  nonhysteretic,  confirming  the  heating- 
related  origin  of  the  hysteresis  which  is  observed  in  most 
superconducting  microbridges . 

At  low  voltages  (0-50  qV)  all  of  our  tin  microbridges 
exhibited  a foot-like  feature  characterized  by  an  initial 
region  of  constant  differential  resistance  followed  by  a 
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rapid  growth  of  the  voltage  witli  little  current  change  at  a 
characteristic  current  I , . This  feature  became  more  pro- 
nounced  at  low  temperatures  unless  heating  effects  became 
significant.  Systematic  evidence  has  been  presented  which 
supports  the  idea  that  this  low  voltage  feature  is  caused  by 
an  enhancement  of  tlie  supercurrent  flowing  through  the  weak 
link.  This  enhancement  is  produced  by  the  nonequilibrium 
quasi-particles  not  being  able  to  relax  because  their  relax- 
ation time  is  slow  compared  to  the  Josephson  period.  Com- 
parisons of  the  experimental  data  with  the  predictions  of  the 
models  of  Golub  and  of  Aslamazov  and  Larkin  have  shown  many 
points  of  agreement  in  the  regimes  where  the  theories  are 
applicable.  Good  agreement  was  found  with  the  magnitude  and 
temperature  dependence  of  the  initial  differential  resistance 
observed,  the  maximum  enhanced  supercurrent  and  the  re- 

lation of  the  excess  supercurrent  I to  I ^ rather  than  I 

s cl  CO 

While  both  the  Golub  and  the  Aslamazov-Larkin  models 
agree  well  with  some  aspects  of  our  experimental  data,  further 
theoretical  work  will  be  required  in  order  to  achieve  a com- 
plete understanding  of  all  of  tlie  observed  phenomena.  The 
models  are  only  complementary,  being  applicable  in  different 
regimes.  Furthermore,  some  of  the  observed  behaviour,  such 
as  the  absence  of  an  initial  steep  rise  in  voltage  at  I 

CO 

and  the  saturation  of  at  low  temperatures,  are  not  pre- 

dicted by  the  models.  Golub's  model  has  been  used  in  the 
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j-ow  voltage  regime  where  approximate  solutions  to  the  non- 
equilibrium Ginzburg-Landau  equations  for  short  weak  links 
are  valid.  Numerical  solutions  of  this  equation  are  required 
to  extend  the  comparison  of  this  theory  to  higher  voltages. 

At  still  higher  voltages  the  assumptions  of  both  models 
would  break  down,  as  the  Josephson  period  Tj  becomes  compar- 
able to  the  diffusion  time  or  the  finite  relaxation  time 
of  the  order  parameter  role  of  these  relaxation 

times  in  determining  the  I-V  characteristics  remains  unclear 
at  this  time  except  in  the  context  of  time-dependent  Ginzburg- 
Landau  theory,  which  might  not  be  justified  in  the  regimes 
of  interest. 

Experimentally,  further  work  is  required  in  order  to 
understand  the  lov/  voltage  I-V  characteristics  of  bridges  of 
other  materials.  We  found  the  characteristics  of  our  lead 
microbridges  to  be  quite  different  at  low  voltages  tiian 
those  of  our  tin  bridges.  It  is  probable  that  these  dif- 
ferences arise  from  the  large  width  of  these  bridges  compared 
to  the  coherence  length,  which  would  allow  single  vortices 
to  flow  through  the  bridges,  making  flux-flow  effects  import- 
ant. Studios  of  other  materials  with  small  coherence  lengths 
might  help  clarify  this  point.  It  would  also  be  quite 
interesting  to  investigate  the  low  voltage  regime  of  aluminum 
microbridges.  For  bridge  lengths  larger  than  those  of  our 
microbridges  one  would  expect  to  observe  similar  enhancement 
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effects  although  at  much  lower  voltages  than  in  tino 

At  higher  voltages  (V  ^ 50  uV)  tlie  I-V  characteristics 
of  tin  microbridges  are  severely  affected  by  heating  effects 
and  modified  by  the  characteristic  features  of  the  gap  and 
subharmonic  gap  structures  at  = 2A(T)/ne  where  n = 1,2, 

...  . We  found  the  actual  voltages  at  which  this  structure 
was  observed  could  be  depressed  significantly  by  heating 
effects.  The  actual  shape  of  the  SGS  structure  was  found 
to  be  broadened  both  in  voltage  and  current  by  heating 
effects,  particularly  when  it  occurred  at  lower  voltages 
in  bridges  with  poorer  high  frequency  performance  and  more 
hysteresis  at  low  temperatures.  This  confirmed  the  correla- 
tion between  heating,  hysteresis  and  the  voltage  limit  for 
high  frequency  performance  in  superconducting  microbridges. 
Using  a simple  heating  model  we  have  confirmed  that  the  tem- 
perature rise  inferred  from  the  gap  and  subharmonic  gap 
structure  depends  on  the  dissipated  power  in  a reasonable 
way.  The  values  of  the  gap  obtained  from  the  zero  dissipa- 
tion extrapolation  were  found  to  be  in  very  good  agreement 
with  the  BCS  theory. 

When  our  variable-thickness  microbridges  were  irradiated 
with  microwaves,  nonequilibrium  effects  were  found  to  sig- 
nificantly modify  their  properties  at  all  voltages.  We 
observed  the  Wyatt-Dayem  effect  only  in  our  tin  miciohriJ-  > . 

and  enhancements  of  the  critical  current  of  as  n.uch  ' 


AO-A049  779 


UNCLASSIFIED 

30^3 

*S)49779 


harvard  UNIV  CAmBrIOM  MA  DIV  of  applied  sciences  F/t  9/s 

NONEOUILIBRIUM  properties  of  variable -thickness  SUPERCONDUCTINS— ETC(U) 
JAN  7t  N OCTAVIO  N00019-77-C-00SS 

TR-13  NL 


4 


179 


were  observed  in  our  lowest  impedance  microbridges.  No 
enhancement  of  the  critical  temperature  was  observed  in  any 
of  our  samples  at  either  10  or  32  GHz.  'Ihe  temperature 
dep>endence  of  the  critical  current  enhancements  differ  from 
those  observed  by  others  in  tin  and  indium  microbridges. 

Our  results  do  not  help  clarify  the  origin  of  the  V^att-Dayem 
effect  in  short  microbridges.  Further  understanding  of  this 
effect  would  require  theoretical  work  more  closely  related 
to  the  particular  experimental  situations  used  and  would 
require  more  extensive  experimental  study  of  the  enhancements 
[ in  variable-thickness  microbridges. 

In  the  presence  of  microwaves,  our  bridges  showed  the 
ac  Josephson  effect,  and  the  maximum  voltage  at  which  it 
could  be  observed  clearly  depended  on  the  extent  of  heating 
effects.  We  found  that  the  improved  cooling  of  the  variable- 
thickness geometry  allowed  us  to  observe  microwave-induced 
steps  up  to  4 mV,  a factor  of  4 improvement  over  any  micro- 
bridge of  any  material  reported  by  other  workers  to  date. 
However,  with  far-infrared  radiation  we  were  unable  to 
observe  a fundamental  n = 1 Josephson  step  at  these  voltages. 

» 

The  limitation  in  this  case  did  not  appear  to  be  heating,  but 
the  low  impedance  of  our  samples  which  did  not  allow  us  to 

i 

efficiently  couple  the  far-infrared  radiation,  Furtlier 
improvements  in  antenna  and  bridge  design  should  make  it  pos- 
sible to  study  the  response  of  variable-thickness  microbridges 
to  far-infrared  radiation. 
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Our  observations  at  high  voltages  led  to  a simple 
analysis  o5  the  effect  of  nonequilibrium  dissipation  on  the 
ac  properties  of  three-dimensional  weak  links.  From  this 
analysis,  we  found  that  the  critical  current  I falls  as 
e:<p(-P/P^),  where  P is  the  power  dissipated  in  the  junction 
and  P^  is  a material-  and  temperature -dependent  parameter 
of  order  10  qW  for  most  typical  materials  far  from  T^.  This 
decrease  of  the  critical  current  with  dissipation  was  found 
to  be  in  excellent  agreement  with  our  observations  of  the 
decrease  of  the  amplitudes  of  the  ac  Josephson  steps  at  high 
voltages.  This  reduction  of  together  with  the  increased 

noise  temperature  of  the  weak  link  in  the  presence  of  dissipa- 
tion leads  to  a maximum  voltage  at  which  ac  steps  can  be 
observed  for  given  junction  parameters,  Hie  values  for  this 
limit  predicted  by  the  model  were  found  to  be  in  very  good 
agreement  with  those  observed  experimentally  in  both  tin  and 
lead  microbridges.  Using  the  experimentally  determined  limits 
at  different  temperatures  we  have  confirmed  the  temperature 
dependence  of  the  characteristic  power  level  for  the  decrease 

of  the  critical  current  P (T) , Near  T , we  have  found  addi- 

o c 

tional  contributions  to  the  decrease  of  I arising  from  the 

CO 

finite  length  of  the  microbridge. 

At  low  voltages  and  near  T^,  the  heating  model  predicts 
that  the  ac  Josephson  steps  can  also  be  modified  strongly  by 
heating  effects.  Comparisons  of  these  predictions  with  our 
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experimental  data  showed  good  qualitative  agreement  with  the 
observed  reduction  of  the  maximum  height  of  the  first  ac  step 
as  well  as  the  observation  by  others  of  a cut-off  frequency 
of  the  oscillatory  Bessel-function  dependence  of  the  step 
widths.  In  our  well-cooled  microbridges,  the  cut-off  fre- 
quency is  comparable  to  those  used  in  this  work  only  at  tem- 
peratures very  close  to  T^;  at  very  low  temperatures  it 
should  be  of  the  order  of  1 THz. 

From  the  heating  model  we  were  able  to  obtain  conditions 
for  the  optimum  design  of  variable-thickness  microbridges 
for  high  frequency  applications.  We  found  that  a microbridge 
made  from  single  homogeneous  film  would  not  yield  any  further 
significant  improvement  over  our  best  microbridges  unless 
the  size  could  be  reduced  significantly.  However,  improve- 
ments should  be  observed  in  a bi-layer  configuration  in  which 
the  bridge  is  made  of  a high  resistance  material  and  the 
banks  of  a cle.m  material.  This  geometry  should  maintain  the 
well-cooled  nature  of  our  bridges  but  would  give  a smaller 
dissipation  at  a given  voltage  because  of  the  higher  resis- 
tance, Initial  experiments  using  an  approximation  to  this 
configuration  support  this  conclusion. 

Finally,  we  believe  that  experiments  using  arrays  of 
variable-thickness  microbridges  could  be  very  useful  in 
improving  our  understanding  of  the  nature  of  these  devices 
as  well  as  their  interactions.  Our  initial  experiments  with 


A 
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series-biased  bridges  have  shown  that  closely  spaced  bridges 
appear  to  have  significant  locking  effects  even  when  not 
independently  biased.  In  one  case  the  observed  Josephson 
steps  were  abnormally  large.  One  advantage  of  the  variable- 
thickness geometry  seems  to  be  that  the  interaction  of  the 
bridges  can  be  studied  over  a wider  temperature  range  without 
simple  heating  effects  overwhelming  more  subtle  nonequili- 
brium interactions  between  the  bridges. 


I 


r 

APPENDIX  A 

, CHARACTERISTIC  PROPERTIES  OF  TIN  AND  LEAD  MICROBRIDGES 

USED  IN  IHIS  REPORT 

« 

In  this  appendix  a sununary  of  the  characteristic  pro- 
perties of  all  tin  and  lead  microbridges  cited  in  the  text 
is  included  for  completeness.  In  many  cases,  a particular 
parameter  is  not  included  for  a microbridge,  either  because 
it  was  not  measured  or  it  was  not  known  with  sufficient 
accuracy. 
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Characteristic  Properties  of  Tin  Microbridges 


Bridge 

^ (a) 

dl 

R CO 

^ dT 
(mV/K) 

(mV) 

Exponent 

of 

CO 

Maximum 

Dayem 

AI  /I 

CO  CO 

1 

0.15 

0.78 

1.65 

1.50 

2 

0.08 

0.80 

3 

0.14 

0.71 

1.75 

1.09 

0.04 

4 

0.24 

1.37 

2.71 

1.07 

0.24 

5 

0.27 

0.96 

4.67 

- 

- 

6 

0.40 

0.46 

1.40 

- 

- 

7 

0.43 

1.02 

1.63 

1.03 

0.06 

8 

0.44 

0.75 

3.30 

1.09 

0.04 

9 

0.36 

0.56 

0.90 

1.07 

0.11 

10 

0.85 

1.27 

1.62 

0.87 

0.01 

11 

0.33 

0.56 

1.32 

0.97 

0.03 

12 

1.14 

1.05 

1.49 

] .05 

0.08 

13 

0.09 

0.68 

1.16 

0.87 

- 

14 

0.20 

1.16 

1.38 

1.00 

- 

15 

0.06 

0.53 

- 

1,10 

- 

16 

0.44 

0.80 

1.19 

0.96 

0.13 

17 

1.25 

0.69 

1.25 

- 

- 

18 

0.13 

0.77 

1.23 

- 

0.27 

19 

0.06 

0.76 

0.88 

0.94 

- 

20 

1.76 

0.42 

1.44 

0.91 

0.09 

21 

0.07 

0.77 

1.10 

1.25 

0.43 

22 

0.62 

0.43 

1.60 

1.03 

- 

23 

2.5  ^ 

1.0 

0.85 

— 

— 

a Microbridge  used  in  far-infrared  experiments 

b Indicates  that  not  enough  power  was  available  to 

determine  V - 

max 

c Bilayer  bridge  as  described  in  Section  5.2.6  ^ 


185 


Characteristic  Properties  of  Tin  Microbridges 


Bridge 

# 

^cl'^^co 

V 

max 

(mV) 

P 

o 

(^iW) 

Gap 
Value 
(low  T) 
(mV) 

Bank 

Thickness 

(lim) 

1 

0.30 

0.10 

0.1 

2 

0.63 

0.63 

1.0 

3 

1.96 

0.70 

0.98^ 

1.0 

0.73 

1.0 

4 

1.29 

0.71 

1.5  ^ 

1.5 

0.56 

3.0 

5 

c 

0.71 

1.4  ^ 

1.2 

0.54 

3.0 

6 

- 

- 

0.70 

0.18 

0.50 

0.1 

7 

1.28 

0.78 

1.7 

1.2 

0.55 

1.0 

8 

1.38 

0.70 

1.3^ 

0.75 

0.50 

1.0 

9 

1.70 

0.72 

2.3 

2.6 

0.84 

1.0 

10 

1.19 

0.65 

2.6 

1.7 

0.84 

1.0 

11 

1.84 

0.72 

3.7^ 

8.1 

0.86 

1.7 

12 

1.44 

0.77 

1.6 

0.5 

0.62 

1.0 

13 

1.67 

0.78 

- 

- 

0.50 

1.0 

14 

1.62 

0.79 

- 

-■ 

0.50 

1.5 

15 

2.00 

0.77 

- 

- 

- 

1.0 

16 

1.72 

0.75 

3.7 

6.3 

0.95 

1.5 

17 

1.59 

0.77  > 

1.5® 

- 

0.90 

1.5 

18 

1.62 

0.72 

2.0 

4.5 

0.73 

1.5 

19 

1.48 

0.75 

1.5 

7.2 

0.59 

1.5 

20 

1.28 

0.73 

2.2 

0.6 

0.56 

1.0 

21 

1.34 

0.74 

- 

- 

0.40 

1.0 

22 

1.20 

0.77 

- 

- 

- 

1.5 

23 

1.40 

0.71 

4.2 

2.2 

0.98 

1.5 
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Characteristic  Properties  of  Lead  Microbridges 


Bridge  # 

R 

N 

(fi) 

I (4.2K) 

CO 

(mA) 

V 

max 

(mV) 

P 

o 

(pW) 

Exponent 
of  I (T) 

CO 

L-1 

0.80 

1.87 

3.0 

2.4 

L-2 

0.16 

33.6 

1.9 

3.8 

1.06 

L-3 

0.60 

3.5 

1.8 

0.98 

1.00 

L-4 

0.08 

37.5 

1.8 

5.9 

- 

L-5 

0.08 

20.5 

1.8 

5.7 

- 

L-6 

0.07 

16.8 

1.6 

5.2 

- 

L-7 

10.0 

0.042 

No  steps 

- 

- 

L-8 

0.10 

4.1 

1.1 

2.7 

_ 
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